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On the Mechanism of the (n, 2n) Reaction 
at 14 MeV Neutron Energy - II. 


K. WINTER (*), B. ToRKI and E. REMY (+) 


Laboratoire de Physique Atomique et Moléculaire du Collège de France - Paris 


(ricevuto il 28 Agosto 1958) 


Summary. — The correlation of the respective directions of emission of 
the two neutrons in the *Be(n, 2n)*Be reaction has been investigated by 
using coincidence detection of the outgoing neutrons. In an azimuthal 
plane, separation angles ranging from 45° to 90° are predominant; a 
maximum correlation of (20+3)% is obtained at 60°. This angular cor- 
relation is attributed to simultaneous emission of two neutrons. Angular 
distribution measurements on the same reaction proved, by using time- 
of-flight techniques, that second neutrons are emitted by a compound 
nucleus, whereas first neutrons are emitted both directly and by a 
compound nucleus. 


1. — Introduction. 


Investigations of the mechanism of three-particle reactions, such as the 
(n, 2n) reaction, are important for the interpretation of anomalies in medium- 
energy nuclear reactions: the nature of the compound systems formed in the 
course of three-particle reactions and the manner of their disintegration can be 
studied separately; whereas, in a two-particle reaction, such as the inelastic 
scattering of neutrons, effects of these two stages could not be distinguished. 

In a previous paper (1) (referred to as I) we have reported on measurements 

(*) Now at CERN, SC-Group, Geneva 15. 

(+) Now at Max-Planck-Institut für Physik, München. 

(1) E. Remy and K. Winter: Nuovo Cimento 10, 664 (1958). 
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of the distribution in energy of neutrons emitted in (n, 2n) reactions from the 
interaction of 14 MeV neutrons with various elements. By using time-of-flight 
techniques, continua of the first neutrons were separated from low energy groups 
due to the subsequent emission of a second neutron in (n, 2n) reactions. It is 
concluded (see I) that the principal mechanism of the (n, 2n) reaction in 
medium weight and heavy nuclei is successive emission of the two neutrons. 
The compound systems emitting the second neutrons have lifetimes which are 
long compared with the transit time of nucleons over nuclear dimensions: 
i.e. there is a well-defined quantum state, a compound nucleus, between the 
emissions of the first and the second neutron. 

Fig. 1 shows a competition diagram for different mechanisms of the (n, 2n) 
reaction. A time sequence has been assumed in the formation of the first 
compound system by 14 MeV neutrons incident on the target nucleus. In the 


1.compound system 2.compound system residual nucleus 


neutron incident 
on target nucleus 


two-particle 
interaction 
stage 


compound 
nucleus 


stage 


lato dl 


first approximation (1), the incident wave is considered as a wave scattered 
by the real part of the complex potential of the target nucleus. In the second 
approximation (2a), the wave inside the potential is modified by the inter- 
action of the incident neutron with a single nucleon of the target nucleus, 
and finally by the interaction with others, until a compound nucleus is for- 
med (3a). However, the two-particle interaction stage (2a), having a life- 
time of the order of the transit time of nucleons over nuclear dimensions 
(10-21 s), can disintegrate, in competition with the final formation of a com- 
pound nucleus, by direct emission of a nucleon, e.g. a neutron, leading to a 
second compound system (2b). Again, there is a possibility of competition 
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between final formation of a compound nucleus (3b) and disintegration by 
direct emission, if energetically possible, of a second neutron to the residual 
nucleus (2¢). 

A sequence of the (n, 2n) reaction in the two-particle interaction stage, 
leading from (2a) to (2b) and (2e). (see Fig. 1), should be considered as simul- 
taneous emission of the two neutrons; whereas, once a compound. nucleus 
has been formed, either (3b) or (8a), the emission of the two neutrons will 
be successive, leading to a distribution in energy which shows two separate 
groups, due to the first and the second neutrons (see I). 

Evidently, there are two different possible courses of the (n, 2n) reaction 
leading to successive emission and to such a distribution in energy: either 
the first neutron is emitted directly by the two-particle interaction stage (2a) 
with subsequent formation of a compound nucleus (3b) via (2b); or it is emitted 
by the compound nucleus (3a). Both possibilities lead to successive emission 
of the two neutrons; they cannot be distinguished by a measurement of the 
energy spectrum. 

This paper will be concerned with the investigation of this alternative. 

As a direct test for the emission of neutrons by the two-particle inter- 
action stage (2a, 2b), we investigated the respective directions of emission 
of the two neutrons in the 14 MeV neutron induced *Be(n, 2n)&Be reaction. 
To decide the above mentioned alternative we investigated separately the 
angular distributions of the first and the second neutron. 


2. — Angular correlation measurements. 


Two neutrons, forming a virtual singlet state in the two-particle inter- 
action stage, may be emitted simultaneously (see Fig. 1). In the center of 
mass system of the virtual neutron-neutron state, the two neutrons will be 
emitted at 180° with respect to each other. In the laboratory system, the 
respective directions of emission will be correlated at angles smaller than or 
equal to 90°. Consequently, a contribution of this pick-up mechanism in (n, 2n) 
reactions can be identified by its angular correlation. 

ANSEROVA et al. (2) investigated the respective directions of emission of the 
neutrons emitted in *Be(n, 2n)’Be reactions, by using Be-loaded photographic 
emulsions (see I). After examination of 84 events they concluded that the 
two neutrons are not emitted in the same direction and with the same energy, 
but in different directions and with different energies. This investigation has 
been repeated recently (*), also with negative result. 


(2) N. S. Anserova, lou. A. NEMÎLOV and B. L. FoUNCHTEÎN: Rec. Trav. Inst. 
Rad. V. G. Chlopin, T, 114 (1956). 
(8) R. CHasrez and Cao-XUAN: private communication. 
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However, no experiments on the angular correlation of neutrons emitted 
in (n, 2n) reactions have been reported using coincidence detection of the out- 
going neutrons. 


21. Experimental method. — A conical beam of 14 MeV neutrons, yielded by 
the T(d, n)*He reaction, has been obtained by using a combined paraffin-lead 
collimator (Fig. 2). The ratio of the neutron fluxes inside and outside the 

collimated beam has been measured to be 


of the order of 10°. 
a Within this cone, at 108 cm of the neu- 
Paraffine | tron source, a beryllium scatterer, 5 cm long 
Ro Io and 3 em in diameter, was placed. It was 


| faced, in an azimuthal plane, i.e. perpendi- 


cular to the incident neutron beam, by two 


r— 40 neutron detectors containing plastic scintil- 
Ecru lators, 5 cm long and 5 cm in diameter. The 


two neutrons emitted in the *Be(n, 2n)*Be 
reaction have been detected in coincidence: 
signals supplied by the neutron detectors are 
fed into a fast-slow coincidence circuit connected with a time-to-pulse-height 
converter. The total time resolution was 5 mus. A complete description of 
the electronic circuits has been given elsewhere (4. Discriminators in the 
slow lines are adjusted in such a way that the slow coincidence output signal 
gates open the amplitude analyser following the time-sorter only for signals 
analysed in the fast circuit which satisfied the required cut-off pulse height. 
The discriminators thus established a minimum neutron energy accepted for 
coincidence detection; this minimum energy has been measured to be 500 keV 
(see I). 

In this way, a time spectrum was measured corresponding to the differences 
in flight times of the two neutrons emitted in the *Be(n, 2n)*Be reaction. The 
full width at half height of this difference time spectrum was 5 mys; the peak 
of the spectrum was shifted, by the help of delay lines, to the center of the 
amplitude analyser. The superposed, uniformly shaped, spectrum of random 
coincidences was simply subtracted. The contribution of (n-y) coincidences 
has been proved to be negligible by inserting absorbers between one of the 
neutron detectors and the beryllium scatterer. The quantity of y-rays ac- 
cepted for coincidence detection has been found to be of the order of 1%, 
by measuring the time-of-flight spectrum (see Fig. 6). The integral of the time- 
spectrum, after subtraction of random coincidences, supplied the number of 


Fig. 2. 


(4) E. Remy and K. Winrer: Journ. Phys. Rad., 18, 112A (1957). 
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coincidences for a given flux of incident neutrons. The neutron flux has been 
monitored by counting «particles associated with 14 MeV neutrons in the 
T(d, n)*He reaction. 


2°2. Results and discussion. — One of, the neutron detectors, kept in a fixed 
position, established the direction of emission of one of the two neutrons. 
By varying, with the 
help of the other de- 
tector, the respective 


directions of emission 
of the two neutrons, 
the number of coinci- 
dences for a given flux N DE 
of incident neutrons 
was measured at an- 
gles of separation of 
45°, 67°, 90°, 135° 1900 
and 180°. The num- 
ber of coincidences as 


T 


a function of the se- 


i ue lo) 7 ay LESS L 1 L 1 
parating angle of the 0 90 780° 


respective directions 
of emission is shown 
In Kol eS. 

In addition to 
the investigated effect 
of angular correla- 
tion, crystal-to-crystal 
scattering of single 
neutrons varies as a 
function of the angle. 
Some contribution of 
this effect on the 
coincidence rate was 
proved by the shape 
of the time spectrum 
which changed with 
the angle. The coin- 
cidence detection of 
erystal-to-crystal scat- 
tering of single neutrons requires a minimum energy loss of 500 keV in 
each of the scintillators. Since the energy loss by elastic neutron scattering 
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in plastic scintillators depends on the angle of scattering, this background 
angular correlation is a function of the distribution in energy of neutrons emit- 
ted in the ‘Be(n, 2n)*Be reaction. Therefore we were searching for an element 
with identical dis- 
tribution in energy 
of neutrons emitted 
under the bombard- 
ment of 14 MeV neu- 
trons and for which 
the (n, 2n) reaction 


N is energetically for- 
bidden at this neu- 
tron energy. These 
two conditions are 

200 fulfilled by, car- 
bon (5%). A contri- 
bution of (n-y) coin- 

0 0° 780° cidences from the 
FN inelastic scattering 

oe. on carbon has no 
angular correlation, 
since the y-rays are emitted by a compound nucleus. The time between the 
emissions of the neutron and of the subsequent y-ray is long compared with 
the transit time of nucleons over nuclear dimensions, 7.e. there will be no 
angular correlation in the respective directions of emission. Consequently, a 
comparison of the coincidence rates as a function of the separation angle, for 
beryllium and for carbon, provides a method to eliminate the background 
correlation of crystal-to-crystal scattering of single neutrons. The rate of back- 
ground coincidences as a function of the separation angle, obtained with a car- 
bon target is shown in Fig. 4. 

The difference of the coincidence rates obtained with beryllium and with 
carbon, after normalization at 180°, is attributed to the angular correlation 
of the two neutrons emitted in the *Be(n, 2n)*Be reaction. The difference 
coincidence rate, read from smooth curves, is shown in Fig. 5, indicating cor- 
relation of the respective directions of emission at separation angles between 
45° and 90°; a maximum correlation of (20 + 3)% is obtained at 60°. 


(5) J. GARG and B. Torkt: private communication. 
(9) For the neutron spectrum of *Be(n, 2n)8Be see I. 
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3. - Angular distributions of the first and the second neutron. 


Since subsequent emission of the second neutron in (n, 2n) reactions leads 
to separate groups in the energy spectra, due to the first and the second 
neutron, angular distributions of the first and the second neutron can be 
measured separately. Such a type of spectra has been obtained with medium 
weight and heavy nuclei (see I). Angular distribution measurements, reported 
recently (7), are indicating that both courses of the (n, 2n) reaction leading 
to subsequent emission of the second neutron (see I, Introduction) occur with 


medium weight and heavy nuclei. 


Be? (n2n) 


1007 


In the energy spectrum observed with *Be, under the bombardment with 
14 MeV neutrons, neutron groups at 1.3 MeV and at 2 MeV are resolved (see I), 
agreeing with calculated transition energies to levels of *Be. Neutron energies 
between 5 and 12 MeV (see Fig. 6) are not resolved, containing either a great 
number of successive transitions of first neutrons or the maximum intensity 
fer simultaneous emission or a mixture of both. In the spectrum observed there 


(7) B. ToRKI: Proc. of the Int. Congress on Nucl. Phys. (Paris, 1958). 
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is no indication of the transition to the 11.3 MeV level in *Be which would 
yield 2.7 MeV neutrons. We could not decide (see I) whether this level does 
not exist or whether this transition does not occur because the time between 
the emission of the first and the second neutron at this excitation energy is 
short, i.e. of the order of the transit time of nucleons over nuclear dimensions, 
so that there is no defined excited quantum state between the two emissions. 
A contribution of simultaneous emissions by the two-particle interaction stage 
(see Fig. 1) has been identified by its angular correlation. 

In addition, we investigated the angular distribution of neutrons emitted 
in the *Be(n, 2n)*Be reaction, to decide whether, in successive emissions, the 
first neutron is emitted directly by the two-particle interaction stage (24) or 
whether it is emitted by a compound nucleus (3a) (see Fig. 1). 

Direct emission by the two-particle interaction stage would give rise to 
anisotropic angular distributions, peaked at small angles (5); whereas, emis- 
sion by a compound nucleus would result in angular distributions which are 
symmetrical with respect to 90°, or even isotropic (*), if the level density of 
a highly excited compound nucleus would be proportional to (251), where j 
designes the level spin; but this additional hypothesis has not been proved 
to be valid in all cases. 


8 1. Experimental method. — The distribution in energy of neutrons emitted in 
14 MeV neutron induced *Be(n, 2n)*Be reactions has been measured at seven 
angles contained between 30° and 120°, by using time-of-flight techniques. 
Two measurements were made on the neutrons to establish the zero as well 
as the arrival time. 

A 100 keV deuteron beam bombarded a tritium-titan target. The reaction 
T(d, n)*He yields a 14 MeV neutron associated with a 3 MeV recoil «-particle 
emitted at 180° with respect to the direction of the neutron. Signals supplied 
by an «-particle detector established a time zero for associated 14 MeV neutrons 
in a conical beam, the width of which was defined by the solid angle of the 
a-particle detector. Within this cone, at 25 cm of the neutron source, a be- 
ryllium scatterer, 5 em long and 3 cm in diameter, was placed. It was faced 
by a neutron detector containing a plastic scintillator, 5 em long and 5 em 
in diameter, placed outside the time-correlated neutron beam at a distance 
of 30 em of the scatterer, this distance representing the flight distance before 
detection of neutrons emitted by the scatterer. The neutron signal thus sup- 
plied the arrival time. The time interval defined by the zero and the arrival 


(*) N. AUSTERN, 8. T. BurLER and H. McManus: Phys. Rev., 92, 350 (1953). 
(9) W. Hauspr and H. FESHBACH: Phys. Rev., 87, 366 (1952); L. WOLFENSTEIN: 
Phys. Rev., 82,:690 (1951). 
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time was measured by using a time-to-pulse-height converter. A complete 
description of the electronic circuits and of their performance has been given 
elsewhere (410), By inserting a combined lead-uranium shield between the 
neutron source and neutron detector, random coincidences were attenuated 
by a factor of 10, varying from 10% at 90° to 40% at 30° and to 35% at 
120° of the real coincidence rate. . 


3°2. Results and discussion. — Time-of-flight spectra have been measured at 7 
angles ranging from (30 +7)° to (120 + 7)°, each for an integrated flux of 
107 «-partieles or 10° neu- 
trons in 47. A neutron spec- 1 to 2 MeV 
trum obtained with *Be at N 


90° is shown in Fig. 6; the 900 Fase 
time scale is linear, the cor- i + 


responding neutron energies 

are indicated. 0 30 60 90 120° 
Angular distributions 

have been deduced by divi- 

ding the measured spectra 


into three regions, ranging from 1 to 2 MeV, from 2 to 4 MeV and from 4 to 
10 MeV, and by taking integrals of these regions as functions of the angle of 
emission. The first two regions, ranging from 1 to 4 MeV contain exclusively 

second neutrons emitted in the 


2 to 4 MeV ‘Be(n, 2n)*Be reaction (see I), 
N the latter contains first neu- 
500 trons; a small contribution 


condoni of subsequently emitted second 


neutrons: at 5.1 MeV and 
jar 30 60 90 —_ 720° 6.2 MeV can be neglected; but 
there is a contribution of neu- 
Fig. 8. trons emitted simultaneously. 
Angular distributions deduced 
for these three energy regions are shown in Figs. 7-9. 

Symmetry with respect to 90° has been found for neutron energies ranging 
from 1 to 2 MeV (Fig. 7) and from 2 to 4 MeV (Fig. 8), indicating the sub- 
sequent emission of second neutrons by a compound nucleus, in agreement 
with our previous results (see I). However, in the energy region from 4 to 
10 MeV (Fig. 9), there is a strong anisotropic contribution peaked at small 
angles as well as another contribution showing symmetry with respect to 90°. 


(9) E. Remy and K. Winter: Nuovo Cimento, 10, 664 (1958). 
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The latter one, showing symmetry at 90°, is attributed to first neutrons emitted 
by a compound nucleus ((3a) in Fig. 1), but the possible area under this con- 
tribution does not account for the total number of subsequently emitted second 
neutrons found in the energy region from 1 to 4 MeV. Consequently, first 
neutrons leading to subsequent emission of second neutrons must also be con- 
tained in the contribution which 
4 to 10 MeV is peaked at small angles, indi- 
cating direct emission of first neu- 
trons by the two-particle interac- 
1000 tion stage with subsequent for- 
= mation of a compound nucleus 

emitting the second neutron. 
The remaining portion of the 
anisotropic contribution, which is 


0 peaked atsmallangles, accountsfor 
(o) à 5 

30 60 30 120 neutrons emitted simultaneously 

Fig. 9. and with almost equal energies. 


4. — Conclusions. 


The principal mechanism of 14 MeV neutrons induced (n, 2n) reactions in 
medium weight and heavy nuclei is the successive emission of the two neutrons. 
The second compound system (see Fig. 1) has a lifetime which is long com- 
pared with the transit time of nucleons over nuclear dimensions: 2.6. a de- 
fined quantum state, a compound nucleus, is formed between the emissions of 
the first and the second neutron. This result, previously deduced from measu- 
rements of energy spectra (see I), has been confirmed by recent angular distri- 
bution measurements (7). 

On the other hand, the first neutron in a (n, 2n) reaction leading to sub- 
sequent emission of the second neutron is emitted both, directly by the two- 
particle interaction stage and by a compound nucleus (’). 

In 14 MeV neutron induced (n,2n) reactions with light nuclei, e.g. 
*Be(n, 2n)*Be, these alternative courses of successive neutron emission have 
been proved by the angular distribution measurement. 

A third possible mechanism, simultaneous emission of two neutrons in 
the *Be(n, 2n)*Be reaction, has been proved by the angular correlation measu- 
rements. This contribution is emitted mainly at small angles with respect to 
the direction of incident 14 MeV neutrons, the two neutrons having almost 
equal energies (see Figs. 6 and 9). 

A possible contribution of simultaneous emission of two neutrons in (n, 2n) 
reactions with medium weight and heavy nuclei cannot be proved since the 


= 
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(n-y) coincidence rate is very important. But there is some positive evidence 
for simultaneous emission in the 55Ni(p, 2p) reaction (11). 


We would like to express our gratitude to Professor FRANCIS PERRIN for 
his hospitality and for his continued and stimulating encouragement. 

We wish to acknowledge gratefully fellowships of the « Commissariat à 
l'Energie Atomique » (K.W.) of the C.N.R.S. (B.T.) and of the « Stifterverband 
für die Deutsche Wissenschaft » (E.R.). 


(4) B. lb. COHEN: Phys. Rev., 108, 768 (1957). 


RAA IS IDENSIRON (ES) 


Si è studiata la correlazione fra le rispettive direzioni di emissione dei due neutroni 
nella reazione °Be(n, 2n)*Be servendosi della rivelazione per coincidenza dei neutroni 
uscenti. Nel piano azimutale predominano angoli di separazione da 45° a 90°; si ha 
un massimo di correlazione di (20+3)% a 60°. Si attribuisce questa correlazione ango- 
lare all'emissione simultanea di due neutroni. Misure della distribuzione angolare nella 
stessa reazione hanno provato, con l’ausilio di tecniche di tempo di volo, che i secondi 
neutroni sono emessi da nuclei composti, mentre i primi neutroni sono emessi sia diret- 
tamente, sia da un nucleo composto. 


(*) Traduzione a cura della Redazione. 
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Zur Interpretation von Jets mit dem Modell 
der angeregten Nukleonen. 


J. BURMEISTER, K..LANIÙS und H. W. MEIER 


Deutsche Akademie der Wissenschaften zu Berlin 
Kernphysikalisches Institut Zeuthen - Berlin 


(ricevuto il 4 Settembre 1958) 


Summary. — A method for the determination of y, in nucleon-nucleon 
collisions of very high energy based on the model of excited nucleons 
in presented. Isotropic angular distribution of the secondary particles 
is the centre of mass systems of the emitting nucleons is assumed. The 
energy of all secondaries is assumed to be the same in both systems. 
Relations for y, and the Lorentz factor y’ of the nucleon after the colli- 
sions in the centre of gravity frame of reference are given. The y, values 
computed for a great number of high energy jets are compared with the 
values following from the statistical methed of Castagnoli et al. A com- 
parison is made with the data given by Koba for the multiplicity in the 
two centre model. The relations recently presented by TSCHERNAWSKT 
for the Heisenberg case of 7-N-collisions are in good agreement with 
our assumptions, 


Ejnleitung. 


Die experimentellen Ergebnisse über die Vielfacherzeugung von Sekundär- 
teilchen bei Wechselwirkung von hochenergetischen Teilchen mit Materie 
wurden in den meisten Untersuchungen mit den Voraussagen der Modelle von 
Heisenberg (*), Fermi (?) und Landau (*) verglichen. Allen drei Modellen ist 
gemeinsam, daf ein und nur ein Raumgebiet hoher Energiedichte entsteht, 
das dann die Ursache fiir die Vielfacherzeugung ist. 


(1) W. HEISENBERG: Zeits. Phys., 133, 65 (1952). 
(2) E. FERMI: Phys. Rev., 81, 683 (1951). 
() S. S. BeLeN'KIJ und L. D. LANDAU: Fortschr. Phys. Wiss., 56, 309 (1955). 
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In jüngster Zeit wurden insbesondere durch die Untersuchungen an Jets 
in Kernemulsionen von CIoK u.a. (+) Hinweise darauf gefunden, daB die 
Abstrahlung der Anregungsenergie unabhängig voneinander aus zwei Zentren 
erfolgt, die nach dem StoB im Schwerpunktsystem in entgegengesetzen Rich- 
tungen auseinanderfliegen (*). Theoretische Arbeiten, die auf dieser Modell- 
vorstellung beruhen, wurden schon vor einiger Zeit publiziert (TAKAGI (7), 
KRAUSHAAR und MARKS (8), KOBA (°)). 

In der vorliegenden Arbeit versuchen wir, eine gròBere Anzahl hochenerge- 
tischer Jets mit dem Zwei-Zentren-Modell zu interpretieren. Dabei nehmen 
wir an, dab die nach dem Sto8 im Schwerpunktsystem auseinanderfliegenden 
Zentren angeregte Nukleonen sind. 

Das in der Arbeit verwandte experimentelle Material wurde zum Teil in 
unserem Laboratorium gewonnen, zu einem anderen Teil entstammt es der 
Literatur (°°) und privaten Mitteilungen (+). 


1. — Berechnung von y, nach dem Zwei-Zentren-Modell. 


Den untenstehenden Beziehungen liegt im Schwerpunktsystem (S-System) 
das folgende Modell zugrunde (Abb. 1): Vor dem StoB bewegen sich beide 
Nukleonen mit gleicher Geschwindigkeit aufeinander zu. Infolge der Wechsel- 


(4) P. Crok, T. CocHEN, J. GIERULA, R. Hozynsxi, A. JURAK, M. Mimsowicz, 
T. SANTEWSKA, O. STANISZ und J. PERNEGR: Nuovo Cimento, 8, 166 (1958). 

(>) P. Crok, I CoGHEen, J. GimRuLa, R. HozyNskt A. Jurak, M. Mimsowicz; 
T. SANIEWSKA and J. PERNEGR: Nuovo Cimento (im Druck). 

(*) Zur gleichen Deutung kam auch Coccont (f) bei der Analyse einiger hochener- 
setischer Jets. 

(6) G. Cocconi: Vorabdruck. 
(eS RAA Cr Prog her AP SERIE 12381952) 
(8) W. L. KRAUSHAAR und L. J. MARKS: Phys. Rev, 93, 326 (1954). 

(9) Z. Kosa: Progr. Theor. Phys., 7, 123 (1952). 

(19) P. Crox, M. Danysz, J. GieruLA, A. JuRAK, M. Mimsowrcz, J. PERNEGR, 
J. VRANA und W. Worter: Nuovo Cimento, 6, 1409 (1957). 

(11) A. DEBENEDETTI, C. M. GARELLI, L. TALLONE und M. Vigone: Nuovo Cimento, 
4, 1142 (1956). 

(12) B. Epwarps, J. Losty, D. W. PeRKINS, K. Prinkav und J. REYNOLDS: Phil. 
Mag., 3, 237 (1958). 

(3) R. G. GLasser, D. M. Haskin und M. ScHEIN: Phys. Rev., 99, 1555 (1955). 

(4) E. G. Boos, A. CH. WINIZKI, J. S. TAKABAJEW und I. J. TSCHASNIKOW : Journ. 
Exp. Theor. Phys., 34, 622 (1958). 

(*) Wir mochten an dieser Stelle Herrn Prof. Dr. M. Murpsowicz, Krakau, 
Dr. J. GreruLA, Warschau, Dr. J.PERNEGR, Prag, sowie Dr. G. ZDANow, FIAN, Moskau, 
und ihren Mitarbeitern fiir die Uberlassung der Winkelverteilungen der Jets ihrer Labo- 
ratorien unseren Dank aussprechen. 
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wirkung findet eine Umwandlung von kinetischer Energie in Anregungsenergie 


der Nukleonen statt. 


Es wird angenommen, daf beide StoBpartner gleich 


stark angeregt werden und demzu- 


S-System L-System folge nach dem StoB gleiche Ge- 
1g) schwindigkeiten besitzen. 
IP) sE O Die Abstrahlung der Anregungs- 
; A % energie erfolgt ‘als Emission von 
vor dem Stop Sekundärteilchen und geschieht fur 
beide Nukleonen unabhängig von- 
| Pe, einander. Die Winkelverteilung der 
<2 | &— oe ae abgestrahlten Teilchen wird in den 
aa PET entsprechenden Ruhsystemen der 
nach dem Sto B strahlenden Nukleonen als isotrop 
Abb. 1. — Schematische Darstellung des vorausgesetzt. Die Energie pro Se- 
StoBes mit den entsprechenden Lorentz- kundarteilchen betrage in den Ruh- 
faktoren der Nukleonen (y = (1 — B?)-4). systemen E" — 0.5 GeV (*). 


Ausgehend vom Energiesatz im 
Laborsystem (L-System) erhalten wir unter Berücksichtigung der Anregungs- 
energien die Gleichung (*) 


1 18 5 } 
(1) Via V5 Vi Ve More (N1Y1 + N2y2)| - 


Hierbei ist m die Nukleonenmasse und n,j=mn,-4+mn, die Zahl der (von den 


beiden Nukleonen) beim StoB erzeugten geladenen Mesonen, y,, y, bedeuten 
die Lorentzfaktoren der Nukleonen nach dem StoB gegenüber dem L-System. 
Mit nani n,/2 folet 


1 


| 1 (3 BY \R 
(1a) Ver RI V3 (2 + V2) |! | —- n.) 


Fier € 
4 Mez 


Für den Lorentzfaktor der Nukleonen nach dem StoB im S-System liefert 
der Energiesatz im S-System die Beziehung 


3 y EE 
2 = ay Aes ae 
(2) Aa /( Lai 4 me? n) j 


(*) Siehe hierzu die Untersuchungen über die Transversalimpulse der Sekundär- 
teilchen von Jets (12:15), 


(5) O. MINnAKAWA, Y. NisHImurA, M. Tzuzuki;! H. !NAMANOUCHI, \H. Arzu, 
H. HaseGAWA, Y. YsHI, 8. ToEUNAGA; Y. Fusimoro, Si HAsEGAwWA, J. NISHIMURA, 


K. Niv, K. Yisnrgawa, K. Imanpa und M. Kazuno: INSJ-7 preprint, Tokyo, 
March (1958). 


(+) Die Indizes I und II beziehen sich auf das Ein- bzw. Zwei-Zentren-Modell. 
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Der Inelastizitatskoeffizient K, ergibt sich zu 


ed vou ! 
Le (me) (dn) 10 


II 


2. — Vergleich der Energiebestimmungen mit Hilfe der beiden Modelle. 


Da sich weder die Gleichung (1) noch (1a) direkt mit den Beziehungen 
von CASTAGNOLI u.a. (15) für da vergleichen lassen, wurde ein Vergleich der 
experimentellen Daten durchgeführt. Hierzu wurde Ve, in erster und zweiter 
Naherung (s. Anhang) sowie 1e nach der Gleichung (1) für alle verfiigbaren 
Jets mit Energien —10!! eV. und N, <5 berechnet. 

Zur Berechnung von y, wurden die Spuren auf die Kegel so aufgeteilt, 
daB bei einer deutlichen Trennung der beiden Kegel die entsprechenden Spuren- 
zahlen genommen wurden. War keine 
klare Unterscheidung môglich, so 
wurden die Spuren im Verhältnis 
1:1 getrennt. In einigen Fallen 
wurde über die Ergebnisse verschie- 
dener Aufteilungen gemittelt. 

Die Aussonderung der Nukleon- 
Kern-St6Be erfolgte näherungsweise 
dadurch, daB nur Ereignisse als 
N-N-StôBe augefaBt wurden, bei 
denen di mit der Beziehung 


Abb. 2. — Vergleich der nach dem Modell 
der angeregten Nukleonen berechneten 
(4) I 2ns — 1 Lorentzfaktoren y, : mit den Werten, die 
= == = I 
. 3 | . . 
2.4 y2 + 1 sich nach CASTAGNOLI u.a. (15) in erster 
Näherung ergeben; e) 1 <2; +) 1>2 


2 = re i 
berechnete Zahl der getroffenen Nu- (1: berechnet mit: Hilfe) von oe 


kleonen / kleiner als 2 war (*). 
In den Abb. 2 und 3 ist Va, in Abhangigkeit von Ve baw. oN dargestellt. 
Wie aus der Abb. 3 hervorgeht, stimmen die mit dem Modell der angeregten 


(16) C. CASTAGNOLI, G. CORTINI, C. FRANZINETTI, A. MANFREDINI und D. MORENO: 
Nuovo Cimento, 10, 1539 (1953). 

(*) Vorläufige Rechnungen hatten gezeigt, daf auch für das Modell der angeregten 
Nukleonen kleine 7-Werte auftreten, wenn der nach (4) berechnete Wert klein ist. 
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Nukleonen erhaltenen y,-Werte gut mit den Daten überein, die man bei 


4024 
7 


I 


[LE 


10} 


Berücksichtigung der Winkel- und 
Energieverteilung nach (!) erhält. 

Die Betrachtungen der folgenden 
Abschnitte wurden mit den nach 
Gleichung (1) berechneten Werten 
für y. durchgefùhrt. 


Abb. 3. — Vergleich von Vers mit den 


Werten y® in zweiter Naherung. 
I 


3. — Vergleich der Multiplizitaten mit den Aussagen der Theorien des Zwei- 


Zentren-Modells. 


Im folgenden seien die experimentellen Multiplizitäten mit den von KoBa (°) 
geforderten Teilchenzahlen verglichen. 


KoBA gibt für die Zahl der erzeugten Teilchen n, den folgenden Aus- 


druck an: 


ye 
ns 


~ | 
(5) Vs ~ } 
[ 


In der Abb. 4 ist lg n) in Ab- 
hängigkeit von lg Vo aufgetragen (*), 
wobei die einzelnen Gruppen durch 
verschiedene nach der Gleichung (3) 
berechnete ,,-Werte charakterisiert 
sind. 

Durch die MeBpunkte der ver- 
schiedenen Gruppen der Abb. 4 wur- 
den mit Hilfe der Ausgleichsrechnung 
Geraden Igy, =a,lg n, +b, gelegt, 
und die GròBen a,, b, auf die GròBen 
m(K,,), C(Kn) der Relation n° — 
=C(K,-)y{ 0° umgerechnet. Den Ver- 
lauf von m mit K,, gibt die Abb. 5 


(*) Hierbei ist n¥ = 


tur, Arad 


7 
[Rae irene Tish 


A 
log ne 
14-09 "8 
3 
13 ay 
+ 
10} + x ° 
x a 
x «& A A 
© 00 a 
0.8} MOIS a 
o 
0.6} log y. 
O Rei ar CURE DUT a 0 
Abb. 4. — Der Zusammenhang zwischen 


der Multiplizitàt und y, für verschiedene 
II 
Inelastizitàskoeffizienten; A) K,,= 0,76— 
+0,80; o) 0,81—0,82; x) 0,83—0,85; 
e) 0,86—0,88; +) 0,89--0,92. 


f n,— 1, falls das Primärteilchen ein n oder p ist, 
| ns — 2, falls das Primärteilchen ein x-Teilchen ist. 
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wieder. Wie hieraus hervorgeht, scheint n, für X,<0,8 von X, nahezu 
unabhangig zu sein, wahrend für wachsende Æ,-Werte ein Anstieg zu verzeich- 
nen ist. Weiterhin hat C(K,,) die A 

o 


gleiche GréBenordnung wie der Pro- MK) 
portionalitàtsfaktor der Relation (5). 06 

Kürzlich erhielten wir Kenntnis Hi 
von einer neuen theoretischen Arbeit aa 
von TSCHERNAWSKI (17) über die Viel- 02 Left 
facherzeugung von Mesonen. Hierin vi x 
behandelt er in dem Abschnitt über 07 i 08 E 09 A) 
die peripheren StòBe das Modell der Abb. 5.— Variation von m mit A,,; ©) der 
angeregten Nukleonen und untersucht von Koba (?) geforderte theoretische Wert 
verschiedene Versionen der Wechsel- von 0,75 für K,, = 1. 


wirkung. Beim Vergleich mit dem 

Experiment findet er, daB die von der Fermi-Landauschen Theorie der 7-N- 
Wechselwirkung geforderten Muliplizitäten gegenüber den anderen Versionen 
die beste Ubereinstimmung mit dem empirischen Material aufweisen. Die Win- 
kelverteilungen der erzeugten Teilchen in den Ruhsystemen der Nukleonen nach 
dem StoB sollten gemäl dieser Theorie stark anisotrop sein. Dies steht jedoch 
in krassem Widerspruch zu dem experimentellen Befund, wonach nahezu iso- 
trope Verteilungen vorliegen (5). 

Der durchgeführte Vergleich der Multiplizitaten scheint uns nicht zulässig 
zu sein, da sich die dabei verwendeten Werte für y’ auf ein Zwei-Zentren- 
Modell (4:5) beziehen, das sich wesentlich von dem Modell der angeregten Nu- 
kleonen unterscheidet. Wir vergleichen nunmehr die Beziehungen der Arbeit (17) 
mit den von uns verwendeten Relationen. 

Lést man die Gleichung (2) der vorliegenden Arbeit nach n, auf, so erhalt 


man (mit Ver = Ve) 


4 me? (ye ; 
5 — - 1 bi 
(6) nN; D JON È | ( ) 


Unter allen von TSCHERNAWSKI behandelten Versionen zeigt die für den 
Heisenbergschen Fall der 7-N-Wechselwirkung gültige Beziehung 


rid me? Ve 


A Nee er y 


4 
=) 


(7) D. TscuerNAwsKI: Vorabdruck. 

(*) Die Gleichung (6) gilt im Gegensatz zu (7) auch fir kleine Werte y,/y’. Im 
Grenzfall des elastischen StoBes ist y,/y' = 1 und somit n, = 0, d.h. es werden keine 
Sekundärteilchen erzeugt. 


=, 2 - Il Nuovo Cimento. 
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die beste numerische Übereinstim- 


453 mung mit (6), wobei gilt: 
20 
: , M Ye 
15F I (8) Ei mc" in ee A : 
Ma y 
I 
10 
Der wesentliche Unterschied besteht, 
sl darin, daB wir in unserer Beziehung 
%y) (die als Bestimmungsgleichung für y’ 
chy! È ie E 
is, os diente) H” = const = 0.5 GeV setzten. 
4 5 ; 
Die Abb. 6 zeigt den Verlauf der 
Abb. 6. — Vergleich der Relationen fiir die beiden Beziehungen (6) und (7) in 
Multiplizitäten; I) Gleichung (7) von dem Energiegebiet der von uns ver- 


Tschernawski (17), II) Gleichung (6), ver- 


i 3 à wendeten Jets. 
wendet in der vorliegenden Arbeit. 


In der Abb. 7 ist weiter der 
Zusammenhang zwischen dem Trans- 
versalimpuls und der GròBe y,/y' dargestellt. 

Es muf schlieBlich bemerkt werden, daf die Heisenbergsche Version der 
r-N-Wechselwirkung auch hinsichtlich der Winkelverteilung der erzeugten 
Teilchen in den Ruhsystemen der Nukleonen nach dem StoB Ubereinstimmung 
mit dem Experiment (5) zeigt, da 
diese Theorie eine nahezu isotrope 


= : a (GeV/c) 
Verteilung fordert (17). a I 

04 ME rei 

II 

02 
Abb. 7. — Vergleich der Relationen für > ee i sep! 
ug, i cor if 6 8 10 
die Transversalimpulse; I) Gleichung (8) | 
von Tschernawski (7), II) p, = const, verwendet in der vorliegenden Arbeit. 


4. — Die Inelastizitat. 


In der Abb. 8 sind die auf Grund des Modells der angeregten Nukleonen 
ermittelten Werte des Inelastizitàtskoeffizienten A,, (Gleichung (3)) in Abhän- 
gigkeit von y.,, aufgetragen. Die Abb. 9 zeigt K, als Funktion von n.. 

Wie aus beiden Abbildungen hervorgeht, liefert das Modell der angeregten 
Nukleonen relativ hohe K-Werte. Dies steht in Widerspruch zu den Ergeb- 
nissen anderer Untersuchungen, in denen niedrigere K-Werte ermittelt wurden. 
Bei der Beurteilung dieses Resultats muf man jedoch beachten, daB die Ver- 
wendung des Inelastizitatskoeffizienten zur Interpretation von Jets sehr pro- 
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blematisch ist. TSCHERNAWSKI kommt in der Arbeit (17?) zu der SchluBfol- 
gerung, daf sich sowohl vom experimentellen als auch theoretischen Stand- 
punkt aus keine eindeutigen Krite- 

rien zur Unterscheidung der Se- ms 


kundär- und Primärteilchen angeben K, 
lassen. r 
ost A 
10 
Ka È 
1 
0.8 x 
0.6} 
0.6 5 
dai 1 heu le à le i je ea ñ iF orden Uses 
6 8 10 2 4 6 80° 10 20 30 40 
Abb. 8. — Variation des Inelastizitäts- Abb. 9. — Variation des Inelastizitits- 
koeffizienten mit y... koeffizienten mit n,. 
* ok % 


AbschlieBend méchten wir Herrn Prof. Dr. M. Migsowicz, Dr. J. GIERULA, 
Dr. I. GRAMENITZKI, Dr. J. PERNEGR und Dr. G. ZpANov sowie Herrn 
T. COGHEN für die wertvollen Diskussionsbemerkungen auf der Tagung in Liblice 
bei Prag unseren Dank aussprechen. 


ANHANG 


Die Bestimmung von y7 (der obere Index weist auf die erste Näherung 


hin) erfolgte mit Hilfe der Beziehung ("*) 


n 


il Di 
lg ya = — ia > lg ted: . 
Rye! 


Bei einigen Jets wurden Spuren beobachtet, die innerhalb der MeBgenauig- 
keit in der Richtung des Primärteilchens lagen. Auf Grund der Symmetrie- 
bedingung im S-System wurde in diesen Fallen bei der Berechnung von y, 
zugleich mit dieser inneren Spur die Spur mit dem gròBten Winkel fort- 
gelassen. 

In zweiter Naherung gilt (1%) 


my gi In Va Are 
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Hierbei berücksichtigt A, die Abweichung der experimentellen von der iso- 
tropen Winkelverteilung im S-System und F die Form des Energiespektrums. 
Geht man von der Verteilung 


dF ~ cos?’ d cos 9 
aus, so gilt 


N,— 2n( < Ÿ) |1 —a? 
Ig | —_—_—_———" | = (2k + 1) lg} ———| =y ted. 
5 Ne ( bin VIS Tg Ver go 


Den so gefundenen Werten von k& wurden die in (1) angegebenen GréBen A, 
zugeordnet. 

Der Berechnung von F lag das Heisenberg-Spektrum zugrunde (!). Dabei 
wurde die « Abschneideenergie » für alle y.-Werte zu €, — 1.2 angenommen. 
Wie aus der Abb. 3 hervorgeht, laBt sich bei geeigneter Variation von a mit 


Ve, eine bessere Übereinstimmung mit der Geraden Tea erreichen. 


RIASSUNTO (*) 


Si presenta un metodo per la determinazione di y, nelle collisioni nucleone-nucleone 
di energia molto alta basato sul modello dei nucleoni eccitati. Si assume la distribuzione 
angolare isotropica delle particelle secondarie nel sistema del centro di massa dei nu- 
cleoni emittenti. L’energia di tutti i secondari si assume uguale nei due sistemi. Si danno 
relazioni per y, e per il fattore di Lorentz y' del nucleone dopo le collisioni nel sistema 
del centro di massa. Si confrontano i valori di y, calcolati per un gran numero di jet 
di alta energia coi valori ottenuti seguendo il metodo statistico di Castagnoli et al. 
Si fa un confronto coi dati di Koba per la molteplicità nel modello a due centri. Le 
relazioni recentemente presentate da TsCHERNAWSKI per il caso di Heisenberg delle 
collisioni 7-N sono in buon accordo con le nostre assunzioni. 


(*) Traduzione a cura della Redazione, 
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IL NUOVO CIMENTO Wioito SU ING a 1° Gennaio 1959 


Investigation of Nuclear Interactions of Energies 
between 10 and 100 GeV. 


E. FENYVES, E. Gompost and P. SURANYI 


Central Research Institute of Physics - Budapest 


(ricevuto il 9 Settembre 1958) 


Summary. — The angular and energy distribution of shower particles 
of seven nuclear interactions, mostly complex collisions, were measured. 
The Lorentz factor y.m of the center of mass system (CMS) was deter- 
mined from the energy of secondaries and the Lorentz factor ysm of the 
symmetry system was obtained from the angular distribution. It was 
found that Ysm/Yem = 0.92 + 0.15, which corresponds to a symmetrical 
emission of shower particles in the CMS. Energy distribution in the 
CMS as well as transversal momentum distribution of shower particles 
were also determined. 


1. — Introduction. 


We have measured the angular distribution of shower particles generated 
in nuclear interactions of cosmic ray particles in emulsion. The nuclear inter- 
actions were found by surface scanning in plates of the I stack irradiated in 
the Po-Valley Expedition, 1955. The only requirement in selecting the events 
was that the energy of almost all shower particles should be determinable 
without respect to the nature of the primary particle and the number of gray 
and black prongs. Thus the majority of these events corresponds to nucleon- 
nucleus or nucleus-nucleus collisions. Our aim was to investigate a) the sym- 
metry of the emission of shower particles in the forward and backward cone 
in the CMS, b) the energy distribution in the CMS, and c) the transversal mo- 
mentum distribution of shower particles. 
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2. — Symmetry of the angular distribution in the CMS. 


The velocity of the CMS relative to the laboratory system (LS) was deter- 
mined in the usual way (1) from 


> pi cos 0; 
DI 


= SE 7 TN) 
SE, 
i 


V= 


where p;, H; and Ÿ, are momentum, energy and angle of emission, respectively, 
of shower particles in the LS (*). From V the Lorentz factor of the CMS can 
be calculated: 

1 


Vem — V1 =: ye . 


As usual, all particles with I<1.5I, were considered as shower par- 
ticles and their energies and momenta determined at lower energies by simple 
scattering and at higher energies by relative scattering measurements. 

We have found seven nuclear interactions in which the energy of almost 
all shower particles could be measured. However, owing to the great statistical 
error in the values of the y,,, determined, the deviation from the forward- 
backward symmetry of the angular distribution indicated large fluctuations 
in individual cases. Therefore, the value of y, which best transforms the 
angular distribution as measured in the LS into a symmetrical one in the 
CMS was estimated for each event and the ratio y,,,,/y., Was built for each 
interaction (Table I). y,,, was determined according to the method of CASTA- 


TABLE I. 


Symbol of Jet Vem YVsym PERDE, 

| 3 eee — 
ace A 1 Baa GE Oss} OT 25 MI 
6 + lla Pte 16 | Ga) ae alee 1250 
5 + 10a Tele Dez | Oi) 46 BY 0.88 + 0.41 
MIL Ali QE Di | Tat te Ne AU = O57 

| gs Cra AS RS He) e iil 0.89 + 0.40 
eo 23,2 + 10.5 DEEE 4 il 1.09 + 0.61 


(*) The above relation holds only if we assume that the conservation of momentum 
is valid for charged particles independently of uncharged particles being emitted; 
this assumption can be made only if the number of shower particles is higb. 

(*) R. G. GLasser, D. M. HasxiN and M. ScHEIN: Phys. Rev., 99, 1555 (1955). 
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GNOLI et al. (?) from 


al 
ii == — = > Inte d, , 
$ i 


where n, is the number of shower particles. 
It can be seen from Table I that the weighted mean of the values of 
Verne: 18 
Vsym 
Vem 


(G2 ee O16 


and the fluctuation of individual values around this mean value is not signi- 
ficantly higher than normal. 

Thus we can conclude that the angular distribution of shower particles 
generated in such nuclear interactions which are mostly complex collisions 
and involving energies between 10 and 1000 GeV, is symmetrical in the forward 
and backward direction within the limits of error given above. This result 
is in good agreement with the hydrodynamical meson production theory of 
Landau (3) which predicts symmetrical emission of shower particles in the 
forward-backward direction in the CMS also for complex collisions. 

Applying the same procedure as above to the jet data published in the 
literature ('‘*) we obtain for these jets 


ESE Tis SEXO IL 


Vem 


which is in good agreement with our result. We must emphasize, however, 
that the data published in the literature are strongly biased, because the 
authors have expressly searched for jets of nucleon-nucleon collision character. 

Our result shows at the same time that symmetry of the angular distri- 
bution of the collisions in the CMS cannot be used as an argument proving the 
collisions to be nucelon-nucleon collisions, aS was done earlier. 


(2) C. CASTAGNOLI, G. CORTINI, G. FRANZINETTI, A. MANFREDINI and A. MORENO: 
Nuovo Cimento, 10, 1539 (1953). 

(3) S. Z. BeLEN'K1J and L. D. LANDAU: Usp. Piz. Nauk., 56, 309 (1955). 

(4) E. FRIEDLANDER: Studii di Cercetari de Fizika (in press). 

(5) E. G. Boos, A. H. Vinicx1y, J. S. TAKIBAEV and I. A. Casnixov: Zu. Éksper. 
Teor. ie, 34, 621 (1953). 
(6) V. D. Hopper, S. Biswas and J. F. DAarBY: Phys. Rev., 84, 457 (1951). 
(7) M. Demeur, C. DirwortH and M. ScHONBERG: Nuovo Cimento, 9, 92 (1952). 
(8) A. DEBENEDETTI, C. M. GARELLI, L. TALLONE and V. VIGONE: Nuovo Cimento, 
+ 


4, 1142 (1956). 
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It is also worth while by note that in the paper of EDWARDS et al. (*) sym- 
metry in complex collisions in the CMS is a basic assumption used in the 
discussion of results on elasticity. Although the energy of their jets is much 
higher than the energy of the jets investigated by us, the symmetry found in 
the present work may support the validity of the basic assumption made by 
EDWARDS et al. 


3. — Energy and transversal momentum distribution of shower particles. 


Assuming that all particles are mesons we have calculated from the ener- 
gies of the shower particles as measured in the LS and from the values obtained 
for y,, the energy distribution in the CMS and the transversal momentum 
distribution for all particles of the seven interactions (Fig. 1 and 2). The 
energy distribution has a maximum at an energy somewhat higher than m_e? 
which was predicted by the Heisenberg theory of multiple meson production (1°). 


wh 
20F 


al 
‘ 
0 - — pis A io lei 
7 04 06 08 12 14 16 18GeV Wir we ae A Me at OU 
Fig. 1. - Transversal momentum distri- Fig. 2. — Energy distribution in the 
bution. C.M. system. 


The transversal momenta have a mean value of about 03. GeV/e which 
is in rough agreement with values found by other authors. At such relatively 
low primary energies, however, the energies of secondary particles are also 
low and the fact that we have measured only shower particles with I < 1.51 
affects the form of the energy and transversal momentum distributions. 

Thus the above distributions are strongly biased by the mode of selection 


of the shower particles. Furthermore, owing to the large spread of the primary 


min 


(9) E. Epwarps, I. Losry, D. A. Perkins, K. Pinkau and I. Reynorps: Phil. 
Mag., 3, 237 (1958). 

(7°) W. HEISENBERG: Vorträge über kosmische Strahlung (Berlin-Gottingen- Heidelberg 
1953). | 
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energies of the individual interactions the building of a composite energy and 
momentum distribution is not justified and thus no definite conclusion can 
be drawn from the distributions found here as to the validity of multiple meson 
production theories. 


RAGS Ss UO NEMO Ce 


È stata misurata la distribuzione angolare ed energetica delle particelle di uno 
sciame di sette interazioni nucleari, prevalentemente collisioni complesse. Il fattore 
di Lorentz y,,, del sistema del centro di massa (CMS) è stato determinato dall'energia 
dei secondari e il fattore di Lorentz del sistema di simmetria y,,,, è stato ottenuto 
dalla distribuzione angolare. Si è trovato che %sm/Yem= 0.92-40.15 corrispondente 
all'emissione simmetrica delle particelle dello sciame nel CMS. Sono anche state deter- 
minate la distribuzione energetica nel CMS, nonchè la distribuzione dell'impulso tra- 
sversale delle particelle dello sciame. 


(*) Traduzione a cura della Redazione. 
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IL NUOVO CIMENTO Moms Sele INS el 1° Gennaio 1959 


On the Interpretation of Cosmic Ray Jets: The “Funnel” Model. 


T. GOZANI and K. SITTE 


Department of Physics, Technion, Israel Institute of Technology - Haifa 


(ricevuto l’11 Settembre 1958) 


Summary. — Arguments are presented in favour of the interpretation 
of cosmie ray jets as composite collisions, in a single act, of the primary 
with a compound of target nucleons. In particular, it is pointed out that 
the identity of «successive » collisions is lost at high energies, because 
the expanding meson cloud enwraps the compound before ejecting the 
bulk of the secondaries. According to this view, the primary energies 
assigned to the jets by the customary procedures are too low in many 
cases, and especially for events with large numbers N, of heavy prongs. 
The resulting difference in the true primary energies of the jets with, 
say, N,<3 and N,>3 accounts for the apparently excessive frequency 
of the events of low N,,. But if one wishes to retain the model of the nucleus 
as an aggregate of corpuscles, the cross section of the avalanche of interact- 
ing nucleons will increase with penetration, and the usual picture of 
the interaction volume shaped lke a cylindrical tunnel must be replaced 
by that of an inverted funnel. This «funnel model» predicts a stronger 
dependence of the multiplicities on the atomic weight of the target nucleus, 
and seems to fit the—rather scanty—experimental data somewhat better. | 
It is then shown that according to this model, the larger part of the 
observed multiplicity spread is not due to the «fundamental fluctu- 
ations » in the emission of secondaries in a nucleon-nucleon collision, but 
to the randomness in the number of particles involved. For nucleon- 
nucleon interactions at about 5-10!2 eV, the « fundamental fluctuations » 
can be approximated by a Gaussian with a width o = 6.5. 


1. — Although so far no completely satisfactory theory of meson production 
in high-energy nuclear collisions has been developed, wide agreement has now 
been reached about a number of general features which must be adopted by 
any such theory. In particular, it is now evident that Fermi’s original as- 
sumption (1) defining as «collision volume » the overlapping part of the col- 


(1) E. FERMI: Progr. Theor. Phys., 5, 570 (1950); Phys. Rev., 81, 683 (1950). 
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liding meson fields, cannot be maintained and must be replaced by a model 
in which collision and ejection of meson secondaries are considered as two 
essentially separate acts. The first only provides the energy necessary to 
« heat » the field, and the second takes place only after this field has expanded. 
to a much larger volume. This picture forms the basis of both Heisenberg’s (2) 
and Landau’s (*) theories, and of their various modifications proposed to meet 
the experimental criteria. 

However, a comparison with the experiment is further complicated by the 
fact that empirical data are available only from the analysis of jets in emul- 
sions or from air showers: both referring not to pure nucleon-nucleon encounters 
but to the more complex phenomenon of the traversal of a nucleus by an 
energetic particle. Cascade theories describing this process have been worked 
out by HEITLER and TERREAUX (‘), and by MOCUSKER and ROESSLER (5). 
A more general approach has been used by Coccont (*), and by BELEN'KIJ 
and MILEHIN (‘). In all these cases, the phenomenon is treated as the col- 
lision of the incident nucleon with all the particles contained in the «tunnel » 
or «tube» which its passage cuts through the nucleus; Coccont viewed the 
problem as a simultaneous interaction of the primary with n nucleons and 
applied the Fermi model to this «composite collision », while BELENK'1J and 
MILEKHIN based their work on Landau’s theory. 

More recently, PERKINS and his collaborators (§) have cast doubt on the 
applicability of this model of composite collisions, claiming from their analysis 
of high-energy jets that the majority of these events should be interpreted 
as nucleon-nucleon collisions. In the present paper, the description of a 
nucleon-nucleus encounter as a single interaction between two unequal partners 
will again be adopted, and used in conjunction with the picture cf the expand- 
ing meson cloud. Accordingly, we shall begin with showing by a very simple 
argument that if the latter view is taken, the former assumption—the com- 
posite-collision model—must necessarily be accepted. 

Let <p,), be the average momentum of the secondaries, in the center-of- 
mass system of the collision of a nucleon with a compound of n nucleons, and 


(2) W. HriseNBERG: Zeits. Phys., 113, 61 (1939); Zeits. Phys., 126, 569 (1949); 
Suppl. Nuovo Cimento, 6, 493 (1949). 

(8) L. Lanpav: Izv. Akad. Nauk SSSR, 51 (1953); S. Z. BELEN'K1J and L. LANDAU: 
Suppl. Nuovo Cimento, 3, 15 (1956). 

(4) W. HerrLER and C. TERREAUX: Proc. Phys. Soc., A 66, 929 (1953); C. TERREAUX: 
Helv. Phys. Acta., 24, 551 (1951). 


(5) F. C. RozssLer and C. B. A. McCusKER: Nuovo Cimento, 10, 127 (1953). 

(6) G. Cocconi: Phys. Rev., 93, 1107 (1954). 

(1) $. Z. BeLen'K1J and G. A. Mivenin: Zurn. Éksp. Teor. Fiz. USSR, 29, 20 
(1955) 

(8) ) 


B. Epwarps, J. Lusty, D. H. PERKINS, K. PINKAU and J. ReyxNoLDS: Phil. 
2201001950) 
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(Yoo cp), their average energy, all expressed in units of nucleon rest 
mass. Meson emission from the heated cloud will then take place only when 
the cloud’s linear extensions have reached a dimension of about h/<p,>. On 
the other hand, the width of the target nucleus, in the C-system, will be of 
the order Ry: A*/(y,),, where (yn), stands for the energy of the n-particle com- 
pound in the C-system, measured in nucleon rest masses: 


(1) (yng = (y + n)/(Qyn + n° +1}. 


y is the energy of the incident nucleon in the laboratory (L-) system. Evi- 
dently, for sufficiently large primary energies y the meson cloud will encompass 
the entire nucleus before emission takes place, and then the concept of a 
sequence of interactions with individual nucleons, or of an intranuclear cascade, 
becomes meaningless since no distinction can be made between successive 
collisions. One can easily verify that in the region of jet energies (y = 500), 
this is true for all emulsion nuclei. 

Another simple argument against the model of an intranuclear cascade 
as a succession of interactions, and for the idea of a «delayed » meson emis- 
sion from the expanding cloud—though not necessarily for the composite- 
collision picture—, follows from the experimental evidence that the interaction 
mean free path in the emulsion of jet secondaries has been found to be very 
nearly geometrical. Hence it is difficult to see why only nucleons, and not all 
shower particles, should participate in the intranuclear cascade. But if this were 
the case, the number of meson-producing encounters would increase very rapidly 
with increasing nuclear dimensions, and their contribution to the total mul- 
tiplicity of the collision would be very appreciable in spite of the comparatively 
lower energy of the x-meson secondaries. The first collision of a 1013 eV 
primary, for instance, would produce about 40 nucleo-active secondaries 
(neutral particles must be included in considerations of the intranuclear cas- 
cade), and even if for the secondaries of the forward cone Heitler’s assumption 
of «plural» collisions (*) offers some reduction, the second-generation collisions 
of the mesonic secondaries would add a high multiple of this number: a most 
unsatisfactory prediction. 

If, on the other hand, the model of meson emission from the expanding 
cloud is accepted, the only real difference between nucleon-nucleon encounters 
and nucleon-nucleus, or composite, collisions lies in the amount of energy 
transferred to the cloud, and hence available for secondary production. In 
other words, one might say that even in a nucleon-nucleon collision the x- 
mesons originated in the interaction contribute to the secondary production 
—in the «cloud stage »—just as in a nucleon-nucleus collision. For example, 
if in accordance with an intranuclear-cascade model distinctions can be made 
between z-meson-nucleon interaction and nucleon-nucleon encounters—pos- 
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Sibly in a distinctive composition of the shower—, events originating in light 
nuclei might differ from those initiated in heavy nuclei where the contribution 
of the z-mesons would be stronger. No such difference can exist according to 
the expanding-cloud picture. 

The following discussion will, therefore, be based on the two assumptions 
of composite collisions and of secondary emission from the expanded meson 
cloud. Together with the general kinematical relations, these two elements 
already determine all the essential features of the method. In particular, the 
dependence of the numbers n, of secondaries produced in a collision with a 
compound of n nucleons on this number n of participants is now fixed, without 
Specific assumptions concerning the details of the underlying theory of meson 
production. For it is logical to conclude from our general assumptions that 
the average energy <y,> of the secondaries should only depend on the energy 
available per nucleon in the C-system, but not on the way in which this energy 
became available. One may also regard this assertion as a plausible extension 
of the conclusion of KAMEKO and OKAZAKI (°) to processes other than nucleon- 
nucleon collisions. If we write <y,>, = g(y/n), it follows from the well-known 
definition of the energy available in the C-system, W., 


(2) W, = (2ny)! (y >n) 
that 
(3) n, = (2ny)*/p(y/n) . 


Since the energy dependence of n, can be approximated by a power law, 
yo, Will, in a similar approximation, also be represented by a relation 


NA 


(4) CV gh A (PM) 
from which one has 


(5) N= beni e 


The empirical y*-law for the multiplicity thus demands a n*-dependence of 
the multiplicity on the number of participants. 

However, this latter statement requires some clarification before it is ap- 
plied to a discussion of the experimental material. It must be pointed out 
that the argument is valid only if the primary energy is properly determined, 
i.e. under the assumption of a composite collision. All the customary methods 
of energy measurement, though, are based on the kinematics of a nucleon- 
nucleon collision. Take, for instance, the simplest case, applicable when the 


— 


(9) Su. Kamexo and M. OKAZAKI: Nuovo Cimento, 8, 1 (1958). 
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relation B_/B,— 1 of the velocities B_ and B, of secondaries and C-system holds: 
symmetry of emission in the C-system then demands that the angle 0, con- 
taining one-half of all jet particles in the L-system obeys the relation 


(6) tg 0, = 1/y, = (2/7). 


Under the conditions of a composite collision, and owing to the new relation 


between the energies y, and y in C- and L-system (always in extreme relati- 


vistic approximation), 


Ye = (yl2nŸ 


I 


the symmetry condition imposes on 0, the demand 
(8) te 0, = (2n/y). 


The true primary energy y is, therefore, larger by a factor n than the energy yy 
computed under the assumption of a nucleon-nucleon collision. This, together 
with equ. (5), defines the relation between the multiplicity n,, the number n 
of participants, and the apparent primary energy yy: 


L—m 


(9) Ne = BN aye 


We shall now present a purely statistical argument in favour of this view. 
It is derived from an analysis of the energy spectrum of 88 high-energy jets 
induced by nucleonic primaries: MCCUSKER and ROESSLER (!°), BRISBOUT 
et al. (11), EDWARS et al. (*), TAKIBAEV et al. (12), ALPERS et al. (3°), GUREVICH : 
et al. (4), divided into two groups according to the number N, of heavy prongs: 
N,<3, and N,>3. It has become customary to consider the first group 
as representing essentially only the nucleon-nucleon collisions among the jets, 
and the second as consisting of the more complex interactions. 


(1°) C. B. A. McCusKER and F. C. RoEssLer: Nuovo Cimento, 5, 1136 (1957). 
(11) F. A. BRISBOUT, C. DAHANAYAKE, A. ENGLER, Y. Fusimoro and D. H. PERKINS: 
Phil. Mag., 1, 605 (1956). 

(1) Z. S. TAKIBAEV, P. A. Ustk and M. G. Anronova: Dokl. Akad. Naulk 
USSR, 111, 341 (1956). 

(13) V. V. ALPERS, I. I. Gurevicn, V. M. Kuruxova, A. B. MiSAKova, B. A. NI- 
KOLSKI and L. V. Surova: Dokl. Akad. Nauk USSR, 112, 33 (1957). 
Gat GUREVICH, A. P. MISAKOVA, B. A. Nrkotsxy and L. V. SurKava: Zurn. 
Hksp. Teor. Fiz., 34, 267 (1958). 
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The results are shown in Fig. 1. It is seen that the spectrum for the two 
groups is practically the same, and it may be noted for reference later on that 
the energy distribution of the events listed differs significantly from the pri- 
mary spectrum of cosmic ra- 


diation at these energies: 
its initial slope corresponds 
to a power law with an ex- 
ponent ~ 0.5, and only at 
the highest energies—where 
unfortunately the statistical 
uncertainties become quite 
large—it appears to ap- 
proach the exponent of the 
primary spectrum, > 1.5. 


Since in this energy region 


the collision cross-section 


2 3 4 
is undoubtedly constant, the 10 10 10 ? 
deviation must result from Fig. 1. — Integral energy spectrum of the 88 jets. 
a strong scanning bias in Empty circles: N, < DE full circles: Ne 5. Aie 


dashed line represents the spectrum of the group 
N,<3, shifted so as to increase the relative fre- 
quency of this group to about 6094. 


favour of events of high 
multiplicity, 4e. in the 
average of high primary 
energies. 

However, at this stage we are concerned with another point: the apparent 
predominance of the «nucleon-nucleon » group N,<3 over the «nucleon- 
nucleus» group N, > 3. About 60% of all the jets quoted in the literature 
belong to the first class, and only about 40% to the second. This is contrary 
to exnectations: given the value of the «fundamental » cross-section © for 
nucleon-nucleon interactions, it is easy to compute from purely geometrical 
arguments the fraction of all jets involving only one target nucleon. AI col- 
lision with the hydrogen of the emulsion belong into this group, and all those 
peripheral collisions with other nuclei in which the length of the « tunnel » is 
too small to permit in the average more than one interaction. A similar estim- 
ate for showers of lower primary energy has been made previously in an 
attempt to evaluate the effect of second-generation collisions (5); for the pre- 
sent study, the method had to be slightly modified since the target nucleon 
is itself capable of further secondary production. Details about the method 
of calculation will be given in the last section; here it may suffice to state that 
its result is a probability of (35--40)% for interactions of the nucleon-nucleon 
type. Moreover, it is interesting to note that this result is very insensitive 


(5) J. G. AskowirH and K. Sirre: Phys. Rev., 97, 159 (1955). 
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to the choice of the parameter o,, as long as no unreasonably small value is 
adopted (less than about 4 the geometrical value, or 30 mb). It is clear 
why this should be so: large 0, «shorten » the path but extend the total col- 
lision area, and small o, permit a longer path, while cutting down the total 
collision area. Over a considerable range of values for the « fundamental cross- 
section », the two effects very nearly cancel each other. 

Admittedly, the relative frequency of jets in the two groups alone would 
not constitute a safe basis for an attempt to re-interpret nucleon-nucleus 
collisions. Surely there is no precisely defined N, which distinguishes nucleon- 
nucleon and more complex collisions. But on the one hand, over one-half 
of the events in the first group has N,= 0, and all of them seem to agree 
so well in their typical properties with the features expected of nucleon-nucleon 
interactions, that it is difficult to ascribe to the group a large admixture of 
complicated collisions. On the other hand, this statistical argument is used 
only as a further support bearing out the conclusions already reached, and 
not as the sole foundation of the composite-collision model. 

The explanation suggested for this strange frequency ratio is the following: 
according to equ. (7) and (8), a jet originated in a collision of a primary of 
energy y with n nucleons, will have in the L-system symmetry properties equal 
to those of a nucleon-nucleon interaction at primary energy (y/n). Since in 
the determination of jet energies the simpler interaction model was generally 
assumed, the data given underestimate the true primary energies by a factor n. 
Correspondingly, in our graph of Fig. 1 the points belonging to the group 
N,>3 should be moved towards the right, to energies higher than those 
given by this factor n representing the average number of participating nucleons. 
One might even be tempted to turn the argument around, and to ask by which 
factor the primary energies must be raised, or the energy spectrum shifted, 
in order to give the «right » frequency ratio. The close similarity of the spectra 
of the two groups permits an answer, although of course its semi-quantitative 
nature as a very crude approximation is self-evident. The choice n & 8 would 
reduce the frequency ratio of the two groups to the expected one between 
2:3 and 1:2. This must be considered as an order-of-magnitude estimate, 
and quantitatively of not more value than a comparison of the multiplicities 
of the two groups. According to equ. (5) and (9), the second group should, 
at the same primary energy, have a multiplicity larger by a factor n. Actually, 
the ratio of the two multiplicities is about 2. It seems fair to say that the 
result thus lends support in a qualitative way to the composite-collision picture, 
while demonstrating the quantitative insufficiency of data and classification. 


2. — Having advocated in general terms the case for the composite-collision 
model, we shall now proceed to criticize its customary formulations. Both in 
McCusker’s cascade theory (°), and in Cocconi’s summary presentation (8), we 
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find the peculiar picture of a nucleus perforated, as it were, by the traversing 
nucleon. But surely this model can be adequate only if the «nuclear matter » 
has the properties of a fluid, divisible without discontinuity or at least built 
up of constituents much smaller than the projectile. If, however, we regard 
the nucleus as a compound of corpuscles of the same kind as the incident 
particle, then we run into troubles with the model of the cylindrical tunnel. 
We cannot have the primary slice off parts of the target nucleons lying inside 
the tunnel, and we cannot have non-integer « numbers of collisions »; but 
neither can we expect the nucleons to line up obediently along the tunnel 
so as to provide head-on collisions with all the particles in the collision 
volume. 

It must be noted that this difficulty is found not only in a model which aims 
at a description of the process as a step-by-step cascade. The simplified treat- 
ment in terms of composite collisions equally cannot ignore encounters: with 
particles not entirely inside the cylindrical tunnel. But once such a particle 
is hit and starts moving alongside the primary, the collision cross-section is 
increased—« the tunnel is widened », we may say—and further collisions will 
occur with an enhanced frequency. Again most of the target nucleons will 
lie near the periphery of the (non-cylindrical) tunnel, and thus after the en- 
counter in their turn add to the effective avalanche cross-section. In other 
words, if we take seriously the corpuscular constitution of the nucleus, we 
must discard the idea of a tunnel-shaped 
collision volume, and replace it by one re- 
sembling an inverted funnel. 

In order to estimate the quantitative 
difference between the «tunnel » and the 
«funnel» model, we have used a crude 
approximation which offers the advantage 
of great simplicity. First of all, we shall 
restrict the considerations to primary ener- 
gies for which the trajectory of the target 


ake i i : Fig. 2. — Schematic diagram of 
nucleons after the collision remains in a MERONI de PEN 
good approximation parallel to that of the D Da on 
primary. The widening of the «funnel» is (ro), and of the approximation 
thus a consequence of the randomness of used (circle of radius 7,4). 


the impact parameter only, and not of the 

divergence of the avalanche of nuclear matter. Secondly, although more pre- 
cisely the funnel should be pictured as a multiple-tube aggregate, we shall 
approximate its cross-section by a circular one of equal area (Fig. 2). Using 
as definition for the «mean free path in nuclear matter» À 


(10) C= 40/5 Vo; 


E 3 - Il Nuovo Cimento. 
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where o is the—variable—funnel cross-section and R, the nucleon radius, we 
can compute the average development of the funnel, i.e., the value of o and A 
as a function of the path length s in the nuclear matter, on the assumption 
of average impact parameters <b) = 2(r+7,)/3. It turns out that for funnel 
lengths not exceeding 6R,, a linear approximation may be used: 


(11) r = ro[1 + a(s/A)], 


y is the funnel radius, 7, the radius connected with the fundamental CTOSS- 
section o) by È = 0, and À, the corresponding mean free path according 
to (10). For 6) —2/3-0,, © 40 mb, the best value for the constant is «= 0.38. 
The deviations of the correct average radii from (11) are nowhere larger than 
about 10%. 

With these radii, we can now compute the funnel volume V, and hence 
the number of participating nucleons n — V/V;. Except for near-central and 
for peripheral collisions, we take as funnel volume that of a cut-off cone of 
a height equal to the length of the primary path. For an impact parameter b 
with respect to the nucleus, one has 


yea aed = 
(12) n = VI", =" {[1 + 20 (Ro Ap) VA — (b/R)2]5 — 1}. 


In the case of central or peripheral collisions, this approximation becomes 
poor and has to be replaced by a somewhat improved calculation. 

Results for a number of nuclei, and for o) = 60, 40 and 30 mb, are shown 
in Fig. 3a-3d. Attention is drawn to the dependence on the atomic weight A 
of the target nucleus which is much stronger than in the case of the tunnel 
model. This fact is more clearly demonstrated in Fig. 4 where the values of 
n(A), obtained at average impact parameters <b> according to tunnel and 
funnel model, are compared. | 

Thus, a check on the dependence of the multiplicities n, in a jet—which, 
according to (5) and (9) reflects the complexity of the collision—should offer 
the simplest test for the relative merits of tunnel and funnel model. Unfor- 
tunately, owing to the paucity of the empirical material available, this test 
remains inconclusive at this time. The only relevant data are those obtained 
in the emulsion chamber work of KAPLON et al. (1); and even of those, only 
the analysis of jets originated in copper yields statistically significant results. 
But unfortunately, copper is one of the elements least suitable for a statistical 


(5) M. F. Kapton, D. M. Rrrson and W. D. WALKER: Phys. Rev., 97, 716 (1953): 
M. F. Kapton and D. M. Rrrson: Phys. Rev., 88, 386 (1952). 
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comparison with emulsion material whose mean atomic weight, properly aver- 
aged, is very close to the 64 of Cu. Moreover, the usual practice of assigning 


n 6 OR, 


Fig. 3. — Number of nucleons participating in a composite collision, as a function of 
the impact parameter b and for three values of 6,, in four elements: a) N; b) Fe; 


c) AgBr (A = 96); 


d) Pb. 


the jets of N, > 3 to heavy target nuclei, and the others to light nuclei, evi- 
dently fails in this case since the group N, <3 includes peripheral collisions 


in Ag and Br as well. 

It follows that while it is pro- 
bably correct to regard jets with 
many heavy prongs as originated 
in a heavy nucleus, these events 


Fig. 4. — Dependence on the atomic 

weight A of the average number of 

nucleons participating in a collision, 

for: a) tunnel model (A*-law); b) fun- 

nel model, o = 30 mb; e) funnel 

model, o, = 40 mb; d) funnel model, 
O) 00: 


10 20 100 200 A 
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do not represent the sum total of collisions with Ag or Br nuclei but only those 
with comparatively small impact parameters. To compare them with an 
unbiassed sample of interactions in Cu nuclei is, therefore, misleading. 

Yet a rough estimate can be made if we exclude peripheral collisions by 
comparing events occurring with an average impact parameter nm = 
= (R,.5,+%)/3 in the heavy emulsion nuclei, with the unbiassed average 
D>, = 2(Re,+%)/3 in copper. In order to do this, the values of B-n = nisi, 
have been computed for the all showers, using the empirical y*-lJaw for the 
multiplicities. The results are: 


Emulsions, all N,: don == BP 
Emulsions, N, <3: no) = 1/1 
Emulsions, N, > 3: AIA = DI, 
Copper: CB DT — 49); 


The ratio of the N, > 3 group to the copper data may be compared with 
the values computed according to the tunnel model, and according to the 
funnel model with three values of o), listed in Table I. While it should be 
emphasized that owing to the arbitrariness of the assumption for <b> no 
undue weight must be given to the calculated values, the much better agree- 
ment of the funnel data with the observations is indicative if not convincing. 
More emulsion chamber data would be highly desirable. 


TABLE I. — Ratios Npm/Ncy of the number of nucleons involved in a collision in heavy 
emulsion nuclei, at b = (Rpm + Y9)/3. and in copper at b = 2(Rca + r,)/3 according to 
tunnel and funnel model. 


| Nem/M Cu 


tunnel (At — law) 1.45 | 
funnel ©, = 60 mb 2.25 

) o, = 40 mb 2.0 

) 6) = 30 mb 1.8 
Experimental 1.8 


Another test for the choice of tunnel or funnel model could, in ‘principle, 
be based on Cocconi’s suggestion (5) to determine the « apparent » n from mul- 
tiplicity and half-angle. But without a much more accurate knowledge con- 
cerning the probabilities of multiplicity fluctuations, no safe conclusions can 
be drawn from such an analysis. The problem of fluctuations will be discussed 
in the next section; here it may suffice to state that with the ‘present material, 
and in particular in view of the scanning bias, no definite answer could be given. 

Two more general consequences of the composite-collision model may still 
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be pointed out. Firstly, it follows from the assertion made that the primary 
energies were, in many cases, underestimated by a large factor, that the 
inelasticities determined for such events are too high by the same factor. No, 
or no appreciable, correction need be made for nucleon-nucleon collisions, 
i.e. essentially for the group N, < 3; but for the others, the estimates reported 
should on the average be reduced in the same ratio as their multiplicity exceeds 
that of the first group. As it is seen from the summary above, for the jets 
originating in emulsion nuclei this ratio of the mean multiplicities is 3.55/1.4 & 
= 2.5. The question of a possible relation between the inelasticity and the 
number of heavy prongs N, is apparently not yet clearly answered by the 
emulsion experiments so far reported; more recent studies have cast some 
doubt on the validity of the affirmative conclusion of BERTOLINO (17) and 
of HOANG (18). 

Secondly, we wish to emphasize that the funnel model of composite col- 
lisions is by no means incompatible with a two-centers model of the kind sug- 
gested by TAKAGI (1), KRAUSHAAR and MARKS (2°), and ZACEPIN (21). Strong 
evidence favouring such a model has quite recently been published by the 
Polish-Czech co-operative emulsion group (22). In fact, the asymmetries re- 
sulting from this model are one of the most useful criteria guiding the analysis 
of individual events on the basis of the composite-collision model. A more 
detailed report concerning this point will be submitted later. 


3. — The occurrence of very large multiplicity fluctuations at equal half- 
angle, and hence presumably equal primary energy, has often been commented 
upon, and has sometimes been quoted as an argument in favour of a purely 
« Statistical » theory of meson production. It is evident that if the composite- 
collision viewpoint is adopted, the argument for a wide spread in a nucleon- 
nucleon collision at fixed primary energy loses strength. In this model, there 
are four reasons for fluctuations in the number #, of shower particles in a jet: 


i) The differences in the path length of the primary, defined by the 
impact parameter b with respect to the nucleus, and described by the prob- 
ability p,(b)db of a collision within the impact parameter interval (b, db). 


ii) The randomness of the collision process along this path s (or the 
« density fluctuations » of the nuclear matter on s), expressed by the proba- 


17) G. BERTOLINO: Nuovo Cimento, 2, 1130 (1955); 3, 141 (1956). 

) T. E. Hoane: Journ. de Phys., 15, 337 (1954). 

Ware tAksGr = Prog. UneormPiys:, 1. 231952): 

) W. L. KRAUSHAAR and L. J. MARKS: Phys. Rev., 93, 326 (1954). 

) V. I. ZAcEPIN: Suppl. Nuovo Cimento, 8, 585 (1958). 

22) P. Ciox, T. CocHan, J. GreruLA, R. Honynsxi, A. JURAK, M. Migsowicz, 
T. SANIEWSKA, 0. STANISZ and J. PERNEGR: Nuovo Cimento, 8, 166 (1958); 10, 741 
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bility p,(s) that while traversing the nucleus, the primary will collide with 
exactly n target nucleons. 

iii) The «fundamental fluctuations » in the emission process, 2e. the 
probabilities p,(n,y) that in a collision of a primary of energy y with n nu- 
cleons, exactly N secondaries are ejected. 

iv) The deviations from the mean charged-to-neutral ratio, characterized 
by the probability p, (N) that exactly n, out of the N secondaries are charged 
particles. 

It then follows for the probability distribution P,(n,, y) of the jet mul- 
tiplicities at energy y: 


bmax 
(13) : Pi(ns,y)= | SY > P.(b)-p(s) Prin) Pn,(N) db. 
È, n=1 N=n, 


0 
Three of the four probabilities are readily determined: 
a) The « geometrical » probability p, is given by 
(14) P,(b) db = 2xb db/x(R,A*+1,)? . 
b) p,(s) can be evaluated, at least approximately, on the basis of the 
determination of the mean free path 2, for the n-th successive « generation » 


of collisions described above. Straightforward application of the elementary 
caleulations gives, first, 


8 A 


(15) AGE lex [— ol Il exp [— Ya/ A2] dYa/ 2 ... 
0 ò 
Sie Una) 
… exp [— Yn/An] dYnfAn'@xp[ (8 — yy —..- — Yn) Anti | - 


0 


The integration gives after some 
transformations 


n +1 
(GS) Sr ES PE 
î=1 
n+1 
‘exp [— 8/44) TT (4 — 2) . 
ki 


Fig. 5 shows some of the p, as a 
function of s. 


Fig. 5. — Probabilities p, of collisions 
with exactly » target nucleons, as a 
function of the path length s/4,. 
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e) The last factor p, (N) is the well-known binomial function 
(17) Dr EN) EN 1 gL gt NL 


where q is the mean value of the charged-to-neutral ratio, for which the fraction 
3 may be used. | 

Hence it is possible, at least in principle, to determine the unknown «fun- 
damental fluctuations » p,(n) if the multiplicity distributions P,(n,, y) are given. 
Practically, however, we meet a formidable obstacle: the effect of the scanning 
bias. 

It has been pointed out above that already the deviation of the jet spectrum 
from the energy distribution of the primaries renders proof of a significant 
Scanning bias in favour of events of high multiplicity. Consequently, the ob- 
served multiplicity distributions P(n,, y) will differ from our theoretical P,: 
large multiplicities will appear enhanced and low multiplicities diminished. 
Moreover, the empirically determined averages », will be higher than the true 
mean multiplicities <n,>, and the error thus made will itself depend on <n,>, 
i.e. on the primary energy. Only for large y’s it appears from the spectrum 
of Fig. 1 that the bias is small. It may here be noted that not even the empir- 
ical y'-law for the average multiplicities can be retained if the bias is indeed 
appreciable. The true energy dependence would be steeper than the one de- 
rived from the present material. 

Since the scanning bias can be assumed to depend on the number of shower 
particles n, only, it is of course theoretically possible to eliminate, or to de- 
termine, the «bias function » b(n,). For if a(y)dy is the differential primary 
spectrum, and s(y)dy the observed jet distribution, then we have, owing to 
the fact that the primary interaction cross-section is not energy-dependent, 


(o 0) 


(18) s(y)dy = > aly): Pons, y):D(n,) dy, 


ns 


or for the observed multiplicity distribution 
(19) P(ns, y) == Poe y)-b(ns) . 


In practice, the experimental material is by far insufficient for such a detailed 
test of the assumptions, and it was not deemed worth while to carry out the 
necessary cumbersome calculations. Instead, an attempt was made to esti- 
mate the effect of the «fundamental fluctuations » by another, less accurate 
but better documented, method. 

One notes, first of all, that for nucleon-nucleon collisions only the third 
and fourth of the fluctuation-inducing effects can occur. This simplifies the 
calculation, but leaves the bias problem still unsolved. But since we are then 


40 T. GOZANI and K. SITTE 


dealing with events of the lowest multiplicities, one may hope that although 
just there the bias is absolutely at its largest, its effect will not vary too much 
over the—reduced—range of multiplicities involved. In a first approximation, 
we may therefore assume a constant bias, and hence for the purpose of esti- 
mating the influence of the «fundamental fluctuations » on the multiplicity 
distribution, identify the observed P(n,, y) with the theoretical P,(n,, y). 

A glance at the experimental data supports this view, and suggests one 
further measure. In order to be as sure as possible of the character of the 
interactions, only events with N,=0 were selected for this study. The in- 
fluence of the energy dependence of the multiplicities was minimized by com- 
puting the values of n,(y/y)*, where y is the average energy of the 20 jets of 
this subgroup. The results are shown in the histogram of Fig. 6a, and they 
reveal an interesting feature: it appears that they consist of two separate 
groups, the larger one of 17 jets centering at an ©, between 10 and 15, and 
a smaller one at n,~ 30. It is most suggestive to consider the smaller group 
as resulting from a collision of a primary with two target nucleons, giving twice 
the number of secondaries. These events will, therefore, be excluded from 
the further analysis which is thus restricted to 17 jets with n,=9.5 at 
y=5290. Its multiplicity spread is no longer very large: the values of 
ns(7/y)* ‘ange from 4.3 to 17.5. The error which derives from neglecting the 
scanning bias is, therefore, unlikely to be very serious. 


(a) 
VA # 9 
30 30 
20 20 
10 "7 
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Fig. 6. — a) Histogram of n,(y/y)* for all jets with N, = 0. 5) The same histogram after 
elimination of the events with the highest multiplicity, and the distribution calculated 
with o = 6:5. 


An attempt was then made to fit the experimental data by assuming a 
Gaussian for the «fundamental fluctuations »: 


(20) Py = (10V2x) exp [— (N —N)/20?]. 
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Together with (17), this leads to a probability P(n,) to find n, charged jet 
secondaries 


(ce) 


| Pet EN. O 
(21) P(n,)=(q"/n.toV2n) | en Oren (NN) Parle 


0 


where for our sample, N —15. The experimental histogram, and the best- 
fitting calculated curve (o = 6.5), are shown in Fig. 6b. The marked deviation 
of the curve from a simple Gaussian illustrates the effect of charge fluctuations. 

If one assumes that the same Gaussian (20), with N a known function 
of y, represents the «fundamental fluctuations » for all energies per nucleon, 
it is possible to work out the complete multiplicity distribution P,(n,, y) ac- 
cording to (13). However, neglecting the bias function over this wide range, 
—indicated in the histogram of Fig. 7b—we can certainly, at best, hope for 


D/o (b) 


Ww 
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SS = 
0 20 0 6005 0 20 550 6005 


Fig. 7. — a) Contribution to the various n, of collisions with n nucleons, and resulting 
distribution. 5) Histogram of the observed distribution of »,(y/y)* for all jets, compared 
with the calculated resulting distribution, and with a Gaussian of width o = 6.5. 


a crude approximation and the calculations are extremely lengthy and tedious. 
The exact calculation was, therefore, not carried through at this time. Instead, 
the two decisive points we wish to emphasize will be demonstrated in an over- 
simplifying model calculation: first, that the «fundamental fluctuations » ac- 
count for only a small part of the multiplicity spread observed in jets, and 
second, that a composite-collision model with a moderate «fundamental flue- 
tuation » can explain the wide multiplicity distribution. 

To prove the first point, a Gaussian with o=6.5 and N=20, corresponding 
to the average of the complete sample, is plotted as a thin line in the histogram 


Zi 
= 
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of Fig. 7b. Clearly, the disagreement is too marked to be attributed to the 
neglected fluctuations in the neutral-to-charged ratio; at those large mul 
tiplicities, their effect is rather small. Next, we have chosen a «typical » fre- 
quency distribution of n, the number of nucleons participating in the compo- 
site collisions averaged over all emulsion elements. For simplicity of the cal- 
culation, the—infrequent—large n corresponding to central collisions in heavy 
elements have been omitted. The following values were chosen for the fre- 
quencies F, of an interaction with i nucleons: F,= 0.40; F,— 0.14; Fs= 0.20; 
F,= 0.14; F;= 0.08; F,—0.04 With this assumption the average multipli 
city (n,> in a nucleon-nucleon collision at the energies of our sample is deter- 
mined: it is related to the observed average ns by 


(22) N = DEAD. 
i 


For our sample, one has <n, & 8. Again disregarding charge fluctuations, 
the approximate multiplicity distribution is then given by 


(23) P*(n,) = (1/ov27)-Y Firexp[— (n, — i<n,>)?/20?] . 


Fig. 7a shows, as an illustration, the relative contributions of the collisions 
with different complexity (thin lines) to the resulting distribution (thick line). 
In Fig. 7b, this resulting distribution is plotted together with the experimental 
histogram. It is seen that the characteristic shape of the experimental curve 
is very well reproduced, and that is all that can be reasonably expected at 
this stage. 

In conclusion, it may therefore be said that the composite-collision model 
of cosmic ray jets seems to give a satisfactory explanation of all the features 
observed, including the multiplicity distribution. It is not necessary to postu- 
late the occurrence of very large «fundamental fluctuations » in a nucleon- 
nucleon interaction; they account for only a fraction of the total spread. 


RIASSUNTO (*) 


Si presentano argomenti a favore dell’interpretazione dei getti di raggi cosmici 
come collisioni composite del primario con un nucleone bersaglio composto che avven- 
gano in un singolo atto. Si rileva in particolare che l’identità di collisioni « successive » 


(*) Traduzione a cura della Redazione. 
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va perduta alle alte energie perchè la nube mesonica che si espande avvolge il com- 
posto prima di emettere la massa dei secondari. Secondo questa concezione le energie 
primarie attribuite ai getti dai procedimenti usuali sono in molti casi troppo basse, 
specialmente per eventi con grande numero N, di rami pesanti. La differenza risultante 
per energie primarie vere dei getti con p.e. N,<3 e N,> 3 spiega la frequenza appa- 
rentemente eccessiva degli eventi con basso N,. Ma se si desidera di conservare il 
modello del nucleo costituito da un aggregato di corpuscoli, la sezione d’urto della 
valanga dei nucleoni interagenti aumenterà con la penetrazione ed il quadro usuale 
del volume d’interazione a forma di tunnel cilindrico deve essere sostituito da quello 
di un imbuto rovesciato. Tale « modello a imbuto » fa prevedere una maggior dipen- 
denza delle molteplicità dal peso atomico del nucleo bersaglio e sembra meglio adatto 
ad accordarsi coi dati sperimentali « alquanto scarsi ). Si dimostra poi che secondo 
questo modello la maggior parte dello sparpagliamento della molteplicità osservato 
non è dovuta alle « fluttuazioni fondamentali » nell’emissione dei secondari in una 
collisione nucleone-nucleone, bensì alla casualità del numero di particelle coinvolte. 
Per interazioni nucleone-nucleone a circa 5:10! eV le «fluttuazioni fondamentali » si 
possono approssimare con una gaussiana di ampiezza 6 = 6.5. 
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Spurious Scattering of 6.2 GeV Protons 
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(ricevuto il 24 Settembre 1958) 


Summary. — Multiple scattering measurements of 6.2 GeV proton tracks 
in emulsion have been carried out for cell lengths up to 10000 microns. It 
was found that the spurious scattering contribution remained constant 
for cell lengths greater than about 3000 microns. Measurements on a 
group of parallel tracks showed that the region of the microscopie disloca- 
tions responsible for spurious scattering had dimensions on the order of 


1 mm to 2mm. 


1. — Introduction. 


In considering the multiple scattering of particles in their passage through 
nuclear emulsions, errors are introduced by apparent scattering from a number 
of sources. Corrections for these errors must be applied to measured tracks. 
One major source of this error, usually called spurious scattering, is suggested 
by Biswas, PETERS, and RAMA (1) to be due to microscopie dislocations within 
the emulsion. This spurious scattering has been confirmed and studied by 
FAY (?), LOHRMANN and TEUCHER (?), BRISBOUT et al. (4), and APOSTOLAKIS 


(*) Assisted by the Joint Program of the ONR and AEC and a grant from the 
National Science Foundation. 

(1) 5. Biswas, B. Perers and H. Rama: Proc. Ind. Acad. Sc., A 44, 154 (1955). 

(?) H. Fay: Zetts. f. Naturfor., 10a, 572 (1955). 

(*) E. Lourmann and M. TEUCHER: Nuovo Cimento, 3, 59 (1956). 

(4) F. A. BRrISBOUT, C. DAHANAYAKE, A. ENGLER, P. H. FowLER and P. B. Jonzs: 
Nuovo Cimento, 3, 1400 (1956). 
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et al. (9). While it is probably true that the dislocations originate during the 
processing of the emulsions it has been surprising that, with but few ex- 
ceptions (*5) all workers obtain about the same magnitude for the scattering. 
The present work was undertaken to extend these measurements (15) to 
larger cell lengths and also to learn something of the character of the emulsion 
dislocations. For this purpose Ilford type G-5 glass backed emulsions 2 in. 4 in. 
in area and 600 um in thickness were exposed to one pulse of the 6.2 GeV 
internal proton beam of the Berkeley Bevatron. The particles entered the emul- 
sion parallel to the surface so that the average track length was about 6 cm. 
The spurious scattering measurements made of these tracks are in general 
agreement with those of other investigators up to the largest cell lengths re- 
ported (15), 4000 um. The results herein reported show that, while the spurious 
scattering increases rapidly to cell lengths of about 3000 um, it appears to 
remain constant to the longest cell lengths used in this experiment, 10 000 um. 


2. — The measurements. 


The multiple scattering measurements were carried out with à standard 
medical type Tiyoda microscope that was attached to a heavy cast iron pre- 
cision stage of local construction. The positions of the track under observation 
were then measured with a filar micrometer eyepiece. 

The multiple scattering of the 6.21 GeV proton tracks in this experiment 
was measured by the method of FOWLER (5) in terms of second differences. 
As shown in Table I the mean of the absolute values of the second differences, 


TaBLE I. — Multiple scattering second difference as a function of cell length. 
| | 
| Cell No. of | Observed second Noise second Theoret. | Spurious | 
length over- differences differences second scattering | 
| (um) lapping D, DEA | diff. | second diff. | 
cells | | | O a. DR | 
| 250 917 1464 .004 | 0964.01 | 024 | .107+ .010 
| 500 |: 502 202+ .008 | .091+.01 .068 | .167-4 .011 
| 750 288 DES CON M o 3 CT 
| 1000 | 670 | 414+ .016 | .117+.01 -205 | .344+ .020 
| 2000 319 1.008 .056 | .110+.01 601 800-4 .070 

5 000 Bl) ie Ole eae sel 10.01 2.47 8424+ 4 | 
| 10000 SU RTE SR LI 42.01 7.24 3 lo, 80 


(5) A. J. APOSTOLAKIS, J. 0. CLARKE and J. V. MAJOR: Nuovo Cimento, 5, 337 
Nie 


(1957 
(6) P. H. FowLeR: Phil. Mag., 41, 169 (1950). 


46 F. W. FISCHER and J. J. LORD 


Dy, was determined for cell lengths from 250 to 10000 pm. As is standard. 
practice, individual second differences greater than 4 times the mean were 
not included in either the measured or theoretical mean values. 

Contributions to the second differences due to errors from thermal ex- 
pansion of the microscope, reading error, grain noise, and stage noise were 
determined by the method given by PETERS (1). Thermal errors were eli- 
minated by shielding and the contribution to D,, the measured second dif- 
ference, by the latter three sources of error is designated by D,,, in Table I. 
In all cases D,,y is quite small compared to D,: The corresponding theoretical 
value of the multiple scattering second difference, D,, was calculated from the 
Molière (7) theory. 

The spurious scattering second differences in this experiment were shown 
to have a normal distribution (+), therefore the mean of their absolute values 
is closely given by: 


cn nin 


KGS * 


di ES ET e e The calculated values of D,, are gi- 
ILL Se i ven in Table I. These values and 
CS Se Eee es those of other workers have been 

ale Toe plotted in Fig. 1. Curve © of this 
= Tite on | IT figure shows that spurious scattering 

I obs fee L e Lut errors become negligible for second 
RS | si differences greater than 1 um. 

La = 44 HH To explore the idea that spurious 
Is LL scattering is due to microscopic 
RS | | HH emulsion dislocations, measurements 
© | Ì were taken of 6 parallel tracks lo- 
QUES E00 cel, length 10 =. cated within a circle of 70 um dia- 


meter. It was then ‘possible to map 
out the local dislocations that give rise 
to spurious scattering. The results 


Fig. 1. — Multiple Coulomb scattering as 
a function of cell length. Curve A from 
Moliére theory. © measured second dif- 


ferences. Spurious scattering second dif- indicated that the regions of local 
ferences: X present work, @ Lohrmann and dislocations were of the order of 
Teucher, & Bristol group. 1 mm to 2 mm in diameter. This is 


consistent with curve © of Fig. 1 
wherein the D,, remains constant for cell lengths greater than about 2 mm 
to 5 mm. 


(7) G. MoLIÈRE: Zeits. f. Naturforsch., 3a, 78 (1948). 
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We wish to thank Professors Y. B. Kim and S. H. NEDDERMEYER for helpful 
discussions. The assistance of J. F. KIRKBRIDE and B. F. RUFFNER is grate- 
fully acknowledged. The exposure of the emulsions was made possible by 
Dr. E. J. LOFGREN and the Staff of the Berkeley Bevatron. 


RIASSUNTO (*) 


Sono state eseguite misure di scattermg multiplo su tracce in emulsione di protoni 
di 6.2 GeV adottando celle fino a 10000 um. Si è trovato che il contributo dello scat- 
tering spurio resta costante per lunghezze di celle superiori a circa 3000 yum. Misure 
eseguite su un gruppo di tracce parallele mostrano che la regione delle dislocazioni 
microscopiche che causano lo scattering spurio ha dimensioni dell’ordine di (1--2) mm. 


(*) Traduzione a cura della Redazione. 
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Summary. — We consider the possibility, suggested by P. ROTHWELL 
and T. GOLD in conversation with the author, that the radiation observed 
by Van ALLEN and co-workers at high altitudes is due to the decay 
electrons and protons from neutrons produced in the earth’s atmosphere 
by cosmic rays and stored in the earth’s magnetic field. Order of mag- 
nitude estimates for the densities to be expected are presented. Only 
scattering loss is considered. Using a lifetime of 3-10°s for loss through 
scattering, an upper limit of 10-? electrons/cm* near the earth and 
-05(R,/")° electrons/em? at large distances », in the equatorial plane 
is obtained. The proton density at large distances is .03(R,/7))? cm-?, 
for a lifetime of 102$. If plasma accelerations are not important, then 
the spectrum of electrons will be that of neutron B-decay. The protons are 
produced by fast neutrons coming directly from nuclear stars and their 
spectrum will be approximately that of the protons from stars. The 
electron density is sufficient to give a counting rate a few times larger 
than is observed, while the proton density is sufficient to give a counting 
rate 10* times higher than the observed lower limit. The lifetimes of 
stored particles are therefore probably much less than those given by 
scattering. There should be a strong latitude effect which is roughly 
estimated as proportional to cost. Reasons are given for believing that 
collective effects will reduce the density below this near the poles. 


1. — Introduction. 


An experiment by VAN ALLEN and co-workers (1) in the American Satellite, 
Explorer III, has demonstrated the existence of an intense zone of radiation 
(1) J. A. Van ALLEN, G. H. Lupwie, E. C. Ray and C. E. McILwaiN: IGY Satellite 


Re port, it oerves no. DI National AAC ademy of Sciences N ational Resea reh Council \\ À- 
5 € 
shineton 29, DAC 
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at altitudes above 1000 km above the earth’s surface. The radiation detector, 
a Geiger counter inside the satellite, does not permit a direct observation of 
the nature of the radiation; however, the sharp lower boundary of the radiation 
zone makes it very unlikely that the primary radiation is X-rays and strongly 
suggests charged particles. The high intensity of the radiation and its relative 
constancy in time suggests that particles are being stored in the earth’s mag- 
netic field. 

We review briefly the results of the perturbation theory of the motion of 
charged particles in a magnetic field. A detailed discussion may be found in 
the books of ALFVEN and SPITZER (**). We consider the motion of a charged 
particle in a region of space containing a magnetic field where the radius of 
curvature of the particle in the magnetic field is very small compared to the 
distances over which the magnetic field varies. We shall want to consider 
only fields which are constant in time, so the velocity of the particle is constant. 
If the magnetic field were uniform, the orbit of the particle would be a spiral 
around the lines of force, and the pitch angle y which the orbit makes with 
the lines of force would be constant. If the magnetic field is not uniform, 
the magnetic moment 


1 mv? sin? y 1 


(1) 2 Vip Bb 


of the particle due to its circular motion is constant. Here v cos y is the com- 
ponent of velocity parallel to the magnetic field, so particles are reflected from 
regions where 


1 sin? Yo 
2) FRS Bae 


If the magnetic field is not uniform around the particle orbit, then the par- 
ticle will be bent more strongly at one side of its orbit than the other and will 
precess in a direction perpendicular to the magnetic field and to the gradient 
of field. The velocity of the guiding center due to this precession is given by: 


me 
(3) = 0” cos? 


dI 
VTT 
Pa = MT ES 


For the rest of this paper we will approximate the magnetic field of the 
earth by a field of a dipole. This is sometimes not a very good approximation, 


(2) H. ALFVÉN: Cosmical Electrodynamies (Oxford, 1950). 
(3) L. Spirzer: Physics of Fully Ionized Gases (1956). 
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particularly when we compute derivatives of the magnetic field. Inclusion 
of higher moments of the earth’s field would be necessary to an accurate treat- 
ment of the shape of the radiation zone and its distortion of the earth’s mag- 
netic field, but these are beyond the scope of this paper. In spherical co- 
ordinates, the components of magnetic field are given by 


2u cos 0 


pe 


u sin 0 
73 


(4) B, , - B= 


E) 


where yw is the earth’s magnetic moment. The equation for a line of magnetic 
field is 


(5) r=, sin? 60, 


where » is the radius of the line at the magnetic equator. The intensity of 
the magnetic field at any point is given by 


(6) B = <V4—3sin?6 


and the intensity of magnetic field along a line can be expressed in terms of 
either co-ordinate by 


B= 


e 
“1 


I) V4— 3 sin? 0 a 1/ r 
re sin’ 6 A; RD, 
Along a line of force the magnetic field is weakest at the equator and increases 
toward the poles; therefore, particles whose pitch angle at the equator is not 
too small will be reflected from regions of stronger magnetic field as they pro- 
ceed toward the poles of the earth and will be in trapped orbits. We would 
now like an estimate of the time a ‘particle would be expected to remain in 
such a trapped orbit. The only process we will consider for removal of par- 
ticles from trapped orbits will be the scattering of particles in close collisions. 
This process is certain to occur and will set an upper limit to the lifetime. Well 
above the earth’s atmosphere it is believed that space is filled with ionized 
hydrogen with a density of about 1000 particles/cm?. In such a medium the 
lifetime for scattering through an angle comparable to 90° is given by (1) 


il 
pa Sr No(e?/2E) In (2 DIE) ) 


= 
co 
— 


t 


(4) L. SPITZER: Physics of Fully Ionized Gases (1956), eq. (5.13). 
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where N is the number of scattering centers per cm? and £, v are the energy 
and velocity of the particle. Table I gives some lifetimes of particles. 


TABLE I 
= e 5 —=# === == === = === 
| 
| t 
| Protons | 0.1 keV | a) CI E 
| 1.0 keV | 1 -10° s 
| 10.0 keV | 2.5-108 8 | 
| 1.0 MeV | 2 -1011g 
pr EIA 
| Electrons 10.0 keV | 6) 108s 
| 100.0 keV | I LOSES 
| | 5 MeV 35109 8 
| | 1.0 MeV TRTOLOES 


In Table II, we present certain characteristics of trapped orbits. 7, is de- 
fined above, R, is the radius of earth and atmosphere, 2 is the latitude at 
which the line through 7, reaches the atmosphere and y,, is the minimum pitch 
angle at the equator for trapped orbits. From Table II we see that the par- 
ticles of quite large pitch angles can be stored so that the scattering through 
90° is a reasonable value to take to give the lifetime. 


TABLE II. 
to/Ba A | Wm 
| 
1.05 ASB" 64.0° 
hath 268 55.0° 
182 | DATES 45) aie 
eS 35.4° DZ 
3.0 | 54.82 | JOE 
20 68.0° | Wee 


The lifetime of a 1 MeV electron against synchrotron radiation in a field 
a + gauss is comparable to, but not much less than, the value in Table I. 

The storage times in Table I are quite long; therefore, it seems likely that 
other scattering processes may be of importance and, in particular, if any plasma 
oscillations are set up they will be much more effective in throwing particles 
out of the storage region than close scatterings. Thorough consideration of 
such plasma oscillations, however, is beyond the scope of this paper since we 
wish merely to demonstrate that Van Allen’s zone of radiation represents a 
place where theories of plasma oscillations can be tested. 


ri 
10 


Or 
LO) 
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We note, however, that the lower boundary of the zone is stable against 
flute instability since the dominant force is the centrifugal force caused by 
the particle’s precession around the earth and motion along the curved line 
of force and this force is into the plasma (5). 

We must now consider how particles get into trapped orbits. The fact 
that the lifetimes are very long means that it is very difficult for charged par- 
ticles to diffuse into this region as well as to diffuse out of it. If the region is 
to be filled by charged particles diffusing from outside, the flux of particles in 
this region can only be time average of the flux incident on the outside of the 
region, and we shall have lost any advantage of storage. Therefore, it seems 
necessary to consider some other means of trapping particles in this region. 
It may be that clouds of particles from the sun can penetrate into the Earth’s 
field and then break up, leaving large numbers of particles in trapped orbits. 
Such a mechanism does probably contribute some particles, but it will not 
be considered in this paper. 

Another mechanism for trapping charged particles evolved in conversation 
between the author and P. ROTHWELL and T. GoLp, and was put in final form 
by GoLp. Cosmic rays produce large numbers of neutrons in the earth’s 
atmosphere. Some of these neutrons will diffuse out of the earth’s atmosphere 
and proceed away from the earth unaffected by further collisions or by the 
earth’s magnetic field. When these neutrons decay the resulting proton 
(~ 1 keV) and electron (~ 1 MeV) have, according to Table IT, a good chance 
of being trapped. In the next section we estimate the flux of neutrons upward 
from the earth due to cosmic rays. 


2. — Estimate of neutron flux. 


We review briefly the history of a neutron in the atmosphere (*). Many 
neutrons are emitted as evaporation partices from air nuclei struck by primary 
cosmic rays or their energetic secondaries and have energies between a few 
and 30 MeV. These energetic neutrons lose energy rapidly through inelastic 
collisions with air nuclei until they drop below the first excited state in ni- 
trogen at 2.3 MeV. After this they lose energy by elastic scattering, losing 
a fraction of their energy, — 2(m/M), per collision on the average, where m 
is the neutron mass, M is the mass of an air nucleus. When their energy 
has dropped approximately to thermal energy, they are absorbed by the ni- 
trogen of the air in the reaction !4N-+n +04 p. Some other reactions are 
important, but this dominates. Deep within the atmosphere the number of 


(5) M. RosenBLUTH and C. L. LONGMIRE: Ann. of Phys., 1, 120 (1957). 
(°) H. A. BeTHE, 8. A. Korrr and G. G. PLACZEK: Phys. Rev., 57, 573 (1940). 
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neutrons between £ and d£ is equal to the rate of production of neutrons 
times the time a neutron takes to lose energy d# and is therefore approximately 
proportional to 


x 
(9) na —— (collision time)= + Les 

(energy loss per collision) — 2(m/M)E v 

where / is the mean free path for the neutrons. Near the top of the atmosphere 
neutrons have à chance of escaping before being slowed down and the spectrum 
is somewhat distorted in that there are more high energy neutrons. To make 
an estimate of the flux of neutrons leaving the atmosphere, we use the formula 
of Bethe, Korff and Placzek (5) for the number of neutrons having velocity 
between v and v-+dv at a depth x in the atmosphere 


(0) see Ole) dr 


m_dv 1 al : (ei)? | (x + &% + (2/V3))?]| 
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Here q is the rate of production of neutrons at a depth x,.and L, given by 


1M f.. dE 
(11) Di AB) IL 


is the mean distance which a neutron of energy £ has diffused in slowing down 
from energy E, which we take to be ~ 2 MeV. 

The integral can be evaluated approximately as follows. At energies well 
above 1 keV, J is constant and several (4--10) times larger than at energies 
well below 1 keV (7). Hence in equation (10) we may take 
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It will prove more convenient if we modify equation (9) to take into ac- 
count capture of neutrons by air nuclei. The fraction of neutrons captured 
per collision is given by o(capt)/o(total) = ¢,/o, and so we may write 


aN M 0, 


dH 2mHo, 


Ta 


(13) 


(7) D. J. Hugues and J. A. Harvey: Neutron Cross Sections (New York, 1955). 
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o(capt) is proportional to 1/V.Z and so we multiply Bethe, Korff and Plazcek’s 
neutron distribution by a factor 
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where x=.05 eV and we have used (7) o, = (.34/V#)4 barns and o,= #(10) + 
+1(3.7) barns. 

Not much is known about the function q. At large depths in the atmosphere 
the slow neutrons in the atmosphere should be in equilibrium with production 
which indicates that q falls off with an absorption length of about 200 g/em?, 
i.e. about ten times longer than a neutron mean free path. Near the top of 
the atmosphere it may be expected to decrease more slowly or even to increase 
with atmospheric depth, because of the greater number of secondaries produced. 
In view of our ignorance concerning q and its probable slow variation we shall 
neglect its variation with «,. According to diffusion theory, the flux of neutrons 
leaving the top of the atmosphere is given by 


it. GN 
15 T= == 3 
(25) 3 Cx 


Then we may evaluate the integral in equation (15) at the top of the atmos- 
phere to give: 


1 M i ade = È fio; 
16 ff —_. = Ww EE al os — I ax =; == ax = / — 
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since //L<1 for most of the range of interest. 

The constant factor (M/m)ql may be evaluated using the results of So- 
BERMAN (5) who measured the counting rate in BF, proportional counters near 
the top of the atmosphere; that is to say, measured 


(17) G(@) -| N(x, v)o,0N,V dv, 


0 


where o, = 650/VE barns is the capture cross-section of 2B, N x is the number 


(5) R. K. SOBERMAN: Phys. Rev., 102, 1399 (1956). 
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of !°B nuclei per cm’, and V is the volume of the counter. In order to put 
this into a form which is as independent of q as possible we compute: 


| le Mg de M 1 
(18) mg = (02 PNA 060) (Patt, 
da |3-0 GF CH a9 m 4/27 L Va 


since //L is constant over the region where o, is large. We evaluate de/dp 
by drawing a straight line between the point ce — 0, pressure — 0 and an aver- 
age of Soberman’s highest measurements. The results are de/dp —0.4 counts/min- 
millibar at the equator; 2.1 counts/min-millibar at St. Paul, Minnesota (~ 55° 
geomagnetic); 4 counts/min-millibar near the pole. All of these measurements 
were made in the period 1952-54, so represent values when the solar activity 
is low or minimum. This gives for (M/m)ql, .014 neutrons/em? s at the equator, 
.-07 at St. Paul and .14 at the pole. 

The total number of neutrons leaving the top of the atmosphere/em? s is 
given for comparison by 


DA GIOIE ae dona 
la = = = Si aa a E RE ’ 
3m | DOT AA Ia, E È 0; O mi 


where the integral is taken from 0 to 2 MeV. We split the integral (19) also 
at H=1keV and put exp[—V«/H]=1 for energies above 1 keV. Then the 
integrals can be done easily using equation (11). The number of neutrons 
leaving the atmosphere is smaller than one might guess for the observed cosmic 
ray flux of one particle/cm? s incident on the atmosphere. The reason seems 
to be the decrease of mean free path at low energies, as discussed in con- 
nection with eq. (11), which allows the neutrons to penetrate deeply, then 
prevents them from leaving after they have been slowed down. We ought, 
therefore, also to consider neutrons which are produced in stars and which 
leave the atmosphere directly wtihout another collision. We review briefly the 
production of particles in stars. These results have been taken from the work 
of R. MCKEAGUE (*), who measured the stars produced in light elements in 
nuclear emulsion by 1 GeV protons. We assume that the production of neu- 
trons is the same as that of protons. In its first collision with an air nucleus, 
a primary proton will produce about .5 neutrons in the backward direction 
with energies in the range (0--30) MeV (black tracks), the maximum number 
being produced at (4--6) MeV. The mean free path for these neutrons is about 
1 that of the cosmic rays, so about & will escape. The upward flux is there- 
fore about .1 per cm? s near the poles, and about + of this near the equator, 
when we take into account the variation of primary flux with geomagnetic 


(9) R. McKeacue: Proc. Roy. Soc., A 236, 104 (1956). 


56 P. J. KELLOGG 


latitude. In addition, about .1 neutrons are produced in the backward di- 
rection in the energy range (30 +400) MeV. The mean free path of these neutrons 
is about the same as that of the cosmic rays, so 3 will escape, giving a flux 
of .05 per cm? s near the poles, and + of this near the equator. Clearly these 
estimates are not as accurate as those for the diffusing neutrons. These fluxes 
are therefore of the same order as the flux of diffusing neutrons. Fewer will 
decay in a given volume because of their greater velocity, but on the other 
hand, the lifetime of the resulting proton will be much longer because of its 
higher energy. The decay electron will have about the same energy as an 
electron from a neutron at rest over most of the energy range, so we do not 
have to consider these electrons further. 


3. — Estimate of particle densities. 


We now consider the decay of the neutrons after they have left the atmos- 
phere. As an approximation we take the neutrons to be travelling straight 
upwards and consider the diffusing neutrons first. 

IEEE 
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is the flux of neutrons leaving the earth with velocity between v and +de 
the number which decay between the distance Z and Z+dZ is given by 


FASI 
(20) CLN = rev VE) I exp per second, 


Z | dZ dv 
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where À is the neutron lifetime of 12 min. 


If we neglect the variation of L with energy, the maximum of the integrand 
in eq. (20) occurs at 


Aw EL AR 
9 hi / 
(21) vB = |va+3)/7| 


From previous discussion we know that l/Z is very slowly varying for 52 MeV 
so we will not make much error if we put 1/Z equal to its value at the maximum. 
The integral in eq. (20) may then be evaluated to give | 


(22) aN; (7) Lu 
"\ Ee} max (Z + (V2a/m) À) ? 


for the number of neutrons decaying per second between Z and Z-+dZ. 
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For a=.05 eV, we have V(2a/m)A=2160km. To get the rate of decay 
per cm? we must multiply by Rp ?, where À, is the radius of the earth’s 
atmosphere. 

The expression for the flux which we have used is based on slowing down 
by elastic collisions and, therefore, fails when the energy of the neutrons is 
greater than about 2 MeV. 2 MeV neutrons travel a mean distance of 
4-107!2 em before decaying. It will be shown later that this limitation is not 
important, but that other considerations probably limit the density of par- 
ticles at somewhat smaller distances than this. 

The neutrons which do not begin the diffusion process, but escape directly 
have a velocity high enough so that we may neglect the exponential factor 
coming from the decay of the neutron, therefore we have: 


dN az 
23 = I: 
( ) dt 1,2 20 b) 


where J,, I, are the fluxes estimated above. 

We now estimate the density of electrons stored in the radiation zone. 
In the lower part of the zone where the density of air atoms determines the 
lifetime of the electrons against scattering, we may expect that the density 
of electrons is inversely proportional to the atmospheric density, and will 
therefore increase exponentially with height. This correspondence will be only 
approximate because the electrons will change altitude somewhat in their 
rapid precession in the earth’s actual field. If this effect proves to be cal- 
culable, then the electron density will be a very sensitive method of determining 
the atmospheric density and scale height at very great altitudes. At present, the 
highest direct measurements of atmospheric density have been made so far below 
1000 km that any value of atmospheric density could appear reasonable, so 
it is possible to calculate lifetimes only in the region above the atmosphere, 
where the lifetimes in Table I are valid. Close to the earth, the density of 
electrons will be given approximately by the product of rate of formation, 
lifetime, and an average fraction of electrons which go into trapped orbits. 
From Table II we see that even for electrons with r,/R,— 1.05, cos 64° = .43 
of ees emitted electrons go into trapped orbits, so we may take some- 
thing like 3 for an average for this fraction. For the rate of formation we 
use eq. (22) neglecting + in the denominator and, using a lifetime of 3 :10° s, 


we obtain: 
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for the density of electrons, where we have used the equatorial value of 
neutron flux. 

Van Allen measures a counting rate in excess of 15 000 counts/s and esti- 
mates that the dead time correction would give a counting rate of more than 
35 000 counts/s in the counter of 20 em? cross-section. We may estimate roughly 
that electrons from neutron decays, therefore approximately .5 MeV electrons, 
will produce about .02 photons per electron in the 1 g/cm? steel shell of the 
satellite. Inside the satellite the X-ray flux ought to be fairly uniform and 
equal to twice the flux of X-rays out of an element of the surface. Using the 
efficiency of .3% which VAN ALLEN gives for the counter in this X-rays range, 
we obtain a flux of 


Sr 


35000 1 photons 1 electrons 
Hy aol ia Se AO: electrons/cm?s , 
20 .003 cem?s .02 photon 2 


~~ 


corresponding to an electron density of 1.5-10’/e=5-10-4/em® or greater. 


We conclude that such neutrons are capable of supplying a sufficient number 
of electrons. 

We now make a rough estimate of the density of electrons far from the 
earth. These electrons will follow paths reaching far out at the equator and 
coming near the earth at the poles, so we must take the variation of dN/dt 
with position into account. On the other hand, particles with small pitch 
angles will still be trapped so we assume all emitted particles are trapped. 
Electrons with guiding centers which cross the equator between 7, and 7,+dr, 
will be those emitted into a volume of revolution generated by the lines of B 
passing through 7, and r,<+dr,. An element of this volume may be written 
2ar sin 0D di, where dl — r d6(B/B,) = dr(B/B,) is an element of length along 
the line of magnetic field, and D is the perpendicular distance between lines. 
The flux between lines is constant so we may write 


(24) 277 sin O.DB = 277, dr Bo 
The total volume is therefore: 
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We now compute the number of electrons emitted into this volume. We 
multiply dN /dt from equation (22) by (R,/r)? to get the volume source density 
for electrons. Neglecting Vx in the denominator of eq. (22) and multiplying 
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by the lifetime to get the density of electrons, we obtain 
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where we have used the neutron flux measured at St. Paul as representing 
an average figure for the earth. 

This equation refers to the density everywhere along a line of B crossing 
the equatorial plane at ». Hence above the atmosphere at a constant height 
we should expect the electron density to show a strong latitude dependence. 
Writing r=r, sin?0 we have, away from the equator 


Rs, sine 0 
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This equation is of doubtful accuracy and we shall show later that there are 
probably even fewer particles near the poles. Nevertheless, it is clear that 
the number of stored particles must be very small near the poles. 

We estimate the density of protons from fast neutrons in the same way. 
In this case (eq. (23)) the production rate varies only as 7-2 so we obtain 


_ Re, Jrarvi em) 
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(28) 0 


The lifetimes for these particles will not be given by eq. (8) because nuclear 
scattering dominates the Coulomb scattering for protons of these energies. 
Over the important part of the energy range under discussion, the cross-section 
is approximately 3/H barns, where Æ is in MeV. This leads to lifetimes of 
5-1021$ for 5 MeV protons and 2-10! s for 50 MeV protons, which energies 
we take as representative of the two groups. Then eq. (28) becomes, for 
(R.,/1o) A 


WE (sal) for 5 MeV protons, 
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.02 (RxTo)° for 50 MeV protons. 
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The actual energy spectrum, of course, will be the same as that from stars, 
but slightly distorted by the dependence of lifetime on energy. The density 
of particles at low altitudes near the equator may be estimated as before and 
the results are collected in Table ITI. 


IAB ETS 


| «5 MeV» | «50 MeV» 


Electrons Total 
| | Protons | Protons 
VER È : Mat; Pa ETA A 
i Particle density 6: 10-3 cm? dello | EURE — | 
i Particle flux 2-108 em? s 6-108 107 == 
Energy flux 100 cies || 5 0 i 800 | 103 


We have already seen that the electron density is sufficient to give a 
counting rate a few times higher than is observed. The «5 MeV » protons will 
not have’ sufficient energy to penetrate the walls of the Geiger counter and 
satellite, but most of the «50 MeV » protons can do so, and would therefore 
give a counting rate about 5000 times higher than what is observed. The 
observations give a lower limit only, of course, but one would guess that col- 
lective processes are importabt in removing particles from storage, so that 
the lifetimes which we have used are too long. 


4, — Certain collective effects. 


We would now like to show that eqs. (26) and (28) are not likely to be 
valid for very large radii, even if the lifetimes we have used turn out to be 
correct, due to plasma instability. In spiralling around the lines of magnetic 
force, the plasma particles have magnetic moments which tend to reduce the 
field in the interiors of their orbits; therefore, the plasma has a tendency to 
be diamagnetic. This tendency is often counteracted by the behaviour near 
the edges of the plasma. Nevertheless, we shall consider only the diamagnetic 
effects here. Then, inside the plasma, an expression like: 


masi 
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will hold where M is the magnetic moment density, B, is the applied field, 
Dì is the average square of the velocity component parallel to the field, and 
o is the density of particles. 


GO 
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Solving for B, we obtain 
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so that we see that we have no solution when the energy density of the plasma 

is greater than one-half of the energy density of the applied magnetic field. 

If we consider a plot of B versus 0, we see that near the limit 

(31) ERO nes 
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the addition of à few particles causes à very great change in the magnetic field 
and we may expect that somehow instabilities develop which prevent the 
limit, eq. (31), being exceeded by throwing particles out of the region. This 
argument is perhaps not completely convincing, because we neglected the 
effects at the edge of the plasma but what we are saying is merely that a 
“plasma probably cannot be confined by a magnetic field whose energy density 
is smaller than that of the plasma. 

The limit (eq. (31)) is exceeded at about 8 earth radii. Beyond 8 earth 
radii, then, we may expect that the energy density falls off at least as fast 
as 1/r6. We should also expect rather complicated and interesting effects in 
this region. We note that the line of B which crosses the equator at 8 earth 
radii reaches the earth at a geomagnetic latitude which is rather close to the 
auroral zone. It may be that the aurora is a manifestation of the instability 
which we expect in this region, and that the observed correlation with solar 
phenomena is due to some triggering process. 

In this section we give a more accurate treatment of the density of stored 
particles. Let n(7o, Yo, u) Sin yo dp de du be the number of particles with 
guiding centers which cross the equatorial plane between 7, and r,+dr, and 
at the equator have pitch angle between y, and y,+dy, and have magnetic 
moment u. We write the rate of change of n as 


a dn (4 (= 
( 2) dt wt dii dt scat 


where dn/dt.,.. is the rate of increase of n due to decays of neutrons and 
dn/dt,,,, is the rate of decrease of n due to scattering of the particles. I the 
density of collision centers were uniform we should write 


(33) CR = novN,, 
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where N, is the number of collision centers per cm’. We expect that the number 
of collision centers will be a function of 7, so we need to average in the fol- 
lowing way, where the integral is taken over the orbit of the particles in 
question from the equator to the turning point 


(34) D) = ner y Pel 


nn i 


where v, is the radial velocity of the guiding center; v,=v cosy(B,/B)(dn/dt) ice 
is somewhat more complicated. Let (AN/A7) be the number of particles pro- 
duced per cm*/s by neutron decays. It may be a function of r and 0. We 
assume these particles are produced isotropically, so (AN/d7') siny dy r, dr sin 0 dé 
is the number of particles produced with pitch angle between y and y+dy 
in the volume element r? dr in 0 d0. Particles with equatorial co-ordinates 7, 
to r,+dr, and y, to w+dyw are those emitted in a volume bounded by the 
lines of B through 7, and r,+dr,, and the turning points corresponding to 
Por To- 

We may express the pitch angle in terms of the pitch angle at the equator 
by means of equation (1), 


(B/B;) ) sin y Yo COS PodPo 
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Let D be the perpendicular distance between two lines of B which cross the 
equator at », and n+d", so D=dr(B/B,). Then, using eq. (24), we have 


(36) 2779 dry By = 2ar sin 0 DB = 2ar sin 6 drB, . 


We have, finally, for the rate of creation of particles 
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In this integral we should express 7’ in terms of 0 by r= » sin?6 
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and the integral is to be taken from along the fluxline between 


the turning 
points corresponding to ro, Wo. 
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At equilibrium we have dn/dt — 0 so that we may solve eqs. (34) and (37) 
for n to obtain 


dr/v, [dN B DO N 
JA a ov(dr/v,) ©) dt Bo 4/1 == (B/B;) sin? Wo 


(wood is the number cf particles of equatorial pitch angle wo which 
cross the equator between » and dr. To obtain the number of particles in 
an element of volume we must multiply by dg/27 and by dl/v cos yT where 
dl = rd6(B/B,) is an element of length along a line of force, v cosy is the ve- 
locity parallel to a line of force and T is the half-period, the time to travel 
from one turning point to another. We express 7,d7, in terms of r by eq. (24) 
and obtain 
B 7 sin 6 dr dé dg 
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for the number of particles having equatorial pitch angle y, in the volume 
element 7? sin 0 dr dû dp. To get the density of particles of all pitch angles, 
we multiply by 


B sin Yo COS Yodo 
Bo V1 — (B/B,) sin? y 


(41) Sin y dy = 


and integrate over pitch angle from zero to yo, where y. 18 the smallest 
pitch angle for which particles reach the volume element in question before 
being reflected. 

The half-period 7 is given by 
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The integral has the value (7/3)V2 & 1.48 for y, ® 2/2 and 2.74 for y, & 0 
so is not a strong function of wp. 

We now attempt a more accurate treatment of the storage problem using 
these equations. We consider only particles with 7, > R,, so these considerations 
do not apply close to the earth near the equator. We take the earth’s atmos- 
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phere to be a perfect absorber with a sharp upper boundary so that 
dr (dr LUE fOr Yo > Pemax 

(43) [x 3 | | | Po Panax 
Ur | CO for Po< Pemax 


where t is given in Table I and y... is the largest pitch angle which allows 
particles to enter the atmosphere. 

We evaluate approximately the integral in equation (39) in the extreme 
cases of large and small y,. For y, & 7/2, 0 will also not be much different 
from 2/2. In the integrand dN/dt ~ dN/dt(7,) so we have 

dn dN 
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where 7 is given by eq. (42). 
For small y, we write 
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where we have changed eq. (37) to make r the variable of integration. 

r, is the value which makes the square root vanish, and we are interested 
in the case 7;<7). Then the principal contribution comes from r ~ 7r,. We 
therefore neglect terms of order r./r, and put B/B,=1 and V4 — r}r, & 2, 
obtaining 
ay i 
dt 
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where we have put «= (r/2*,) sin #w. The upper limit can be made infinite 
without serious error and then the value of the integral is about .86. For 
comparison with eq. (44) we put (dN/dt)(R,)R? = (dN/dt)(n)r} and have 
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For 7, very large, the distribution at the equator will be strongly peaked near 
the critical value y,,. 
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At large radii then, near the equator, there are large numbers of particles 
moving nearly perpendicular to the equatorial plane, in both directions. Such 
a situation arises because of the 1/r° dependence of the source density which 
implies that most of the particles were created in trapped orbits near the earth, 
and so have small pitch angles when they cross the equator. Such a situation, 
with charged particles moving in opposite directions through each other, is 
unstable against plasma oscillations. But in this case it seems likely that the 
plasma oscillatorions act at least as much to increase the pitch angles as other- 
wise and so do not throw the particles out of the trapping zone. In any case, 
the particles will not return to near the earth, but will be thrown out, or con- 
fined nearer the equatorial plane. 


5. — Summary. 


The earth is a source of neutrons which are produced in the atmosphere 
by the cosmic rays. Neutrons which leave the atmosphere and decay will 
often leave the decay proton and electron in trapped orbits in the earth’s 
magnetic field. We compute an upper limit of the density of trapped par- 
ticles and show these may be the particles observed by VAN ALLEN and co- 
workers in satellites. Reasons are given for believing that there may be few 
or no particles near the magnetic poles, the boundary of the distribution oc- 
curring near the auroral zone. If the auroral zone is actually the boundary 
of the radiation zone, then it seems possible that the aurora is a manifestation 
of instabilities occuring at the outside of the zone, triggered in some way by 
solar phenomena. 


The author would like to thank T. GoLp, Harvard University, and PAMELA 
ROTHWELL, Imperial College, London, who are responsible for the central idea 
of this paper. His debt to them has been mentioned more fully in the text. 
He would also like to thank E. P. NEy, C. J. WADDINGTON and J. R. WINCKLER 
for helpful discussions, and the members of the cosmic ray group at the State 
University of Iowa, particularly C. E. MCILWAIN, for communicating preli- 
minary results in Explorer IV. 


Supplementary note. 


As this paper was being submitted for publication, the author received Phys. Rev. 
Letters, 1. no. 5 containing two letters from 8. F. SINGER dealing with the same subject. 
Although the present paper would have been somewhat different if Singer’s letters 
had been read earlier, still the author feels that the topics treated in the present paper 
are sufficiently different from those treated by Singer to warrant its publication. 
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Note added in proof. 


Recently results of observations in Explorer IV (satellite 1958 e) have become 
available (*). This satellite carries four different instrument for measuring radiation. 
A geiger counter shows that the flux of 50 MeV protons, at the maximum altitude 
attained, must be at least one thousand times less than that calculated in Table IIT. 
On the other hand, a scintillator measures an energy flux of 80 erg/em?/s of particles 
of range greater than 1 mg/em?. Since it is unlikely that the satellite is high enough 
so that the values of Table III really apply, and since it would in any case be difficult 
to account for the great reduction in proton flux with no reduction in electron flux, 
we must consider this energy flux too great to be accounted for on the present hypoth- 
esis (the relevant flux from Table III would then be the 100 erg/cm?/s due to electrons). 
It seems likely, therefore, that the penetrating (range greater than 5 g/cm?) part of 
the Van Allen radiation may be protons due to neutron decay but that some other 
source is necessary to account for the less penetrating part. 


(*) J. A. Van ALLEN, C. E. McILwaIN and G. H. Lupwic: submitted for publi- 
cation in Journ. of Geophys. Res. 


RAS tS) (INP) (>) 


Discutiamo la possibilità suggerita da P. RoTHwEeLL e T. GOLD in conversazioni 
con l’autore che la radiazione osservata da VAN ALLEN e collaboratori alle grandi alti- 
tudini sia dovuta agli elettroni e protoni di decadimento di neutroni prodotti nel- 
l’atmosfera terrestre dai raggi cosmici ed immagazzinati nel campo magnetico terrestre. 
Si presentano valutazioni dell’ordine di grandezza delle densità prevedibili. Si consi- 
derano solo le perdite per scattering. Assumendo una vita media di 3-10° s per scat- 
tering, si ottiene un limite superiore di 10? elettroni/em? vicino alla terra e .05(R,/r,)° elet- 
troni/em* a grandi distanze 7) nel piano equatoriale. La densità protonica a grande 
distanza è .03(k,/r,)° cm per una vita media di 10!? s. Se le accelerazioni nel plasma 
non sono rilevanti, lo spettro elettronico sarà approssimativamente uguale a quello 
del decadimento 8 dei neutroni. I protoni sono prodotti da neutroni veloci provenienti 
direttamente da stelle nucleari e il loro spettro sarà approssimativamente uguale a 
quello dei protoni provenienti dalle stelle. La densità elettronica è sufficiente a dare 
un tasso di conteggio alcune volte maggiore di quello osservato, mentre la densità 
protonica è sufficiente a dare un tasso di conteggio 101 volte maggiore del limite minimo 
osservato. Le vite medie delle particelle immagazzinate sono pertanto probabilmente 
molto minori di quelle risultanti dallo scattering. Si dovrebbe rilevare un forte effetto 
di latitudine che si stima grossolanamente proporzionale a così 2. Si danno ragioni 
che fanno ritenere che effetti collettivi riducano in vicinanza dei poli la densità al 
disotto di tale valore. 


(*) Traduzione a cura della Redazione. 
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Elastic Differential Cross-Section 
for 60 MeV K*-Scattering in Emulsion. 


T. G. Lim and J. P. VAN DER LINDEN 


Natuurkundig Laboratorium - Universiteit van Amsterdam, The Netherlands 


(ricevuto il 2 Ottobre 1958) 


Summary. — By the method of the phase shift optical model analysis, 
the differential cross-sections for 60 MeV K'-mesons elastically scat- 
tered by complex nuclei in Ilford G-5 nuclear emulsion are computed. 
The calculations are based on the assumption that the elastic scattering 
is due to a symmetrical square well potential of +13 MeV, superposed 
on the Coulomb potential. The resultant emulsion cross section do/dw=f(#) 
curve is obtained from the separate contributions of the cross sections 
for the K*-meson scattering by C, Ag, Br and O. Previously the resultant 
cross sections were computed by Costa and PATERGNANI (1) on the 
assumption that the K°-mesons were scattered by the hypothetical 
elements $}E and !!E, representing the average of the heavy and the 
light elements in the emulsion, respectively. It is found that the do/dw=f(#) 
curves obtained by their method and by the present ones do not coincide 
but agree fairly well. 


1. — Introduction. 


The successful application of the optical model of the nucleus to the analysis 
of the m-meson scattering (?) is well known. Assuming also the applicability 
of this optical model to the K*-scattering, the real part of the nuclear potential 
can be obtained by the analysis of the elastic K*-scattering against complex 
nuclei. From the analysis of the inelastic K*-meson scattering, it is found 


(1) G. Costa and G. PATERGNANI: Nuovo Cimento, 5, 448 (1957). 
(2) A. PEVSNER and J. RAINWATER: Phys. Rev., 100, 1431 (1955). 
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experimentally that the imaginary part is small, of about a few MeV (1°). 
The inclusion of the imaginary part of the nuclear potential has little influence 
on the differential cross-section in the small angle region (4. With other 
authors (25) we therefore only consider the real part of the complex potential 
in the calculations of the elastic differential cross-sections in this angle region. 
The real part of the nuclear potential is assumed to be represented by a con- 
stant potential throughout the nuclear volume. We further assume a homo- 
geneous charge distribution in the nucleus, so that the total potential for the 
elastic K'-scattering is envisaged as a Coulomb potential on which is super- 
posed a spherical and constant real potential. The interference of the electro- 
magnetic and the nuclear scattering is expected to be strong in the region of 
small scattering angles. A refinement of the shape of the nuclear potential, 
as e.g. the use of a more realistic one (*) influences the do/dw = f(Ÿ) curve 
most in the region of large angles. It may be noted that through the obser- 
vation of the scattering events just in this large angle region, one can try to 
determine the form of the potential. However, this is more meaningful, if a 
separate scattering experiment could be set up using only one element (7). 
So, as long as we try to fix the sign and to get an estimate of the average 
depth value of the real part of the nuclear potential in this region, the omis- 
sion of the imaginary part in the calculation seems to be justified. 


2. — Computation. 


The differential cross-section is computed by the method of partial waves. 
According to this procedure the differential cross-section in the center of mass 
system, due to a pure Coulomb field is (8:°) 


GI 4 a : ? 
(1) He =" 19 SP [— 46 In (sin? 40) + in + Qing] | . 


The differential cross-section due to a symmetrical square well potential only, 


À 


(5) G. Cocconr, G. PuPPI, G. QUARENI and A. STANGHELLINI: Nuovo Cimento, 
D 172 (1957). 

(4) G..IGo, D. G. RAVENHALL, J. J. Tiemann, W. W. CHuPP, G. GOLDHABER. 
. GOLDHABER, J. E. LANNUTTI and R. M. THALER: Phys. Rev., 109, 2133 (1958). 

(5) T. F. HoanG, M. F. KAPLON and R. Custer: Phys. Rev., 107, 1698 (1957). 

(9) J. E. LANNUTTI, S. GOLDHABER, G. GOLDHABER, W. W. CHuPP, S. GIAMBUZZI, 
C. MARCHI, G. QUARENI and A. WATAGHIN: Phys. Rev., 109, 2121 (1958). 

() L. T. Kerra, T. F. Kycra and L. van Rossum: Phys. Rev., 109, 1784 (1958). 

(5) N. F. Morr and H.S. W. Massey: The Theory of Atomic Collisions (Oxford, 1949) 

(9) W. Gorpon: Zeîts. f. Phys., 51, 180 (1928). 
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can be written as (1) 


d 1 ae (i 
(2) = E D (224 1) exp[id,] sin 6, P,(cos 9) 
A 1 


We assume the elastic K*-scattering to be due to a modified Coulomb field 
(viz. a field due to a homogeneous charge distribution on which a constant 
real potential is superposed), so that the differential cross-section becomes (1) 


do | S Di Lon : A 
(3) or “Ok aC fa el - 16 In (sin? 39) + ia + 2ino] + 
+5 Se + 1) exp[iô, + 2iy,] sin 6,P,(cos d) |, 
l 


where ¢ = Ze?/hv, Z: nuclear charge number, e: elementary charge, #: Planck’s 
constant, divided by 2x, v: relative velocity of the colliding particles, 9: scat- 
tering angle; 7,: phase difference between the undisturbed /-th partial wave 
and the corresponding wave, disturbed by a Coulomb field; 6,: phase diffe- 
rence between the /-th partial wave distorted both by Coulomb potential 
and the symmetrical square well potential, and the corresponding wave dis- 
torted by the Coulomb field only. 

With Costa and PATERGNANI (') the evaluation of 6, is carried out using 
the W.K.B. approximation (1:11). 

The composition of the G-5 nuclear emulsion at 58% humidity (1?) is given 
in Table I, in which the elements are arranged according to their relative 
abundance. 

in E 


16 14y7 32Q 126,6 
80 7N 16> 531 


ERA | 
129 107,9 79,9 a 
20 | 10849 | "99Br 


| | | 


-- 
7 Element las 


l'a 
| 
Tu 
| Pace tt cumin 
sei of atoms 3.2102) 1.3401 1.0221 1.0169 0.9373 0.3141 0.0131 0.0057 


022/em? | | | | | 


Scattering of K*-mesons by free protons are rare in emulsion and more- 
over these can be distinguished from those events which are due to collisions 
against complex nuclei. As we only consider K*-scattering by complex nuclei, 

(19) L. I. Scurr: Quantum mechanics (New York, 1949). 

(11) P. M. Morse and H. FesHBAcH: Methods of Theoretical Physics, part II 
(New York, 1953), p. 1101. 

(2) A. RosENFELD, M. Backus, J. FRIEDMAN, W. F. Fry, D. HASKIN, J. LACH, 
R. Lux, M. Orans, J. OREAR, E. SILVERSTEIN, W. SLATER, F. SOLMITZ, R. SWANSON 
and H. Tarr: How to devilup emulsion (Chicago, 1955), ch. IT. 
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according to Table I it is clear that a good approximation to the real situation 
is obtained by taking only into account the scattering by C-, Ag-, Br-, and 
O-nuclei, being the most abundant elements. We therefore compute the elastic 
differential cross-sections of the K*-mesons against these nuclei separately and 
call them (do/dm),, (do/do),,, (do/dw),, and (do/dw),, respectively. The re- 
sultant emulsion cross-section do/dw per nucleus is the mean of the separate 
cross-section and can be obtained by the relation 


oe To sli (E): 


where «; is the relative number of atoms of the i-th element à — C, Ag, Br, O. 
vie 
i 


Previously Costa and PATERGNANI (!) computed the differential cross- 
sections starting from the hypothetical elements A= 95, Z=42 (57%), re- 
presenting the average heavy element and A=14, Z=7 (43%), representing 
the average light element. As expression (3) is not a linear one—it namely 
depends in an intricate manner on the nuclear parameters—we compute the 

resultant differential cross-section 


10° : = a i 7 5 by means of (4). 
St | To have a comparison of the 
| asian noto Gomme! 
rp by the two mentioned methods 
10 +! CON in LÉ A ati are carri 
| È our computations are carried out 
Sy De for the values of the nuclear po- 
n tential and K*-energy, also used 
10 by CostA and PATERGNANI (1). 
sl We chose for the depth of the 
square well nuclear ‘potential: 
À ue MeV and for the K*-energy: 
10° + 60 MeV. 
5+ 
Fig. 1. — Differential cross section 
10 Fa mbe. ‘ per nucleus versus scattering angle, 
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From (4) it follows that the value of the resultant cross-section is strongly 


influenced by the greatest «,;(do/dm),;, i.e. 
high constituting differential cross-section. 
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Fig. 2. — Differential cross section per nucleus 

versus scattering angle for 60 MeV K*-mesons, 

elastically scattered by C-, Ag-, Br- and O- 
nuclei, respectively. V: +13 MeV. 


By way of example we varied the value 
of the Ag-radius from 6.3 fermis (from 
the neutron experiment) with its highest 
I-number= 5 to 9.0 fermis (from the 
x-experiment) with its highest /-num- 
ber=8. The corresponding (do/do) ,,=f(?) 
curves are drawn in Fig. 1. 

The effect of the variation of the Ag- 
radius is reflected in the variation of the 
resultant curves (Fig. 3). 

To demonstrate of which components 


by the most abundant element of 
This leads us to investigate how 
the differential cross-section do/dm 
varies, when the Ag-nuclear radius 
varies, keeping V: +13 MeV and 
the K*-energy: 60 MeV. In this 
way we can also obtain an idea 
of the spread of the theoretically 
expected value, as we in fact do 
not know the appropriate range 
of the nuclear potential. 


Per as 
| 


3 
10 rr ye. 


L 
| 
10° Scattering angle 

0° PP 160 AP 2 À 
Fig. 3.— Resultant curves do/do={(9). 
K+-energy: 60 MeV. V: +13 MeV. 
For curves I and II, the K*-mesons 
scatter against C-, Ag-, Br-, and O- 
nuclei; for curve I the Ag-radius is 
9.0 fermis and for curve II the Ag- 
radius is 6.3 fermis. Curve III is 
derived from the computations of 
Costa and Patergnani. 


the resultant do/dw = f(#) curve is built up, we draw as an example the com- 
ponents of curve I of Fig. 3. These (do/dw);— f(Ÿ) curves (î= €, Ag, Br, O) 
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are drawn separately in Fig. 2. The Ag-radius is 9.0 fermis with the highest 
l-number = 8; the radii of C, Br and O are taken from the neutron exper- 
iments and are 3.9 fermis, 6.0 fermis and 3.4 fermis respectively with their 
highest /-numbers 3, 6 and 3, respectively. The oscillatory character of the 
curves depends on the form of the nuclear potential, which is superposed on 
the Coulomb potential. This is most pronounced in the large angle region, 
whereas in the small angle region the oscillations are smoothed out. 

In Fig. 3 the resultant curves I, II and III of do/dw = f(Ÿ) are drawn for 
the small angle region. Curves I and II are obtained by considering the scat- 
tering of K+-mesons to take place against C-, Ag-, Br-, and O-nuclei; for 
curve I the radius of Ag is taken 9.0 fermis and for II it is taken 6.3 fermis. 
Curve III is derived from the computation of COSTA and PATERGNANI, where 
the K+-scatterings are assumed to take place against the hypothetical ele- 
ments A=95, Z=42 and A=14, Z=7 only (7). We note that curve III 
does not differ very much from the curves I and IT and that it lies within the 
spread of the resultant cross-sections with the corresponding Ag-radius 9.0 
and 6.3 fermis, respectively. 


This work is part of the Research Program of the Netherlands Institution 
of Fundamental Research of Matter (F.0.M.), financially supported by the 
Netherlands Foundation for Pure Scientific Research (Z.W.0.). Many thanks 
are due to Professor Dr. G. W. RATHENAU for his stimulating suggestions and 
criticisms. The valuable discussions with Mr. F. ARTMANN, Dr. F. BRUIN (*) 
and Dr. M. BRUIN (*) are here gratefully acknowledged. Miss C. J. VAN ALBADA 
and Mr. 8. BosGRA are thanked for their help in the tedious computations. 


(*) Now at the American University, Beirut. 


RIASSUNTO (*) 


Si calcolano col metodo dell’analisi degli spostamenti di fase nel modello ottico le 
sezioni d’urto differenziali pei mesoni K* da 60 MeV diffusi elasticamente dai nuclei 
complessi dell’emulsione nucleare Ilford G-5. I calcoli sono basati sull’assunzione che 
lo scattering elastico sia dovuto ad una buca di potenziale quadrata simmetrica di 
+13 MeV sovrapposta al potenziale coulombiano. La curva che rappresenta la sezione 
d’urto do/do=f(9) risultante per l’emulsione si ottiene dai contributi separati delle 
sezioni d’urto dello scattering del mesone K* su C, Ag, Br e O. Le sezioni d’urto risul- 
tanti erano state precedentemente calcolate da Costa e PATERGNANI nell'ipotesi che 
i mesoni K' fossero diffusi dagli elementi ipotetici DE e E, rappresentanti, rispetti- 
vamente, le medie degli elementi pesanti e leggeri dell’emulsione. Si trova che le 
curve do/do=f(9) ottenute col loro metodo e col presente non coincidono ma si accor- 
dano soddisfacentemente. 


(*) Traduzione a cura della Redazione. 


IL NUOVO CIMENTO Woo SIN dI 1° Gennaio 1959 


Three Body Forces in Hypernuclei (*). 


G. G. BACH 


Istituto di Pisica dell Università - Roma 


(ricevuto il 15 Ottobre 1958) 


Summary. — An investigation is made to ascertain the relative importance 
of the two-body and three-body A-nucleon potentials which have the 
same order of coupling constant. The three-body potential, which is 
derived using perturbation theory from a pseudoscalar interaction, is 
found to be weakly singular and hence, much smaller than the two-body 
potential. In addition, the largest terms are non-central. In the hyper- 
triton, where correlations are large, it is found that three-body forces 
are repulsive for S states and may amount to 4%. Suitable D state 
admixture could lead to a repulsive three-body force of nearly 6%. In 
heavy hypernuclei, the relative importance of three-body forces is about 3% 
if correlations are neglected, and may be as large as 73% il correlations 
are as strong as they are in the hypertriton. It is also found that by omit- 
ting the Feynman graphs with bare lines the three-body contributions 
are negligible. 


1. — It was first pointed out by DALITZ (1) that the hypernuclear force 
between a A-particle and a nucleon could be brought about through a double 
pion exchange. If this is the case, it then seems likely that the three-body 
hypernuclear potential may be important since it would have the same order 
of coupling constant as the two-body hypernuclear potential. Calculations 
have been made by SPITZER (?) on the relative importance of this three-body 
potential. 

In this paper, a quantitative estimate is made of the relative importance 
of the three-body force for the case of the hypertriton, where one would expect 


(*) This paper forms part of a thesis, which was presented at McGill University 
in April, 1958, in partial fulfilment of the requirements for the degree of Doctor of 
Philosophy. 

(1) R. H. Dazrrz: Phys. Rev., 99, 1475 (1955). 

(2) R. Spirzer: Phys. Rev., 110, 1190 (1958). 
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correlations to be quite strong, and in heavy hypernuclei, where one would 
expect correlations to be quite weak. In Sect. 2, the three body-potential is 
derived and examined for its relative magnitude. It is found that without 
correlations, the more singular parts vanish since they are non-central. A self 
consistant calculation is then carried out in Sect. 3 on the hypertriton to obtain 
a wave function with correlations as well as a two-body A-nucleon potential 
which will give a A separation energy (binding energy) near zero. Then, in 
Section 4, using this wave function and two-body potential as well as the 
three-body potential, which is derived in Section 2, a direct calculation of 
<V (3-body)>/<V(2-body)> is made. The sensitivity of <V(3-body)> with D 
state admixture is also calculated. Following these considerations, heavy 
hypernuclei are treated in Section 5. 


2. — If one assumes pseudoscalar coupling with the spins of all baryons $ 
and the A and & particles with the same relative parity, the lowest order 
terms of the interaction hamiltonian are 

9 ; 35 a) LV | 
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where we have set M.— M,. Here 2, À, y and o are the wave functions 
for the three baryons and the 7-meson. 

To evaluate the three-body potentials, we use ordinary perturbation theory 
and assume that the baryon wave functions are 6 functions. In the energy 
denominators, we neglect the quantity 4 = M.,— M, in comparison with o,, 
the energy of the meson, in the intermediate states corresponding to graphs 
in Fig. la. The calculation of the three-body potentials is then straight- 
forward and the results for the three type of graphs in Fig. 1 are as follows: 
From Fig. la: 
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In these formulas, units of pion Compton wavelength are 
used and S,,, and S,,, are the usual tensor operators formed oat 
with the vector x or y. | ; 
Fig. lc 
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The factor 2 multiplies V, in accordance with BRUECKNER and WATSON (?) 
who introduced this damping factor for pair terms. 

It is interesting to notice that the more singular parts of the potentials 
are non-central and with no correlations between spin and direction or between 
directions, only the weakly singular central parts survive. Also V, then aver- 
ages to zero. Comparing V,, and V,, we see that V, is larger than V,, by 
a factor of approximately zu/A= 6. To get orders of magnitude from these 
potentials then, we examine V, with no correlations and consider the two 
nucleons in a spatially symmetric state so that <(t,-T,)(6,;-6,)> =—3. Then 
setting (92/47)(u*/4M?)= 0.08 and g,=g, we obtain 


9— 2—Y 


(5) V, (no correlations) = 0.74 — 
ry 


MeV. 


Examining equation (5), we see that V, with no correlations is very small 
and can be neglected in comparison with the two-body potential. Also, V,,, 


(3) K. A. BrueckNER and K. W. Watson: Phys. Rev., 92, 1023 (1953). 
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which is about } 
‘annot be reconciled with Spitzer’s (2) where V,, was calculated and found 
to have a comparable magnitude with the two-body potential when there 


were no correlations. 


3. — The form of the three-body potentials derived in Section 2 show that 
without correlations in the wave function of the hypernucleus, one would not 
expect that they would be very important. However, since the non-central 
parts are more singular, it is worth-while to consider a particular hypernucleus 
in some detail to see if correlations can make three-body forces important. 
The simplest hyperfragment which can be considered for this purpose is the 
hypertriton. A complete calculation on any heavier hyperfragment would be 
much more complex. 

To calculate the expectation value of the three-body A-nucleon potential 
for the hypertriton, a wave function with correlations is required as well as 
a knowledge of the core, e, of this three-body potential. It is also necessary 
to have a A-nucleon two-body potential so that one may calculate the ratio 
<V (3-body)>/<V(2-body)>. Finally, one requires that the value of B,, the 
A separation energy, be near zero. 

We therefore consider first a trial wave function 


(6) Ws= As EXP = (x + y) Pe 


9 ? 
where we let x be the variational parameter and choose f — 1.4 inverse meson 
units, 


(7) Xs = — Vaa(l)a(2)B(A) + VG (a(1)B(2) + B(1)a(2)) (A), 


N, which is a spin wave function with J=4 and M ee 


This choice of x, is made so that the two nucleons 
are in a mutual triplet spin state and also in accor- 


x dance with evidence that the A-nucleon interaction 
e is more attractive in the singlet state (4). 
x We next may write the hypertriton hamiltonian 
= (omitting three-body forces) as 
Na 
Fig. 2. (3) He ABC VW + V(A, N) + LAVO N;) ) 


where T is the total kinetic energy of the system excluding the motion of the 
center of mass, V, is the two-body nucleon-nucleon potential and V(A, N,), 


(4) R. H. Darirz and B. W. Downs: Phys. Rev., 111, 967 (1958). 
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of this can certainly be neglected for this case. This result 
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V(A, N,) are the A-nucleon two-body potentials. If one calculates the expec- 
tation value of H and adds to it the binding energy of the deuteron (2.226 MeV), 
one obtains approximately the binding energy B, of the A-particle. We choose 
a form of the A-nucleon two-body potential such that it requires a core e’. 
We then write 


(9) B,(a, €') = <H) + 2.226 MeV 

and setting 0B,/ca=0 when B,=0 determines an « and a core e’. 
need a core e for the A-nucleon three-body potential for the calculation of 
<V(3-body)>/<V(2-body)>, the A-nucleon two-body potential is chosen which 
‘an be derived in a similar manner as the three-body potentials so that we 
can assume e= e’. 

DALLAPORTA and FERRARI (5) and LICHTENBERG and Ross (5) have de- 
rived A-nucleon two-body potentials in a similar manner as the three-body 
potentials are derived here. DF considered only graphs where there were no 
bare lines while LR considered only bare line graphs. Bare lines (refers to the 
nucleon lines) means a graph where an intermediate state exists with no 
mesons present; in other words a graph where no meson lines are allowed 
to cross. Of these two potentials, we note that the one derived by DF is 
more attractive in the singlet state. This potential is given by equation (5) 
in that paper. We therefore select the one derived by DF. 

V,, the nucleon-nucleon potential, is taken to be (?) 


Since we 


e Te : eve À 
PET: — 9 V: : > .283 
(10) | vg ee co, Qe 
| 0 02590 


where V,=100.7 MeV, V,= 257 MeV and y and y’ are 1.735 and 1.696 in- 
verse meson units. We see that equation (10) differs from that mentioned in 
reference 16 in the paper by GAMMEL and THALER (’) in that the latter has 
a hard core at 90 =.283. It is modified by putting V, = 0 so that it will be 
more suitable for the wave function y,. The most suitable wave function 
would be zero at the hard core. This is further discussed in the conclusions. 
Equation (9) is plotted in Fig. 3. It is seen that the minimum for B, with 
respect to « is very close to zero when ¢=.388. From this, we take 


{ «= .386 meson units 
(11) i . I 
| nm = LS inverse meson units 


(5) N. Datuaporta and F. FERRARI: Nuovo Cimento, 5, 111 (1957). 
(6) D. B. LicarenBERG and M. Ross: Phys. Rev., 103, 1131 (1956). 


no 


(7) J. L. GammeL and R. M. THALER: Phys. Rev., 107, 1337 (1957). 
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and using these values of « and e, we get 


| TR 
(12) | <V,> =— 14.2 MeV, 
| <V(AN>)+ V(AN:)) =—46 MeV. 


It is interesting to notice that the size of the parameter «= 1.8 is larger 
than the deuteron parameter (B— 1.4). This may seem objectionable, since 
it infers that the A-nucleon parts 
of the wave function are more 


concentrated at the origin than 
the deuteron part is, which ap- 
pears to contradict the assump- 
tion that the A-particle is loosely 
bound. To explain this, we note 
that the true A-nucleon and deu- 
teron parts of the wave function 
are themselves composed of two 
parts; one is in the region of the 
potential and the other is in the 
-30 tail region. The latter is the part 
on which depends how loosely 
bound the particles are. The trial 
wave function given by equation 


B, MeV > 


50 
(6) corresponds more to the part 
-60 in the region of the potential since 
di See it is there that the parameters are 
Mics 13: 


determined. The value obtained 

for « thus reflects the fact that the 
effective range of the A-nucleon potential is smaller than the nucleon-nu- 
cleon potential. 


4. — It is of interest to examine the sensitivity of <V(3-body)) with D 
state admixture in the wave function. We rewrite the wave function for the 
hypertriton as 

Yo 
VN,” 


where y, is given by equation (6), and N, and N, are normalization constants. 
y,, the D state admixture is assumed to be 


(13) yp = cos = + sin g 
à 


(14) Yo = Xp EXP 
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where 

R Mt 12 K + 428 4 
(15) Xp cuni (w Dive Ì Y DET 


and the percentage of D state is given by 100 sin? g. The D state part of 
the spin wave function is chosen as shown since it is felt that the one spatially 
symmetric between the two nucleons will give the largest contribution to the 
binding energy. 


We now calculate <V(3-body)) using this wave function obtaining 


(16 CV.) = .108 cos? m + .119 sing cos mw —.106 sinz@ MeV, 
NB (2 P (2 (2 


(17) Va + Vp = 2.00 cos? p— 3.09 sin y cos g —.719 sin MeV, 
where V,, and V, are given by equations (2) and (3). Note that <V,>=0 
for this case, 

{Vy)/<V(2-body)> (no bare lines) 
and 

Vie t+ V,>/<V(2-body)> (bare lines: included) 

are plotted against percent D state in Fig. 4 <V(2-body)> is taken to be 
— 46 MeV (equation (12)) so that the ratios are valid for small g. For the 


T— T T T iS T T T | 
5h JI 
Ji | 
al tive È 

2 


0 NE ON en {= 


a0 AO GO CC, 70 ee 
Positive D Percent D state 


<V(3-Body)> /<v(2-Body)> * 100 
L 
L 


ij 
N 


larger values of y, the ratios only show the sensitivity of the three-body force 
with D state. 
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5. — The results of Section 4 (see Fig. 4) lead us to expect that three-body 
forces may be as large as 6% and repulsive when compared to the two body 
A-nucleon forces for the case of the hypertriton. In a heavy hypernucleus 
with one A-particle present, one expects that the three-body A-nuclear forces 
would be more important since the A-particle will not be governed by the 
Pauli principle and could be in the vicinity of a large number of nucleons. 
This might lead one to expect that the relative importance of three-body forces 
would be proportional to A with large A. 

If we again look at the heavy hypernucleus problem from another point 
of view, we can see that three-body A-nuclear forces may not be important 
since the more singular non-central parts of the potential would tendto average 
out to zero coherently when the A-particle interacts with many pairs of neigh- 
boring nucleons. In other words, one would expect that the effect of cor- 
relations would be small. 

The following consideration of the three-body problem follows in many 
respects the work reported by DRELL and HUANG (8). We consider a large 
hypernucleus consisting of A nucleons and one A-particle contained in a vo- 
lume v. We write the volume as 


7 3 
(18) n di Lo 


where u!= 1.42-10-1% cm, which is a meson Compton wavelength, and 7 is 
a parameter which determines the nuclear density. 

We then assume that the wave function for the A nucleons is given by 
the Slater determinant 


| fey... ot) 
(19) US 


Wye ea UE 
pal) DOLO Pu(A) 


where ;(j) = (1/40) exp [ik,;-r;}y;(0;)v;(t;) which are all mutually orthogonal. 
zi(0;) and »,(t;) are the spin and isotopic spin variables assigning the i-th spin 
and isotopic spin state to the j-th nucleon in the spatial state k,. 

Since the A-particle is not governed by the Pauli principle here, we write 
the total wave function for the system as 


(20) P=y-y(A), 
where we put 


(21) OA Ve ein [i VE 
V/V EN NOH 


(8) S. D. DreLL and K. Huane: Phys. Rev., 91, 1527 (1953). 
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For the A-nucleon two-body potential, we again use the one derived by DF. 
The expectation value is 


(22) ‘V(2-body)> = > [Oort bares, 


t=1 


Notice that the (6,°6,) and tensor terms of V,,. vanish so that only the spin 
independent part need be evaluated. We therefore get 


A = — ee 
V (2-body)> = > {ety i ( 1) V; DEF (1A); (1 )p(A) dra Cha 5 


where the factor 1/4! comes from the normalization and (A —1)! from the 
fact that every term in the sum contains (A — 1)! permutations of the re- 
maining A — 1 factors which integrate to unity. 

Therefore, since 


AT 
| <V (2-body)> = È | Vor(1A) ar dry 5 
(23) LE 


| = - | "pr(14) dt 


In the calculation of the three-body expectation values, only the (6,- 6,)(t, : T.) 
term survives so that we get 


ED body) = ¥ [OAI 0 de, 


A\(A=2) ENT. 
E fig er Gt n) 


-#(4)(@;:(1)p;(2) —- p,(1)@:(2)) dt: AT: ÂTA . 


Since the three-body potentials do not contain the spin of the A-particle, we 
can put |y(A)|?=1/e. 

The spin-isotopic spin part is done as follows: It has been assumed that 
the system consists of 4/2 protons and 4/2 neutrons each having equal num- 
bers of spin up and down. Therefore, 5/8 of the nucleon pairs are anti- 
symmetric and 3/8 are symmetric in interchange of their spatial co-ordinates. 
The table below gives the values for {(o: :6:)(t,: t,)> for the various possibilities. 


Space Spin triplet 1 Spin singlet 

antisymmetric I-Spin triplet I-Spin singlet 
Space Spin triplet 3 Spin Singlet _ 

symmetric I-Spin singlet I-Spin Triplet 
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From this, equation (24) becomes 


OT | , 
(25) <V(8-body)>=— > X jexpli(ki— kj):(ri=r)]V024) dr, dz, dt, ; 

OÙ i,)=1 

ORAL 

ds v3 m 124)dr dr, dre ; 
where 

By Ie 
(26) D?(k, ne) Tg FE fe i | b) 
ml12 


and k,, —1.52u/n or the largest momentum of the highest filled level for the 
free particle states in the nuclear well. Here, we use the same notation as 
DRELL and HUANG (8). 

From equations (23) and (25), <V(3-body)>/<V(2-body)> can be evaluated. 
Using the three-body potential with no bare lines or V,, (equation (2)), this 
ratio is 
JD? (Ao) (e-*/ay) (ky (y) — y ko(y)) dx d'y 


yu ZE 


~ 120 [as (4 + 4a + 2?) (Fy (a) /a) + (2+ Qe+ 2°) ky (a) | dx ? 


x>e 


and the result, using the potential V, (equation (3)), is 


[D*(Ao)(e-*-¥/ary) dx doy 


mw ob À rue 
2404107 f(e*/a?)f (4+4e+¢@ £2) (Ke, (a ae) fac) + (2 + Qa + x?) (a r)] déxr 
2>e 


where we have changed to meson units in equations (27) and (28) and put 
A=Kk,|u. Since |R,| will be larger than R,,, we evaluate R,. The nume- 
rator of equation (28) becomes, apart from the factor — (x/24)(u/1A)(4/v), 


a Cea (ae ee 
LE | j1(40) oy d°x dey = — | Îi(40) (o — 2e)do, 


This integral and the denominator of equation (28) were evaluated nume- 
rically for various values of 7, using «= 0.386 (equation (11)). The results 
are shown in the table below, 

n 0.8 1.0 1.2 


Ri, — 0.005 1 — 0.0047 — 0.0043 


oF 
D 
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Since R,,, which is given by equation (27), is of opposite sign, it will tend 
to reduce the ratio if included. However, the calculations in Section 4 (see 
equations (16) and (17)) show that R,, will be much smaller than his Ae 

It is interesting to notice here that the ratio R, does not vary as A as one 
would expect it would, but rather increases in a complicated manner as the 
nuclear density increases. The reason for this can be clearly seen as follows: 
First consider the A-particle and a particular nucleon N,. The two-body 
potential will be described by V(AN,) and the three-body potential with these 


A 
two baryons will be > V(NN;). If we now assume that A and v are infinite, 
i =2 


a measure of the relative importance of the three-body to the two-body po- 
tential is 


which we see is finite since the range of the three-body force is finite. In other 
words, the numerator of equation (29) will only consist of a finite number of 
terms, the number of which will depend on the nuclear density and the range 
of the A-nucleon three-body potential. 

The above table therefore shows that three-body forces are not likely to 
contribute more than about 3% in heavy hypernuclei if the correlations are 
not present. From that table, a qualitative upper estimate of the size of 
three-body forces in heavy hypernuclei can be obtained if correlations are 
present to the extent that they are in the hypertriton. This can be done by 
estimating the values of 7, corresponding to a heavy hypernucleus and the 
hypertriton, and then by obtaining values of À, from the table for these values 
of 7. The ratio of these values of R, multiplied by 6%, which is an upper 
estimate of the relative importance of three-body forces in the hypertriton, 
gives then a qualitative estimate of the size of three-body forces in heavy 
hypernuclei if correlations are as strong as they are in the hypertriton. 

For heavy nuclei and the deuteron, it is estimated that 7 is 0.94 and 1.48 
respectively. We therefore get 


R,(n = 0.94) = — 0.0048, 
R,(n = 1.48) =— 0.0038 . 


From this, we conclude that with correlations in heavy hypernuclei, the rela- 
tive importance of three body forces may be as large as 75% and repulsive. 


6. — It must be admitted that there are many weak points in the foregoing 
derivations and subsequent calculations. To begin with, if the spins of A and X 
were not } or if their relative parities were different, the form of the inter- 
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action hamiltonian would change. Also, if the K-meson plays the more im- 
portant role in hypernuclear forces, we would expect that three-body forces 
would be about as important as they are in ordinary nuclei. Another question 
is the use of perturbation theory. Furthermore, approximations were made 
in the perturbation theory, the most important of which was the assumption 
that the baryon wave functions were 0-functions and neglecting the recoil. 
These assumptions will be in error for mesons of large momentum. On the 
other hand, disregarding the mass difference 4 = M,— M, in the energy de- 
nominators, introduced errors for mesons of small momentum. If we examine 
all of these items, it seems likely that the most critical factor is the form of 
the interaction hamiltonian, since once it has been chosen the qualitative 
form of the potentials is fixed. 

The next item open to criticism is the method of determining the wave 
function for the hypertriton and core e associated with the A-nucleon three- 
body potential. To begin with, the choice of the two-body potential V,, may 
be criticized. Also the trial wave function (equation (6)) is not suitable for 
use with the phenomenological nucleon-nucleon potential since the deuteron 
part of the wave function does not vanish at the hard core. In addition, the 
value chosen for 6 may be poor. 

To make the comparison of the two and three-body forces, one must have 
a two-body potential which has all of the desirable features; i.e. it should be 
more attractive in the singlet state. Further, since the core e, associated with 
the three-body potential is desired, a two-body A-nucleon potential is re- 
quired which is derived in a similar manner as the three-body potential; i.e. it 
also must have a core associated with it which, for consistency, we assume is 
equal to the core associated with the three-body potential. These conside- 
rations lead to the choice of V,p. 

As for the choice of the trial wave function, we see that since it does not 
vanish at the hard core, errors are introduced into the value of <H). It tends 
to make the kinetic energy term too small since the trial wave function is too 
smooth near the core. Again, since the trial wave function is finite and not 
zero at the core boundary, it tends to make the value of |< V(NLN,.)>| too 
large. Both of these errors which are in the same direction are estimated to 
make <H) too small by about 4.5 MeV. If we examine the graph in piges, 
it is seen that x would not be greatly changed since the minimum is quite 
shallow. e, the A-nucleon core radius, will be reduced somewhat, approximately 
from 0.386 to 0.380 meson units. 

From the foregoing remarks, it is of interest to estimate the changes in 

V(3-body)> with small variations in the quantities «, B and e. By writing 


E (x, B,e) = <V(3-body)), we have for small variations in the parameters 
i a i OE, < OA S oH 
(30) Fac: B, e) = Ey (a, Bo, 80) + a Om + Fe SI dp + = A Se ; 
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where %, f and e, equal 1.8, 1.4 and 0.386 and the derivatives are evaluated 
at these values of the parameters. If we consider only the S state trial wave 
funetion (equation (6)) and the three-body potential V, (equation (2)), the 
result is 


(31) Ey y, = 0.108 + 0.205 da + 0.096 38 — 0.977 de MeV 
and when bare line graphs are included (equations (2) and (3)), we get 
(32) yt Enng= 2.00 + 3.610x + .2.4936 — 8.463e MeV. 


The variations expected in the parameters are small. For example, a better 
value for e may be 0.33 and perhaps / should be nearer 1.52 say. If « remains 
unchanged, we see from equation (32) that the three-body expectation value 
would go from 2 to 2.35 MeV which is a change of 17%. It is further expected 
that the two-body expectation value would also vary so that the ratio 
<V(3-body)>/<V(2-body)> would not be greatly affected. 

In Section 5, three-body forces are considered in heavy hypernuclei. An 
obvious assumption made there is that only one A-particle is present and 
effects due to the heavier X-particle are ignored. The -particle probably 
plays an important role in heavy hypernuclei because the binding energies 
involved become more comparable to A= M,— M,. However, since a com- 
plete calculation of that nature was not within the scope of this work, a qualil- 
tative estimate of the importance of three-body forces was made considering 
only one A-particle in a large nucleus. 

We may summarize the conclusions as follows: 


a) Since the A-nucleon potential can be obtained by a consideration of 
a double pion exchange, one would expect that the two and three-body 
A-nucleon potentials, which have the same order of coupling constant, would 
have comparable magnitudes. It is found, however, that the three-body po- 
tential is much smaller than the two-body potential due to its weak singu- 
larity. 


b) The largest terms in the A-nucleon three-body potential are non- 
central whereas the central part is small. 


c) For the hypertriton, where correlations are expected to be large, 
three-body forces are repulsive for S states and have a relative importance 
of about 4%. Suitable D state admixtures could make them further repulsive 
and perhaps large as 6%. 

d) In heavy hypernuclei, two factors play an important part in deter- 
mining the magnitude of three-body A-nucleon forces. The first factor is the 
increased nuclear density which enhances the effect of the three-body force 
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and the second factor is the reduced amount of correlations in the nucleus. 
The latter effect is the more important since the most singular parts of the 
A-nucleon three-body potential are non-central. It is expected that with no 
correlations, the relative importance of three-body to two-body A-nucleon 
forces is only about 4%. With correlations as strong as in the hypertriton, 


it may be as large as 75%. In either case the contribution is repulsive. 


e) If the three-body potential can only be represented by graphs with 
no bare lines, then the relative importance of three to two-body forces in 
hypernuclei is less than $°% even with correlations present. 
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and rendering valuable assistance in the preliminary stages thereof. Special 
thanks are due to Dr. E. L. Lomon who, in Dr. SHARP’s absence, directed the 
major part of the work. Finally, financial assistance in the form of Studentships 
for the years 1956-57 and 1957-58, given by the Canadian National Research 
Council, is gratefully acknowledged. 


RIASSUNTO 


Si studia l’importanza relativa che hanno le forze tra due e tre corpi nel poten- 
ziale A-nucleone. Il potenziale proveniente dalle forze tra tre corpi, è ottenuto con 
un calcolo perturbativo con interazione pseudoscalare all’ordine più basso. Si trova 
che esso è poco singolare e quindi assai più debole di quello generato dalle forze tra 
due corpi. Nell’ipertritio, dove le correlazioni sono forti, si trova che le forze tra tre 
corpi sono repulsive per stati S e possono dare fino il 4% del totale delle forze; una 
conveniente introduzione di stato D può portare a forze repulsive del 6%. Negli iper- 
nuclei pesanti, l’importanza relativa delle forze tra tre corpi è di circa 3% se si tra- 
scurano le correlazioni e può giungere fino a 7.5% se le correlazioni sono forti come 
nell'ipertritio. Si trova anche che le forze tra tre corpi sono trascurabili e si omettono 
i grafici di Feynman che contengono stati intermedi senza mesoni x. 
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On the Renormalization of a Parity Non-Conserving Interaction. 


K. SEKINE 


Department of Physics, University of Tokyo - Tokyo 


(ricevuto il 10 Novembre 1958) 


Summary. An exact formulation is given of the renormalization method 
in a parity non-conserving theory. Equations determining the renor- 
malization constants are obtained and, among other things, the internal 
consistency of the theory is examined. 


Introduction. 


Although the parity non-conservation is well established experimentally (*), 
it does not seem that its theoretical implications have been fully exhausted. 

As is well known, the present quantum theory of fields is based on 
certain general assumptions such as Lorentz invariance, the existence of the 
energy-momentum four vector, asymptotic conditions, and so on. So far the 
invariance of the theory under space inversion has been one of them. However, 
now that the parity non-conservation is established experimentally, the in- 
variance mentioned must be ruled out from the number of fundamental postu- 
lates of quantum field theory. 

In this way, a problem is to be put forward to reformulate the theory without 
assuming its invariance under space inversion. This problem was partly taken 
up by SOLOVIEV, GUPTA and Iso (?), who have attempted to classify inter- 


(1) T. D. Lez and C. N. YANG: Phys. Rev., 104, 254 (1957); ©. S. Wu et al.: Phys. 
Rev., 105, 1413 (1957); 106, 1361 (1957); R. L. Garwin, L. M. LEDERMAN and M. 
WEINRICH: Phys. Rev., 105, 1415 (1957); J. I. FRIEDMAN and V. L. TeLEGDI: Phys. 
Rev., 105, 1681 (1957); F. S. Crawrorp, M. Cresti, M. L. Goop, K. GOTTSTEIN, 
E. M. Lyman, F. T. SoLmirz, M. L. Stevenson and H. K. TicHo: Phys. Rev., 108, 
1102 (1957); F. Ester, R. PLANO, A. PRopELL, N. Santos, M. SCHWARTZ, J. STEIN- 
BERGER, P. Bassi, V. BoreLLI, G. Puppr, G. TANAKA, P. WALOSCHEK, V. ZOBOLI, 
M. Converst, P. FRANZINI, I. MANNELLI, R. SANTANGELO and V. SILVESTRINI: Phys. 
Rev., 108, 1353 (1957); and other papers. 

(2) V. G. SoLoviev: Zu. Éksp. Teor. Fiz., 33, 537, 796 (1957); Nucl. Phys., 6, 
618 (1958); S. N. Gupra: Can. Journ. Phys., 35, 1309 (1957); C. Iso: Progr. Theor. 
Phys. 20, 410 (1958). 
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actions assuming only the CP invariance of the theory. These authors, however, 
have not made a penetrating analysis of the mathematical structure of quantum 
field theory. 

A question seems to be open as to whether the usual theory can still remain 
as a consistent formalism when we have taken away the requirement of its 
invariance under space inversion. Especially, it will be of great interest to 
investigate how the renormalization is carried through for parity non-conserving 
interactions in the framework of the conventional S-matrix theory. 

D’ESPAGNAT and PRENTKI (*) have first investigated the problem of re- 
normalizing parity non-conserving theories. These authors carried out the 
renormalization of a CP invariant theory according to Dyson’s method (+), 
and concluded that the theory is renormalizable in the sense of removing 
divergencies. They also have pointed out that the ways of renormalization 
are not unique. However, the equations determining the renormalization 
constants have not been given explicitly, nor examined the question as to 
whether the equations permit a consistent solution. 

It is the aim of this paper to investigate systematically the above men- 
tioned problems which seem to be of great importance from the point of view 
of principles. 

In the first section we start with an unrenormalized Lagrangian which 
contains a parity non-conserving interaction. Without the use of power series 
expansion, we try to obtain the equations defining renormalization constants 
and, find that the result is given without any ambiguities. In terms of these 
constants, we can formally carry out the renormalization according to KAL- 
LEN (5). In Section 2 the Lagrangian will be rewritten using only renormalized 
quantities. Formulae are given which express the renormalization constants 
in terms of Lehmann’s spectral functions (°). 

Finally, the problem is solved explicitly in the lowest order approximation 
of perturbation theory and it will be shown that the above obtained equations 
permit a consistent solution only when the theory is parity conserving. 


1. — Equations for the renormalization constants. 


Let us consider a system with three kinds of particles which we call nucleons, 
hyperons and mesons. To each kind of particles corresponds a field which 


(?) B. D’Espacnar and J. PRENTKI: Nuovo Cimento, 6, 2989 (1957). 

(4) F. J. Dyson: Phys. Rev., 75, 1736 (1948); A. SALAM: Phys. Rev., 79, 910 (1950); 
J. C. WARD: Phys. Rev., 84, 897 (1951); P. T. MaTtTHEWS and A. SALAM: Phys. Rev., 
94, 185 (1954); and other papers. 

(5) G. KALLEN: Helv. Phys. Acta, 25, 417 (1952); Physica, 19, 850 (1953). 
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will be denoted by y, y and q, respectively. We assume the system to be 
described by the following unrenormalized Lagrangian: 


(1.1) fey = Jen ere 
(1.2) Lo = — ¥(&) (ye + m1) p(x) + my (x)p(x) — 4(@) (yo + mez) y(x) + 
+ dm, 7 (x) x(0) + p*(x) (DO — pw?) p(x) + Ju?p*(x)œ(x) , 


(1.3) 2, = — gl pla) (x + by) x(w) p* (a) + h. c.]. 


From the above Lagrangian we obtain the following equations of motion for 
the field operators in the Heisenberg representation 


(1.4) | (V0+ m,)p(x) = gla + Bys) x(x) p* (x) + Sm p(x) = f(x 


| (yO+ ms) x(x) = g(a*— B*ys)p(x) p(w) + dm, x(a) = g(a) . 
The condition of reproducing the right and left asymmetry: 
(1.5) Re (a*B) 0, 


is to be fulfilled. For reasons of simplicity, however, we shall discuss exclusively 
a special case of CP invariance, that is, both x and f are assumed to be real in 
what follows. Furthermore, to avoid inessential complications, we put m,—m,. 
It is easy to introduce the mass difference into the formalism if required. 

We now solve the equation (1.4) in view of the following asymptotic con- 
ditions: 


PT. 
(1.6) 
| sim, 4 (0) a as = ra x). 


The incoming fields are to satisfy the free field equations 
HE) 


and the covariant commutation relations 


{p(a), par") == — iS (ax = 10) 


{rx}, Fe) = — Sa). 
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On the other hand, as for the operators y and y of interacting fields, it 
can be deduced from the Lagrangian that the anticommutators for equal times 
have the following form: 


ro) A) VAI 
(1.9) 


mo SS 
| 
| 
2 
2 


{x(a), X(0')}, <9 = Ys dla x’). 


It is well known that by virtue of the equations (1.7) and (1.8) it becomes 
possible to define the «particle number». We consider the complete set of 
state vectors |) as eigenstates of the number of incoming particles. Finally 
we assume the existence of the energy-momentum four vector 17 which has 
the property 


(1.10) [P,, F()] = 12, F(a) 


and of the vacuum |0> characterized by the lowest eigenstate of Po. 
Under the above assumptions we can prove the general formula 


(1.11) £0| F(x) |q> = 
= i 9x — x')<0|{F (a), f(a’)}|0> da’ <0 | p(x") |q>+<0| cote |O><0 lp (a) |q>, 
, cpl 


where F(x) is an arbitrary operator and |q> denotes a state with one particle with 
energy-momentum vector q. We shall not describe here the proof explicitly, 
which is quite similar to that in parity conserving theories (7). 

Putting F(x)=w(x) in (1.11) gives 


(1.12) O|y(x)|q> = 


©) i] Oe — æ)<0 | {y(x), f(æ')} 10) da’ <0 | p(x") | > + <0 | pa) |q>, 
where f(x — 2')0|{wp(x), f(x')}|0> can be expanded as 


Lee] 


= | dp exp [ip(x — a')][A(p) + ty,p,B(p) + iy:y,p,B'(p)] . 


(7) Ge KAULEN: Hand ds Pins: 5, 169) (1958). 
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Clearly the last term results from the violation of the invariance under space 
inversion. The absence of the term A’y; is due to the CP invariance. 
Substituting (1.13) into (1.12) leads to 


(1.14) <0 |w(a)|q@> = (N + N'y5)<0 px) |q> , 
where N and N’ are the constants defined by 


N=1+A—mB, 
N= —— mB, 
The fact that A(q), ete., are universal constants independent of x and q, is 


easily verified in a way similar to that in parity conserving cases. 
Similarly, for the hyperon field, 


(1.14a) <0| x(x) |q> = (N — N'y5)<0| x(a) | où . 


It should be remarked that the reality of the constants N and N’ are verified 
by virtue of the CP invariance. 
Next we put F(x)= f(x) in (1.11). Considering the equation 


<O| f(x) |a> = (ye + m)<0 Iy(x)g» , 
= 2mN'y;<0|p(a) |@ , 


we obtain 


(1215) 070, CO | (x) | g> = 
ee INC — 2") 0 |{f(x), f(æ')}|0 da’ <0 | p(x’) |q> + 2mN'y,<0 | p(x) | g>, 


where the expansion analogous to (1.13) is possible: 


(1.16) <0]{f(2), f(a')}|0> = 


= aay [2 exp [ip(w — a')]e(p)[Zx(p®) + (typ + m)22(p*) + tys(yp)*2(p?)] - 


It can be proved that the expansion coefficients are allreal. Substituting (1.16) 
into (1.15) and performing the x integration, we have 


(1.17) Sm, = Dm) + ysm[2N'— F,(— m*)], 
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with 
foe} 
CA 
= DI a 
2,(=m°) = P| tee) a 
JTE 
0 
Similarly 
(1.17a) dm, = E,(—m?) — ysm[2N'— alan?) | 


The above equations implicitly define the mass renormalization. 

To obtain the definitions for the constants N and N’, we utilize the fact 
that the anticommutator for equal times is one and the same for the free field 
and for the interacting field. 

Integrating the fundamental equation (1.4) in view of the initial condi- 
tion (1.6) follows 


and 


+ {y(x) — p(w), pa") + | Sale — 2"){f(a'), {x")} Sa a) dx” dan”. 


We now put %x,=%, and take the vacuum expectation values of both sides. 
Considering the definition of a matrix product and Eq. (1.14), we obtain 


Va O(x — x") = (2N —1)y, d(x — x) + LE (— m2) + 2m ZT) (— m*)] y, 6(x — x’) + 


+ m[2N'— Zi m?)|ysy,6(x — x!) 
and henee 


(1.18) N =1—3[2,(— m?) + 2md\(— m*)], 
(1.19) Nie Em), 
where 
— N na 
Zi(-m?)=—|da asi 
(a — m?)? 
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These equations define N and N’ respectively. Using the result (1.19), Has. 
(1.17) and (1.17a) can be reduced to a simpler expression 


(1.20) Om, = Z;(- m?) r= 1,2 


from which we can see that the mass correction dm is a real c-number not 
involving y;. It is remarked that the equality of the constants for the nucleon 
and for the hyperon is merely a consequence of the assumption m, = my. 

Thus, the equations which determine the renormalization constants are 
given without any ambiguities. The equations for N and dm are exactly the 
Same if the parity non-conserving terms were absent. On the other hand, the 
third renormalization constant introduced on account of the parity non-con- 
servation, is uniquely determined from the requirement that the anticommu- 
tators for equal times of the interacting field operators are to be identical with 
those for the free field, that is to say, they cannot contain any y; terms. This 
requirement is necessary, as KALLEN has emphasized (7), for being able to 
define consistently the «particle number » at arbitrary time. It is of interest 
that this very condition assures of obtaining dm as a e-number. 

Before closing this section, it should be mentioned that, without solving 
explicitly, we cannot be sure of the internal consisteney of the above obtained 
equations. 


2. Renormalized field equations and commutation relations. 


Once the renormalization constants have been determined, the renormali- 
zation can be performed very simply. We construct the renormalized field 
operators 


| p(x) = (N + N'y) (a) , 
(2.1) { 
[| x @) = NY — Wy) x(a) , 


which, as is easily seen, fulfil the renormalization conditions of KALLEN (°°) 


[ <Olp'(w)q|> = <O|p(a)|@ , 
| <9] x(n) |a> = <0 x(a) |a 


Here, it is to be remarked that 


(2.3) N2—N”’ 40 


is necessary for the existence of the inverse operators in (2.1), which can be 
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explicitly given by 


2A N 
(2.4) (1 Et Ys) saat N? A N”2 
Tf we introduce new notations: 

(2.5) Zo = N° NV, 

(2.6) a = NZ", b= N'Z?, 


then the expression (2.1) is brought into the form 


yp (a) = Z,*(a + bys)p(x) , 


| 
(257) 
| x) = 2, (a — bys) x(@) . 


Let us now rewrite the original Lagrangian (1.1-3) in terms of the renor- 
malized field operators. Then, the following expressions are obtained 


(2.8) Ceo ee, 
(2.9) Ly = — Zp(A+ bys)(y8 + m)(a — bys)p" — 


= VAGA = by;)(ye + m)(a + bys) x al 4.0 (LI mu) > 
(2.10) 2, = Z,3m[y'y' + xxv] — 92.23 (a + By) xe +h. c.], 


where 


(2.11) 


We have used, in the above, the well known expression for the renormalized 
meson operator: 


(2.12) p(x) = Zi o(x). 
The definitions for the constants Z, and du will not be described here, since 
they are exactly the same as those in the usual theories. 


Putting 


(2.13) g = 4,'2,2%g, 


D+ 


ON THE RENORMALIZATION OF A PARITY NON-CONSERVING INTERACTION 95 
leads to the coupling constant renormalization 


gx > Zigia', 
(2.14) 


gp > Zig". 


It is to be noted that, considering a and b are real, the transformation (2.11) 
is unimodular but not unitary. 

In what follows we drop the dash on the renormalized operator, since the 
corresponding unrenormalized operator will not be used again. 

The Lagrangian (2.8-10) follows the equations of motion for the renor- 
malized. fields 


(a+ bys)(70 + m) (a — bys)y = Imp — Zi, gla+ By) xp" , 
| (a — by;)(y0 + 


m)(a + by;)x = my — Z,g(~— Bys)y@ , 
and the anticommutators for equal times 


(2.16) I (9(2), PO Ne, a = Za (a + by)yp(a—Dy;) 0(x — x’) , 
2.16 | | 

[ {X(2), Moi = Z, (a = bys)yala op bys) Ô(x — x’). 
If we now put 

(2100) ra == (a + bys)y,(4 + bys) , 


it can be verified that the va Satisfy Dirac’s relations as well as Va 


(2.18) Dune 4 DAT = 26.,, 
and 
(2.19) o = Vini = 5. 


With the use of these expressions, (2.15) and (2.16) can be rewritten as 


À (y 0 + my = Smy — Zig(a + Bys)x@™ , 

(2.154) ; 

| (y0+ m)y = my — Z;g(x — Bys)yy , 
I {p(x), G(@')}, a = 2, dea), 

2.16 jar, | 

fe | {y(æ), %(a let, ie OX — x). 
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From the above, we can conclude as follows: 

In the renormalization of a system where parity non-conserving inter- 
actions are present, i) there appears a certain difference in the y-represen- 
tation between the physical and the bare particles, apart from the mass and 
the coupling constant renormalizations carried out in the conventional way. 
ii) As far as the y-representation is taken to be equal for the two kinds of 
bare particles, then occurs a definite relative difference in the y-representation 
between the physical nucleons and the physical hyperons (*). iii) These dif- 
ferences in the representation is non-unitary because of the reality of a and b, 
and hence it is impossible to give y_,, va and yi unitary representations at 
the same time. | 

The renormalized propagator is now given by 


(2.20) <0| Ty(a)p(x')|0) =—48 


According to LEHMANN (°), we introduce the spectral representation 


(2.21)  8,(p) =— Zi [dee 


0 


iy po, (% (522) 0» (x? eS wey Dos(x?) 


p? + x? — ie ; 


th 


from which we get the expression 


(2.22) O|{y(x), p(a’)}|0>, ao = | fo 2) dx? + ve] 03(%?) de vi o(a— x’), 


for the anticommutator for equal times. Comparing this with (2.164) leads to 


(2.23) Zi =[evayaes, 
‘o 
and 
(2.24) = o| 05(242) da 
0 
An inequality 
(2.25) Oa Fo | 


(*) On the contrary, to obtain the equal representation for the two kinds of physical 
particles, it is necessary to take certain shifted representations in the corresponding 
unrenormalized quantities. This means, however, that the free part of the unrenor- 
malized Lagrangian also violates parity. 
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can also be proved in a way similar to that in parity conserving theories. 

Here, N? > N” has been assumed. The alternative case N?Z N'° must be 

excluded because the negative sign appears on the right hand side of the anti- 

commutator. This is a situation similar to the ghost in Lee’s model (8). 
The mass renormalization is given by 


(co) 


(2.26) om = Za move) + 02(x?)| dx? . 


0 


And it is to be noted that here use has been made of the condition (2.24). 


3. — Perturbation theory in the lowest order approximation. 


So far we have discussed the problem without the use of expansions in 
the coupling constant. In this section, we shall solve it explicitly in the lowest 
order approximation of perturbation theory to see in more detail the pro- 
perties of the renormalized field. 

Now, in the equation 


(3: L) (48 + m)(y 70 — m)<0 y(x)y(x')|0> = 


= g? (a ae Bys) 0 | (ar) PAC) 0 (oo — Bys) 0 | p* (ar) p(x’) LOS i 


we expand. 
(3.2) <0 | w(x)p(a’)|0> = 


=| ‘KO ko) | AeA iy oye!) + (x?) + psy Pkos(x?)] Ok? + x? ] , 


È 
(Gee) <0] x(a) ya’) | 0> =| d4qO(qo) (yg — m) ol g? + m?] . 


Hence, the Fourier transform of (3.1) is 


(3.4) 0(ko) (OR + m) (iy? ho, + 02 + Wysy kos) (Pk + m) = 
n. 


= (og | UG Allo — do) do) AL + MIA (E — ge + 1] 


(ode Bs) (gm) Bys) : 


(8) T. D. Lee: Phys. Rev., 95, 1329 (1954); G. KALLÉN and W. PauLi: Dan. Mat. 
Fys. Medd., 30 no. 7 (1955). 
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We multiply this first by unity, by ik from the right, and then by y; from 
. the left. Taking trace in each step gives respectively 


2mr o, + (m?+ x?)e, = — max? — B*)I 
(3.5) (m2 + x?) Ao,+ 2mAo.+ (m°— x?) Bos = (a? + b?)J 
(m?-- x?)Bo,+ 2mBo,+ (m?— x?)Ao; = 2afpd , 
where 
na 
Ge) T= — pl [arg hy — IMMA MATE — a)? + = 
A Cent en A 
| 16772 x? , 
Ped ber anor — qo) 9 (qo) (kg) OL g? + m?] oL (& — ga)? + u?] = 
à 1 


“(4a Ma) Ee 


[| A = (a? + b*)?+ (246), 
(357) 4 
| B = 2(2ab)(a? + b?). 


Solving the simultaneous equations (3.5), we have 


(3.8) re a JCA)T2a BA — (a? + f?)B 


m2 — 


In order to satisfy the condition (2.24), it must be 


(3.9) 2aBA = (a? + B?)B 
1.€. 
(3.10) o=a+b?, f= ab. 


By using this, the other spectral functions 0,(%?) and o,() are completely 
determined, from which the renormalization constants Z, and 3m can be cal- 
culated with the aid of the formulae (2.23) and (2.26). 
for them will not be given here. 

The point to be remarked is rather the very relation (3.10). This shows 
that, not only the integral of the form (2.24), but the function 03(x?) itself 
is identically zero at least as far as the lowest order approximation is concerned. 


The explicit expression 


co 
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In other words, there is no room in the spectral decomposition of the renor- 
malized propagator to introduce consistently the third o function charac- 
teristic of the violation of parity. 

On the other hand, the renormalized Lagrangian 


(3.11) 2 = — Lp(a + bys)(y0 + m)(a— by,)yp 


— Z,4(a — by:)(y0 + m)(a + bys)x 
— Z, g(x + Bys)xp* — h. c. 


is rewritten, in view of (3.10), 


(Le) L = — Z,p(ye + m)y 


— 2:%(y0 + m)x 


— Z19p 19" — h. c. 


in terms of the conventional canonical variables 


(3.13) 
which satisfy 


(ie 
(3.14) | 


As can be seen, the above obtained Lagrangian does not contain any terms 
which violate parity. We have found a parity conserving Lagrangian by re- 
normalizing the parity non-conserving one. This means, however, that the 
original Lagrangian also cannot contain parity violating terms in order that 
the theory be consistent. The theory is thus reduced to the usual parity con- 
serving one. 

Here may be raised an argument that we shall be able to have a renor- 
malized parity non-conserving Lagrangian, if we take, at the outset, certain 
shifted representations between the two kinds of unrenormalized particles 
(see the footnote of the previous section). To examine this possibility, sup- 
pose that, instead of (3.11), a renormalized Lagrangian of the form 


(3.15) L = — Z,:y(y0 + m)y 


— Z,4 (yo + m)x 


— Zagp(x+ Bys)x — h. c. 
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be given. This Lagrangian follows the anticommutator for equal times 
a \= _ 7-1 Dr 

(3.15) {p(a), Pe')}, yi = Za Vada — x): 


On the other hand, the renormalized propagator 


à q' De farre PQ?) + a(x?) + dy5yDoax?) 
QE Sale) == zi fan pra ON 


0 


(3217) fut), pa)h=s = | face dx? + va | eal) de yaÔ(x — x’). 


~ 


Comparing this with (3.15) gives 


© 


(3.18) Jostera 2 


0 


Now, let us calculate 0, in the lowest order approximation of perturbation 
theory. Then we find 


(3.19) CE) = JT) 


which is inconsistent with the condition (3.18), unless «6 = 0. 


4. — Concluding remarks. 


Opinions will be divided on the interpretation of the revealed phase of the 
theory. First of all, it might be possible to consider that the inconsisteney 
mentioned is only a result of the unjustified use of perturbational calculations, 
and there will exist an exact solution which cannot be expanded in powers 
of the coupling constant. 

If, contrary to this expectation, it turns out that the inconsistency is inherent 
to the theory, then we shall conclude that the present formalism of quantum 
field theory is to be drastically changed in order to reconcile the notion of 
renormalization with the fact of parity non-conservation. 

The third point of view is also possible which asserts that the notion of 
renormalization cannot be applied to the parity non-conserving interaction. 
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In view of the fact that parity violates only in weak interactions associated 
with the decay of elementary particles, it may be reasonable to presume that 
there exists a complementarity between the concept of the dressed particle 
and the fact that the particle decays spontaneously. 


The author wishes to thank Mr. T. YosHimura for stimulating discussions 
and criticisms. 


RIASSUNTO (*) 


Si da una formulazione esatta del metodo di rinormalizzazione in una teoria non 
conservante la parità. Si ottengono equazioni che determinano le costanti di rinorma- 
lizzazione, e, fra l’altro, si esamina la compatibilità interna della teoria. 


(*) Traduzione a cura della Redazione, 
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NOTE TECNICHE 


Detection of High Energy u-Mesons 
by an Air Cerenkov Counter (*). 


R. GIACCONI (+), W. BLUM (+) and G. T. REYNOLDS 


Palmer Physical Laboratory, Princeton University - Princeton, N. J. 


(ricevuto il 16 Settembre 1958) 


Summary. — An investigation of the properties of air Cerenkov counters 
is described. A quantitative interpretation of the results is given. 


1. — Introduction. 


The interest in high energy particles selectors has led various groups (1?) 
in the past to study the properties of Cerenkov counters filled with gases. 
In order to obtain a better understanding of the possibilities of these counters 
as reliable ultrarelativistic particle detectors, it was felt necessary to repeat, 
with the improved photomultipliers now available, some of the experiments 
already performed and to present a more detailed interpretation of the data 
obtained. The purpose of this paper is to describe such an experiment and 
to present a quantitative analysis of the results. 


2. — Experimental arrangement and procedure. 


2°1. The experiment was performed at sea level with u-mesons (energy 
> 280 MeV) of the cosmic radiation. The apparatus consisted of a light tight 
box with a parabolic mirror at the bottom. The Cerenkov light produced 
by single particles in the air contained in the box was focused on the cathode 
of a photomultiplier. The geometry of the box and mirror and the dispo- 


(*) Sponsored by the joint program of the U.S. Office of Naval Research and 
the U.S. Atomic Energy Commission. 

(+) Fulbright fellows. 

(1) A. Ascott BALZANELLI and R. AscoLi: Nuovo Cimento, 11, 562 (1954). 

(2) F. R. BarcLay and J. V. JELLEY: Nuovo Cimento, 11, 26 (1955). 
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sition of the G.M. counters, absorber and phototube are shown in Fig. 1. 
A Dumont 6263 photomultiplier with 12.5 cm cathode diameter was used. 
The amplitude of the solid angle and the height at which the phototube was 
placed were selected according to the results of an empirical focusing pro- 
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Fig. 1. — Essential features of the apparatus. 


cedure in order to assure complete light collection for all possible paths of 
the incoming particles. The 18x20 cm? G.M. tray (G) on top of the photo- 
tube prevents the acceptance of pulses produced by particles passing through 
the photomultiplier itself, 


22. The experiment consisted in recording simultaneously the rates of 
coincidences of the type 0,C,PGA,A,A;=N, and 0,0,GA;A4,A43=N,. The 
resolving time of both coincidences was approximately 2-10-*s. Time con- 
stants of 2-10-7s were used in the circuit of the phototube amplifier. The 
signal from this amplifier was fed through a single channel differential pulse 
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height analyzer to the coincidence circuit. For each setting of the discrimi- 
nator the N, rate was measured, once with the photocathode surface unco- 
vered, (N,),, and the second time with the photocathode surface covered by 
a shutter, (N,),, to obtain the rate of background counts. This latter rate 
was of the order of 30% of the total counting rate and was in agreement 
with the expected chance coincidence value due to the noise of the photo- 
tube and amplifier and to the rather poor time resolution. All measurements 
were started from some minimum value of the discriminator setting below 
which the statistical fluctuations of the number of background counts became 
large compared to the difference between the total and background counts. 
Since, even at larger values of the discriminator bias, the background was a 
considerable fraction of the total number of counts, care has been taken to 
avoid the effects of slow changes of the noise level by taking the counts with 
shutter open and closed, for the same pulse height interval, one immedia- 
tely after the other. In order to measure the possible influence of fast chan- 
ges of the noise, various counts for a given discriminator setting were repeated. 
The difference between total number of counts and background were found 
to be constant for the same pulse height interval within the statistical 
errors. 


23. Between measurements, periodic checks of the overall amplification 
of the photomultiplier circuit were performed using a Cs and a plastic 
scintillator with constant geometry. The peak given by the 624 keV electrons 
was determined with 20%, resolution. 


24. The pulse height of a single photoelectron pulse was calculated from 
the experimentally measured noise spectrum of the phototube (*). This value 
depends on the overall amplification of the photomultiplier and associated 
circuit. Since the experiment was repeated with two phototubes and it was 
felt desirable to be able to compare directly the results, the voltage applied 
to the tubes was adjusted so that both would give the same pulse height for 
the electron peak of the !Cs source. Under these conditions the amplitudes 
of the single photoelectron pulse height, V,, were found to be V,, = 3.5 and 
Vos =1.5 V. It can be shown that Vo,/Vor = Epn2/fn1, Where em represents the 
quantum efficiency of the phototube. 


3. — Results. 


Because of the coincidence requirements and the thickness of the lead 
absorber over tray 0,, we feel confident that our N, counting rate represents 
the flux of single particles and more specifically of u-mesons with energy 
greater than 280 MeV through the detector. The corrected counting rate 
of 1075-436 per hour is in good agreement with the value that can be cal- 
culated from Rossi’s data (*). The total running time with each phototube 
was approximately 200 h. For each pulse height interval the measurements 
with shutter open and closed were extended for a period of time sufficient 
to collect a statistically significant number of events. The ratios (NN) 
with shutter open and (N,/N,). with shutter closed were calculated for each 


(#) G. A. Morton: ROA Rev., 10, 525 (1949). 
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discriminator setting. The 


difference, (N,/N.,).— (NJN,), represents 


fraction of u’s emitting a quantity of light corresponding to the given pulse height. 


In Figs. 2 and 3 the 
differential spectrum of 
the difference, (N,/N,,), — 
(NN), versus pulse 
height is given for the two 
runs with different photo- 
tubes. The errors to be 
associated with the exper- 
imental values are of the 
order of 10%, at all points. 
The electron scales atta- 
ched are derived from 
the measured one-electron 
pulse height. 


Fig. 2. — Differential pulse 
height distribution. (Photo- 
multiplier 1). 
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The sum of (N,/N,).— (N,/N,) on all pulse height intervals represents 
the measured value of the fraction of p’s producing Cerenkov radiation. This 
value was found to be S,= (19.9+2.0)% and S, = (25.5+3.5)% for the two 
phototubes respectively. According to the known y spectrum (1), 27% of the us 
which pass through the counter (H > 280 MeV) should produce Cerenkoy light. 


4. — Discussion. 


A discussion of the general principles of the experiment is given in ref. (?). 
We are interested here in the derivation of the expected pulse height distribution. 

From the classical treatment of Frank and Tamm, the average number 
of photons emitted in the wave length interval for which the photomultiplier 
is sensitive (approximately (3 700 --5 200) À) can be calculated as à function 
of the energy of the incident particles. Assuming, further, a given energy 
spectrum for the u-mesons, namely W(H)dH = E-*°d£E, the distribution of 
pulses of a given number of photons can be found. If one also knows the 
efficiency of collection of light, &.n, and the quantum efficiency of the photo- 
tube, en, the pulse height distribution can be directly calculated. However, 
since the average amount of Cerenkov light created by a single particle in 
air is generally small, even for distances of the order of meters, one has to 
take into account the statistical character of photon creation, photoelectric 
emission, and secondary electron emission at all dynodes. 

Since the number of photoelectrons at the cathode is an order of mag- 
nitude smaller than the number of photons from the Cerenkov effect, and 
smaller than the number of secondary electrons at the dynodes, we have used, 
as a first approximation, only the statistical fluctuations in the emission of 
photoelectrons. For these we have assumed a Poisson distribution (3). It can 
be shown that the resulting differential distribution expressed as a function 
of the number of photoelectrons, À, is given by: 


7 1 K à © °K n il 
Wie eo 5 ee er 
A tt! 5 n! 2(n—1) —1 
where: ¢ is the product ewn'em (both quantities previously 


defined); 
K=4H(1-1/N?) is a constant which represents the maximum average 
number of photons emitted (approximately 80); 


H is a constant of the order of 105; 

N is the index of refraction of air at 1 atm, 20 °C 
CN == 1.00027): 

A is a normalization constant; 


and where the approximation H — Ka H has been used. 

To compare the calculated distribution with the experimental results, the 
fit with various values of eA was tried using a least square method. The cal- 
culated histograms associated with the experimental results, Fig. 2 and Fig. 3, 


are the ones for which the best fit was obtained. The corresponding values 
of eK are (ef), = 3.0 and (ek), = 11.0 for the first and second photomultiplier. 
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Assuming en to remain constant for the two different runs, it follows 
that (eK),/(e)) = Esni/Emns = .27. This value should be compared with the 
experimentally determined e,m1/€n0 = Vos] Vos = -43 which has an estimated error 
of 40%. Using the calculated histograms, the experimental values of the fraction 
of u's producing Cerenkov light, 8; = (19.9+-2.0)% and S, = (25.5+3.5)%, can 
be corrected to take into account the finite value of the minimum bias level. 
The corrections yield $,=(21.0+2.1)% and S,= (25.6+3.5)% to be com- 
pared with the expected 27% 


5. — Conclusions. 


The agreement between calculated and experimental values seems to justify 
the approximations used and to indicate no appreciable discrepancy from 
what would be expected under the assumption that most of our light pulses 
were produced by Cerenkov radiation of single u-mesons of energy > 4 GeV. 
We feel confident that with minor improvements in the electronics to reduce 
the background noise, our apparatus could be used as an efficient detector 
of particles with y=1/(1— 62)? greater than 40. 


RIASSUNTO (*) 


Si presenta un esame delle proprietà dei contatori Cerenkov in aria. Si da un’inter- 
pretazione quantitativa dei risultati. 


(*) Traduzione a cura della Redazione. 
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TL NUOVO CIMENTO NOT Bells IN L 1° Gennaio 1959 


Grain Size and Grain Density in Ilford K-5 and L-4 Emulsions. 


R. C. Kumar (*) 


Pure Physics Division, National Research Council - Ottawa 


(ricevuto il 20 Settembre 1958) 


Summary. — Electron micrographs of grains in unprocessed and processed 
Ilford K-5 and L-4 emulsions are shown. The silver bromide grain sizes 
in the unprocessed K-5 and L-4 plates have been measured and found 
to be normally distributed with mean values of (0.210+.007) um and 
(0.134.004) um and standard deviations of 0.035 ym and 0.021 um 
respectively. The electron micrographs of the silver grains in processed 
emulsions reveal absence of spherical shape or solid structure in case 
of most of these grains. The grain densities on plateau ionization tracks 
have been found to be approximately the same in K-5 and L-4 plates 
as in G-5 plates processed under similar conditions. The L-4 plates 
however are subject to more serious fading with delay in processing 
after exposure than the K-5 or G-5 ones. 


1. — Grain size in unprocessed emulsion. 


The size of the silver bromide grains in nuclear emulsions and its distribution 
is of importance in the estimation of loss of energy by ionization of charged 
particles in these emulsions. Electron microscope photographs of these grains 
have been published for Ilford G-5 emulsions (42). Recently Ilford Ltd. have 
introduced two new types of emulsions — K-5 and L-4 — with smaller grain 
sizes, and it was thought useful to investigate grains in unprocessed and 
processed emulsion by electron microscopy. 

The specimens for the electron microscope were prepared by dissolving 
a piece of dry unprocessed emulsion with warm de-ionized water, the solutions 
being well stirred. The emulsions used were from Ilford batch No. Z711 
(April, 1958). Specimens were collected on standard holders (copper grids 


(*) National Research Council Post-Doctorate Fellow; now at University College, 
London. 

(*) G. BARONI and C. CASTAGNOLI: Nuovo Cimento, 7, 394 (1950). 

(?) E. Pickup: Can. Journ. Phys., 31, 898 (1953). 


108 


R. C. KUMAR 


a) b) 


c) 


Fig. 1. — Electron microscope photographs: a) L-4; b) K-5; e) latex grains of known 
uniform diameter (0.340-+.005) um; all at same magnification. 


R. C. KUMAR 


b) 


Bio. 2. — Silver grains in processed emulsions: a) L-4; b) 


K-5: magnification approxi- 
mately same as for Fie. 1. 
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covered with Formvar films) for use in a 50kV RCA electron microscope. 
Figs. 1a and 1b show reproductions of the silver bromide grains in L-4 and K-5 
emulsions respectively. Fig. 1e shows some latex grains of known uniform 
size (0.340 +.005) um under the same magnification. Measurements on 300 grains 
gave mean grain diameters of (0.134+.004) um and (0.210.007) um for L-4 
and K-5 emulsions respectively, where the errors quoted are standard errors. 
The distributions of the grain sizes shown in Figs. 3a and 3b were found to be 
normal with standard deviations of 0.021 um and 0.035 um for L-4 and K-5 
emulsions respectively. 


100 v 
* —— Observed 80 = 
800 2 --- Normal = — Observed 
a curve fit BO ra --- Normal 
60 E curve fit 
5 40-5 
40F & 2 
$ 20 S 2 
20F § = 


Ee shah PEN PP) 0 3909082 705 
Grain size in um Grain size in ym 
Fig. 3. — Grain size distributions: a) L-4; b) K-5. 


It will be noticed that in Figs. 1a and b there are groups of grains in addition 
to what appear to be single resolved grains. It is not certain whether these 
groups represent grains as they were originally packed in the emulsion or 
whether they were formed while preparing the specimens. Measurements were 
made only on resolved grains. Besides, some of the single grains showed some 
structure at the edges. This is due to the action of the electron beam on 
the silver bromide grains causing migration of silver ions in the grains; similar 
phenomena have been observed for Ilford C-2 and G-5 emulsions (**). In making 
measurements of the diameter, such grains showing any structure at the edges 
were ignored. Care was also taken, while making the electron microscope 
exposures, not to use too intense an electron beam, or a long time of exposure 
such that visual observation indicated structural changes in the grains. 


2. Grain density. 


Two stacks of stripped emulsions, each containing G-5, K-5 and L-4 pellicles, 
were exposed to electrons and positrons of energy (300 —700) MeV at the Cornell 
synchrotron. These stacks were processed at different times, the first two 
days and the second fourteen days after the exposure. Prior to development 
and after exposure, the emulsions were transported and stored at low tem- 
perature, about 32°F. The usual temperature development method was em- 
ployed, the temperature of the hot bath being 78°F and the time of develop- 
ment 40 min. 


(3) L. WINAND: Fundamental Mechanisms of Photographic Sensitivity (London, 
1951), p. 286. 
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The grain density on the electron (positron) tracks of plateau ionization 
was found to be approximately the same in G-5, K-5 and L-4 emulsions of the 
first stack, being about 26 blobs/100 um. The second stack had undergone 
obvious fading; the grain density on electron (positron) tracks of plateau 
ionization was about 18 blobs/100 um in the G-5 and K-5 plates and the 
L-4 plates showed no tracks at all. 


3. — Grains in processed emulsion. 


Figs. 2a and 2b show electron micrographs of silver grains in processed. 
L-4 and K-5 emulsions at approximately the same magnification as in Fig. 1. 
In contrast with the undeveloped silver bromide crystals, most of the silver 
erains do not possess either a regular spherical shape or a solid structure, 
This agrees with previous findings (4) that the developed silver grain consists 
of a mass of tangled filaments. However, it is not certain that the structure 
of the grains was not somewhat modified in the process of centrifuging that 
was found necessary in preparing specimens for the electron microscope due 
to the relatively sparse concentration of grains in the processed emulsion. 
The mean grain size was not measured for the processed grains because of 
the difficulty of assigning a parameter common to all (such as the diameter 
in case of spherical grains). Visually the size of the developed silver grains 
appeared to be larger than that of the silver bromide grains in the unprocessed 
emulsion as a comparison of Figs. 1 and 2 will show. The distribution in size 
of the silver grains was less regular than that of the silver bromide grains. 
This may be due to the presence of the smaller background grains, in addition 
to those forming tracks in the emulsion. 

AR 


The author is greatly indebted to Dr. E. Pickup for valuable advice and 
to Dr. S. T. BAyLEY of the Applied Biology Division for kindly making the 
electron microscope exposures. Thanks are due to Drs. B. JUDEK and 
J. G. MCEWEN for help in processing the plates. 


(*) C. E. Hatt and A. L. SCHOEN: Journ. Opt. Soc. Amer., 31, 281 (1941). 


RIASSUNTO (*) 


Si presentano micrografie elettroniche di granuli di emulsioni Ilford K-5 e L-4 
sviluppate e non sviluppate. Si sono misurate le dimensioni dei granuli di bromuro 
d’argento delle lastre K-5 e L-4 non sviluppate e si sono trovate normalmente distri- 
buite con valori medi di (0.210+.007) um e (0.134+.004) um e scostamenti di 0.035 um 
e 0.021 um rispettivamente. Le micrografie elettroniche dei granuli d’argento delle 
lastre sviluppate rivelano per la maggioranza dei granuli assenza di forma sferica o 
di struttura solida. Si è trovato che la densità dei granuli delle tracce mostranti ioniz- 
zazione massima è nelle lastre K-5 e L-4 approssimativamente uguale a quella delle 
lastre G-5 sviluppate in condizioni simili. Le lastre L-4, tuttavia, sono soggette, tardando 
a svilupparle dopo l'esposizione, a maggior sbiadimento che le K-5 o le G-5. 


(*) Traduzione a cura della Redazione. 
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(La responsabilita scientifica degli scritti inseriti in questa rubrica è completamente lasciata 
dalla Direzione del periodico ai singoli autori) 


Presence of °Co in the Atmosphere (*). 


L. Marquez, N. L. Costa and I. G. ALMEIDA 


Centro Brasileiro de Pesquisas Fisicas - Rio de Janeiro 


(ricevuto il 29 Novembre 1957) 


We have established the presence of 
the radioisotope 5°Co in the atmosphere. 
It has been done by collecting a few 
thousand liter of rain water, adding Co 
carrier to it, and passing it through an 
ion exchange column to retain the ca- 
tions. The technique is similar to the 
one we used to determine ??Na in the 
atmosphere and it has been described 
already (1). i 

The cations with the Co were eluted 
from the column and the Co precipitated 
with HS in basic solution. The Co was 
dissolved, then scavenging was done with 
sulphides of Cu, Sb, Te, Se, Bi, and Pb; 
with sulphate of Ba; and with hydroxides 
of Fe and La. The Co extracted from 
thiocyanate solution with a mixture of 
amyl alcohol and ethyl ether, returned 
to H,0 phase, precipitated in solution 
of HCl with 1-Nitroso-2-Naphtol, and 
ignited to the oxide. The procedure 
described is specific for Co, and the chem- 
ical yield ranged from 15 to 60%. 

The Co was counted in a well crystal 
scintillation spectrometer Model 516, 
built by Baird-Atomic, Inc., Cambridge, 


(*) This work was done under the auspices 
of the Conselho Nacional de Pesquisas of Brazil. 

() L. MARQUEZ, N. L. Costa and I. G. AL- 
MEIDA: Nuovo Cimento (to be published). 
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Mass., which we calibrated for energy 
measurements and for efficiency. The 
spectrum obtained showed one peak with 
an energy of (12545) keV corresponding 
to the tabulated y-rays of Co of 123 keV 
and 173 keV. We observed, of course, 
one peak only due to the poor resolution 
of the spectrometer. The counting effi- 
ciency at the photoelectric peak was 
a9 foe 

We carried out three experiments and 
the results are shown in Table I. The 
average activity of 77Co for these three 
experiments is 0.020 dpm/liter. 


TABLE I. 
Amount Absolute | 
Date of rain specific | 

| 1957 water activity 
| in liters | in dpm/liter 
| | | 
| ci — Sp Se ———— ——— n 
| June 8 500 0.016 | 
| September 3 500 0.011 | 
| October 3 000 0.040 


One does not expect that 57Co will 
be formed with appreciable yield in 
nuclear tests, whether they derive their 
energy from fusion or from fission. To 
see that, we will write the reactions that 
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can produce 57Co in decreasing order of 
likelihood. They are ’Fe(p,n)5Co; 
56Fe(p, y)®’Co; *Fe(d,n)*’Co; Fe(p, 2n) 
Co; 8Ni(n, An) Ni; 5°Co(n, 3n)°’Co; 
etc. 

Examining these reactions, one ar- 
rives at the conclusion that nuclear 
tests having mostly fission could not 
account for the 57Co observed, since all 
the reactions indicated produced by neut- 
rons require threshold of 12 MeV or 
greater, and one can calculate that this 
is not compatible with the amount of 
the fission product 137Cs in rain water 
which is 1.2 dpm/liter (?). 

The two first reactions written are 
the most likely to produce **Co in a test 
having mostly fusion. They should be 
hindered by the potential barrier. The 
only thermonuclear reaction giving pro- 
tons with energy greater than the barrier 
in Fe is the reaction *He(d, «)p; but in 
a device where this reaction takes place, 
the reaction d(d,*He)n will be more 
prominent and give a large number of 
neutrons. In any case, in a fusion nuclear 
test where protons with a few MeV are 
produced, it is likely that more neutrons 
will be produced, and since there might 
be iron in the device, we decided to look 
for the radioisotopes formed in iron by 
neutron capture, namely Fe and 5Fe. 
One expects that (n,y) will be more 


(2?) L. MARQUEZ, N. L. Costa and I. G. AL- 
MEIDA: An. Acad. Bras. Ciencias (to be published. 
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probable than (p,n) and (p, y) in iron 
in a nuclear test. 

A sample of 3000 liters of rain water 
was processed, and the iron purified with 
a chemical yield of 50%. The sample 
was counted with a Geiger counter and 
no detectable activity was found, giving 
an upper limit of 1epm. We corrected 
this upper limit for absorption, self- 
absorption, back-scattering, geometry, 
efficiency and yield and it gives an upper 
limit for **Fe of 0.010 dpm/liter. The 
corresponding upper limit for 5°Fe is 
0.003 dpm/liter. Both upper limits are 
lower than the activity of °’Co. 

In conclusion, the presence of °’Co 
in the atmosphere has been proven. The 
present evidence from theoretical con- 
siderations and from the absence of 
radioisotopes of iron seems to indicate 
that it is not likely to come from a 
nuclear test, although this point is not 
definite yet. The systematic observa- 
tion of “Co in rain water and the radio- 
isotopes in its neighborhood during three 
or four years will help to decide its 
origin. 


Note. — We have made more mea- 
surements from January to July, 1958 
on 5Co and 5*Mn. These measurements 
indicate that the bulk of 5Co is pro- 
duced in thermonuclear tests (3). 


@) Le MARQUEZ: No in. Cosma and TC 
ALMEIDA: Proceedings of the Second International 
Conference on the Peaceful Uses of Atomic Energy. 
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A New Type of Particle Detector: the « Discharge Chamber». 


S. Fuxur and S. Mryamoro 


Department of Physics, Osaka University - Osaka 


(ricevuto il 26 Settembre 1958) 


Recently, the devices based on gas- 
eous discharges which can determine the 
trajectory of ionizing particles, such as 
the hodoscope chamber (1) or the trig- 
gered spark counter (2) have been re- 
ported. However, these have some in- 
herent imperfections, in spite of their 
wide usefulness as detectors. Regarding 
the former, even taking into account 
the considerations of F. ASHTON et 
al. (3), one cannot expect a very pre- 
cise determination of the trajectory of 
the particle, owing to the finite dimen- 
sions of the glass tubes and their dead 
spaces, and in the latter it is difficult 


@) M. Converst and A. GOZZINI: Nuovo 
Cimento, 2, 189 (1955); M. CONVERSI, S. FOCARDI, 
C. FRANZINETTI, A. Gozzini and P. MURTAS: 
Suppl. Nuovo Cimento, 4, 234 (1956); G. BAR- 
SANTI, M. CONVERSI, P. MURTAS, C. RUBBIA and 
G. TORELLI: Report at the Geneva International 
Conjerence (1956); M. GARDENER, S. KISDNASAMI, 
E. Rosste and A. W. WOLFENDALE: Proc. 
Phys. Soc., B70, 687 (1957); S. FUKUI and S. 
Mryamoro: INS-TCA, 10 (1957) and 11 (1958), 
(Technical reports at the Institute for Nuclear 
Study, University of Tokyo). 

(2) T. E. CransHaw and J. F. DE BEER: 
Nuovo Cimento, 5, 1107 (1957). 

() F. ASTON, S. KISDNASAMY and A. NV. 


WOLFENDALE: Nuovo Cimento, 8, 615 (1958). 
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to determine all trajectories when more 
than two particles pass through simul- 
taneously, because the discharge, which 
grows in advance of others, decreases 
the intensity of the electric field very 
rapidly. We describe here the possibility 
of constructing a new type of detector 
also based on the gaseous discharge. 
This detector secures the precision of 
the trajectory of the particle and works 
for a group of particles. 

The chamber consists of glass boxes 
filled with a mixture of neon and 
argon gases at atmospheric pressure. 
Each glass box is constructed of glass 
plates glued with Araldite adhesives. 
The inner size of the box used for the 
test is 8.5 cm x 13 cm X 2 cm or xlem 
and the thickness of the glass plates is 
0.5 mm or 2 mm. Metal plates or con- 
ductive glass plates (Nesa) are used as 
electrodes. We also constructed the 
boxes with Nesa whose outer surfaces 
were conductive. 

A pulsed, intense electric field is 
applied across the spaces immediately 
after the passage of the ionizing particle. 
If the intensity and the duration of the 
electric field are properly chosen, tiny 
luminous discharges occur along the 
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trajectory of the particle and can be 
photographed. When the electric field 
is kept longer than a certain time at an 
intensity sufficient to cause discharges, 
successive discharges oecur at random 
places. It seems that these spurious 
discharges are caused by the photo- 
electrons freed by photons emitted from 


excited atoms inside the streamers ini- 


tially produced. Therefore, the time 
width of the high voltage pulse must be 
very short in order to prevent spurious 
discharges. 

The arrangements of the chamber for 
the test are shown schematically in 
Fig. 1 and the circuits of the triggered 
high voltage pulse generator is shown 
in Fig. 2. As seen in Fig 2, the pulse 
has an exponential shape with time 


constant CR. We can obtain good 
HV 
Trigger 
pulse 


Chamber 


Fig. 1. — The arrangements of the chamber for 
test are shown schematically. In Arr. A, tracks 
are photographed through the conductive glass 
and in Arr. B, through the side of the box. 


tracks when we employ the value of 
10-7s for the time constant, i.e. 
6 =0.001 uF and R=100Q. From 
the results of the present investigation 
the value of OUR must be shorter than 
10-68. Some typical examples of photo- 
eraphs are presented in Fig. 3. Fig. 3a 
shows the two photographs taken with 
Arr. A and Fig. 3b is the one taken 
with Arr. B of Fig. 1, respectively. 


FUKUI and S. 
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Specifications of these chambers for 
the present test are as follows. The 
time constant of the voltage pulse is 


“ou 


Arrangement ‘A" 
Side view 
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Fig. 2. — The circuit of triggered high voltage 


pulse generator. 


fixed at the value of 10-7s throughout 
the investigation. The efficiency for the 
formation of the track is not changed 
by the value of the applied voltage 
from 12.5 kV/2cm to 15 kV/2 cm in the 
case of the box of 2 em width and from 
8 kV/1 em to 12 kV/1 cm in the case of 
the box of lem width. The efficiency 
of the chamber placed in Arr. A is 
almost 100% except at the spaces near 
the glass plates when the time delay 
between the passage of the particle and 
the application of the voltage pulse is 
shorter than 2 us. When the time delay 
is about 10 us, discharges showing the 
track do not occur and randomly sepa- 
rated discharges occur at a few places 
in the chamber. Therefore, the maximum 
sensitive time of the chamber is expected 
to be 10 us. The inefficient space close 
to the glass plate, which is somewhat 
larger near the positive electrode, is 
caused by the discharge mechanism and 
depends slightly on the pulse height. 
Its width is about 3mm along the 
surface of the plate of the positive 
electrode in the case of Arr. A. The 
width of the track is ~2mm and 
the average number of discharge columns 
(this quantity is like the drop density 
of the cloud chamber track) is 2.6/em. 
The width of the track in the case of 
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Fig. 3. — Examples of photographs taken with the discharge chamber. (a) shows the two taken 

in Arr. A: one is a track of a penetrating particle and the other shows tracks of an electron shower. 

(b) is a photograph taken in Arr. PB, a track of a penetrating particle. These were taken under the 
condition of CR= 1077 s, Ereak = 13 KV/2.cm tgelay = 1 US. 


A NEW TYPE OF PARTICLE DETECTOR: THE « DISCHARGE CHAMBER ) 


Arr. B is ~3mm. The recovery time 
was determined experimentally to be 
0.1 s. As one sees in Fig. 3a, this chamber 
can detect many particles passing the 
same box simultaneously and the spatial 
resolution is expected to be ~ 2mm. 

The details of this investigation will 
be reported elsewhere in the near future. 
As described above, this detector is 
simple in construction and stable in 


19 
rd 


115 


operation, and also has a very short 
resolving time (10 us) and a short re- 
covery time (0.18). As it can be made 
easily in a large size, it may be a useful 
instrument in cosmic ray experiments as 
well as in research with accelerators. 

This work was performed in one of 
the branches of the Air Shower Project 
of the Institute for Nuclear Study, Uni- 
versity of Tokyo. 
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Simmetry Properties of Racah’s Coefficients. 


40, lkisecas 


Istituto Nazionale di Fisica Nucleare, Sezione di Torino - Torino 


(ricevuto il 9 Ottobre 1958) 


We have shown in a previous letter (!) that the true symmetry of Clebsh-Gordan 
coefficients is much higher that is was before believed. A similar result has been 
now obtained for Racah’s coefficients. Although no direct connection has been 
established between these wider symmetries it seems very probable that it will 
be found in the future. We shall merely state here the results which can be checked 
very easily with the help of the well known Racah’s formula: 


(a+ b—e)! (a+e¢—b)! (b+c—a)! (d+b—f)! (A+f—D)! (f+b—d)!- 
fa b c\ -(d+ e — ce)! (e+e—d)! (c+d—e)! (a+e—f)! (a+ f—e)! (f +e—a)! 


ld e if ri (a+b+c+ 1)! (ae e A y + 1)! (A+e+ce+ 1)! (b+d ay; + 1)! 


D (—1} — = =e age ee = SIL nn 
(a+b+d+e—a2)!(b+e+e+f—2)!(at+e+d+f—z)!- 
*(g—a—b—e)! (e —a—e—f)! (e—b —d—f)! (e —d—e — à) 


iO sc 
From the usual tethrahedral symmetry group of È 1 we know already that: 
de 


9 fa b c| _ fb a  _faef I 1 
(2) ire eg ld bd i È i b À = eto.. 


(1) T. REGGE: Nuovo Cimento, 10, 544 (1958). 
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SYMMETRY PROPERTIES OF RACAH’S COEFFICIENTS Hayy 
Our results can be put into the following form: 


| tt 
a 
Al 


(3) fa b 2 2 
DA die REESE 
d — 5 - à 
a+c+f—d at+e+td—f a+b-te—d at+b+d—e 
— a Di 2 2 2 ; 
i e+d+f—a at+td+tf—e a b+dte—a atd+te—b = 
SCR a ir go à i 9 > f 


bt+te+toe—f a+f+e—-4d a+d+b—e | 


2 2 2 


7 b+fte—e c+d+f—a 4 


D DI DI 


a ZI 


Lita ate+td—f HE tai 


cte+tf—b a+d+f—e ue 


2 2 pani 


4 a 


Only the first of these symmetries is essentially new, the others can be obtained 
from it and (2). We see therefore that there are 144 identical Racah’s coefficients. 
These new symmetries should reduce by a factor 6 the space required for the tabu- 
lation of W. It should be pointed out that this wider 144-group is isomorphic to 
the direct product of the permutation groups of 3 and 4 objects. 
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Institutul de Fizica Atomica - Bucuresti 


(ricevuto il 3 Novembre 1958) 


It is well known that an ordered alloy has a smaller residual resistivity than 
a disordered one. Various workers have discussed the effect of ordering on resist- 
ivity, namely Muro (!) has considered long range order, MURAKAMY (?) and GIBSON 
have taken into account short range order, KRIVOGLAZ and Marysina (4) have dis- 
cussed the effects of both long and short range order. There exists also a large 
experimental material, which deals with this problem (?). 

The purpose of the present letter is to point out a different method of treatment 
of this effect, in which the basic and periodic potential U,(r) is not that of Nord- 
heim. We are dealing with a binary alloy, with equal atom fractions of consti- 
tuents A and B. We have supposed that the lattice of the alloy can be separated 
into two impenetrating sublattices, with equal number of points, and that each 
atom contributes a single electron to the conduction band. The potential of a 
conduction electron is U, in the field of an A ion, and U, in the field of a B ion. 

For complete ordering, and in the particular case of equal number of A and 
B ions, the A ions are occupying the points j of the first sublattice, and the B ions 
are occupying the points k of the second sublattice. In this manner we can choose 
the following basic, static, and periodic potential U,(r), for the electron in the 
r position: 


(1) U(r) = > UL (rr) > UBT) 
i k 


where r;, r, are the vectors of the j, k points. Resolving the Schrédinger equation 
of a conduction electron with the potential U,(r), we can find the solutions in the 
well known form 


(2) Brite) 99 Lier |wy (1) - 


(1) T. Muto: Sci. Papers Inst. Phys. Chem. Research (Tokio) 30, 99 (1936). 

(©) T. MURAKAMI: Journ. Phys. Soc. Japan, 8, 458 (1953), 

(*) J. B. GIBSON: Journ. Phys. Chem. Solids, 1, 27 (1956). 

(4) M. A. KRIVOGLAZ and Z. A. MATYSINA: Journ. Exp. Theor. Phys., 28, 61 (1955). 
(@) T. Muro and Y. TAKAGI: Solid State Physics, 1, p. 193. 
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Restricting us to a single band, we can omit the index 1. In this band we can 
write approximately 


(2°) yy(r) = exp [ixr]u,(r 


Now we introduce the parameters 4; and y, as follows 

when j is occupied by an — atom, 
B 

when k is occupied by an | atom. 


In this manner the potential of a conduction electron, in the field of the ions of 
the lattice is 


(4) = =4> (01 SF lt) NU a(T° I; U;) (1 Lt) UE tr a fi u;)] ae 


Dal è LE py) Dr — re = Uy) + (1+ ur) Ugr-r,- ul, 
x 


where we have taken into account the fact that the ions are in movement with 
respect of the lattice points, u; and uw, being the corresponding displacements. The 
difference between the actual potential U(r), and the basic potential U,(r), used 
to determine the zero order states (2') is the perturbing potential H?*, which will 
produce transitions between the basic states. Knowing the matrix of the perturbing 
potential we can calculate the resistivity of the material, using the actual formalism 
of the theory of electric resistivity (°). 
With the notations 


(5) VOTE Vie Wr), 40,—U,=0, NUS W., 
and performing the series expansions with respect to w;, u, the calculation leads to 


(6) AP* = H"+ H+ H"=|/-YuW4r-r)-YwWr-r)]+ 
j Ie 
sn [> = y) OG — r;) Dal Hr) U(r — r;)| ate 


k 


a8 > Lu Wir cr, — > ut, — 1) rer.) 


k 


We see that the perturbing potential is the sum of a vibrational term H°, a disordering 

term H, and a mixed term H™. We note that the resistivity is proportional to an 

integral in which appears | Hy,’ |? (*). In this way it is clear that the correct treat- 

ment of the problem must take into account not only the contribution of | Hy,’ |? 

which gives a term o, (the thermal resistivity) and |H “|? which gives a term 097 
2 


(the residual resistivity), but also of |H%,|? and the other mixed terms, which all 
give a correcting term 0, Thus 


(7) 0 = Qy + Cat Om: 


(5) A. H. Witson: The Theory of Metals (Cambridge, 1953), p. 193. 
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The calculation of these terms is very difficult. We must begin with an evalu- 
ation of the terms which appear in |} |°. We have performed it, making some 
approximations. We have done some averages on the vibrations of the lattice and 
on the order parameters. It is not here the place to discuss the calculus (*). We 


id 


have shown that 0,> 0, and also that 0,, 04 > 0,. Thus the formula (7) shows that 
o> o0,+ 04; in accordance with Konter (7), the term o,, giving the deviation from 
the Mathiessen law. Complete quantitative calculations have not been performed 
MDP Opie 

We note that in 0, and o,, appear the mean values of u and uz, quantities which 
vanish when the temperature is very small compared to the Curie temperature of 
the alloy. Thus for small temperatures only the term |H%,'|? contributes to the 


resistivity. Its mean expression, for equal number of À and B ions, is 


(8) | ee’ |}= N] Cel R?) — 2p, (0 R?)} 5 


where U,g = | yxtr) U(r) yy (r) dr, the integration being performed in one cell, N is 
te) 


the number of lattice points in the alloy, R and o are the long respectively the 
short range order (*), 7z—(1/2) » explikr;], k = x — x', the sum being performed 


TE 
on the < nearest neighbours r; of an arbitrary ion. The residual resistivity at low 


temperatures will be 


à fe ST, at Ee 


where o, is the residual resistivity calculated by NorpHEIM (°) for complete disorder 
of the alloy. The quantity zy is usually smaller than unity. 


(*) Details will appear in «Studii si Cercetàri de Fizica» 10, (1959). 

(7) M. KOHLER: Zeits. Phys., 126, 495 (1949). 

(5) TER HAAR: Hlements of Statistical Mechanics (Reinehart, 1955), p. 251. 
() L. NORDHEIM: Ann. Phys. Lpz. (5), 9, 607 (1931). 
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(ricevuto il 3 Novembre 1958) 


As known the number of electrons N, which have an angular momentum of 
l(l+1)h? for the Thomas-Fermi field (!) can be written as 


sig 
foe} xD?x? 
7 n 


dx 


(1) N;= 2(20+ nf = faut 0 3 


In the last formula for N, the symbol J,,, denotes the Bessel function, ® is the 
Thomas-Fermi function for the free neutral atom and « is given by 4«—(%77); where 
Z is the atomic number. GOLDEN (?) has shown that the integral of eq. (1) can 
be evaluated in terms of Bessel functions, his result for N, is 


(IM) Za (MT (M 


IN == (USE Dew 
1 ( )æ I y f 


— 
bo 
— 


where y appearing in eq. (2) is given by 


œ 


(3) ce ot | dr Bir h( Day" = 22" (Bar) ar - 


RI 


0 


The integral in eq. (3) was carried out numerically by Golden for n=1, 2,3, 4 with 
the Thomas-Fermi field values given by GomBas (3). In this note we evaluate anal- 


(1) J. H. D. JENSEN and J. M. LUTTINGER: Phys. Rev., 86, 107 (1952); T. A. OLIPHANT: Phys. 
Rev., 100, 954 (1956); S. GOLDEN: Phys. Rev., 110, 1349 (1958). See also P. GomBaAs: Die statistische 
Theorie des Atoms und ihre Anwendungen (Wien, 1944). 

(2) See reference (1). 

(3) See reference (1). 
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ytically the above mentioned integral with the approximate solution of the Thomas- 
Fermi function for the free neutral atom given by the author (1). 


(4) Dix) = (1+ ax +bx?) #, a=0.7105, b= 0.03919. 


In order to compare the results for N, in our case with the numerical results of Golden 
we evaluate the integral of eq. (3) with the help of eq. (4) only for n=4, 3, 3, T, £: 

It is necessary to stress that for every n the above mentioned complicated integral 
may be evaluated analytically by means of our approximated @(r) given by eq. (4). 
The values n= +, 3, 3, 7, 3, were chosen because then the integral of eq. (3) is 
given by elementary functions suitable for practical calculations. The results for 


(co) 


+ 


if = | drier and (TRUE 


0 
are given in Table I. 
TABLE I. — Values of the integral I, and (1,)/?" for several n. 
n 2 3 È | 3 | > | 
x —— 
The. 0.60360 | 0.233795 0.097 294 0.041773 | 0.019260 | 
TEST 0.60360 | 0.616 044 0.627505 | 0.635 309 0.644 77 


Table I shows that a relatively uniform value of y= RCD yay pe is obtained 
for n<3. For the mean value y = 0.62545« we calculate for Z = 6, 12, 18, 24 the 
number of electrons N, given by Golden’s formula eq. (2). In Table II we 
have a comparison of the numerical results of Golden with our results and the 
experimental values for N,. 

Table II shows that the results for N, calculated with our mean value according 
to eq. (3) agree very well with the numerical results for N, of Golden. This good 
agreement may be explained by the fact that our approximate solution for ®(x) 
eq. (4) is very accurate as also that the Thomas-Fermi field values used by Golden 
are not the best ones, hitherto known (5). The steps of argument x in the numer- 
ical table used by Golden are little convenient for the numerical integration and 
therefore the error is quite a big one. Further Table II shows that the Thomas- 
Fermi model explains roughly in a simple way the experimental facts concerning 
the angular momentum distributions (5). 


(4) T. Trerz: Nuovo Cimento, 4, 1192 (1956). 

(°) S. Kopayasui and T. Tara: Table of the exact values of T.F. function. Memoirs of the 
Faculty of Liberal Arts and Educations Kagawa University. Part III, no. 33. 

(°) G. HERZBERG: Atomic Spectra and Atomic Structure (New York, 1944), pp. 140-142. 


n 


TABLE II. — A comparison of angular momentum distributions for the Thomas-Fermi 
field in our case with the numerical results and experimental resuits. 


ANGULAR MOMENTUM DISTRIBUTIONS FOR THE THOMAS-FERMI FIELD 


123 


AME 


cn 


Jd, = (i) = NE Z = 24 
s = ME a = = | iI 
| numerical results of Golden 3.8 5.8 6.9 io 
Ny | our approximation eq. (4) 3.78 5.76 6.90 1.54 
| experimental 6 i 6 7 
numerical results of Golden 2.0 Doll iees.6 11.9 
N, | our approximation eq. (4) 1.97 5.14 8.55 11.85 
| experimental he 4 6 12 12 
| 
| | o 
| numerical results of Golden 0.2 1.0 DIS 3.9 
N, | our approximation eq. (4) 0.22 1.00 DAL 3.95 
| experimental 0 0 0 D 
| 
numerical results of Golden | 0.0 0.1 0 0.6 
N, | our approximation eq. (4) 0.0 0.1 052 0.6 
experimental 0 0 0 0 
| numerical results of Golden | 0.0 0.0 0.0 0.1 
N, | our approximation eq. (4) 0.0 0.0 0.0 0.1 
experimental 0 0 0 0 
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Compound Model and the Relative Parity of Baryons. 
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(ricevuto il 15 Novembre 1958) 


Compound models of the strongly 
interacting particles (baryon, ~ and K) 
have been already discussed by many 
authors. 

One of the common and basic as- 
sumptions of them is that the strange- 
ness of the compound is simply the sum 
of the strangenesses of their component 
particles. This implicitly postulates the 
charge independence of the interaction 
inside the compounds. 

Here, we wish to explore the pos- 
sibility of obtaining the outlook of the 
structure of the strong interaction by 
assuming that all compounds are coupled 
in S-state and the intrinsic parities (*) 
of the compounds are simply the pro- 
ducts of the intrinsic parities of the 
component particles. 

Of course, it is impossible to give a 
reliable base for such a crude assumption, 
but it is most reasonable if all ba- 
ryons have spin 3 as is commonly as- 
sumed. 

Though the type of the compound 
model has many 
consider only the 
which consider only 


varieties, here we 
compound models 
the Fermions as 


(*) The word “intrinsic parity’ is used in 
the sense of P., T. MATTHEWS: 
6, 642 (1957). 


Nuovo Cimento, 


basic particles along the line of the 
original idea of Fermi and Yang (1). 

In the following considerations, we 
further assume that the spinors of the 
basic particles are transformed generally 
as y—>—+ y,y for space reflection, for it 
is suggested from the quaternion for- 
malism of Gürsey (?) that the trans- 
formations such as y _>—+iy,y might be 
attributed only to the space reflection 
combined with the inversion through 
origin in the isotopic spin space (so to 
say, combined reflection). 

First, we consider the compound 
model that has been given by SAKATA (8), 
which treats the nucleon and the A as 
basic particles. In this model the 7-meson 
is regarded as the compound of a nucleon 
and an antinucleon, as in the theory 
of Fermi and Yang, and the pseudo- 
scalar nature of this particle immediately 
follows from the above assumption. 

The N-& and A-» relative parities 
become both odd. And furthermore the 
K-A relative parity becomes odd relative 


() C. N. YANG and H. FERMI: Phys. Rev., 
76, 1739 (1949). 

(*) F. GURSEY: Nuovo Cimento, 3, 988 (1956); 
7, 411 (1958). 

(*) S. 
686 (1956). 


SAKATA: Progr. Theor. Phys., 16, 


COMPOUND MODEL AND THE RELATIVE PARITY OF 


to the nucleon and the K-X relative 
parity becomes even relative to the 
nucleon, which just coincides with Bar- 
shay’s conjecture (*) concerning K-nucleon 
scattering and associated production. 
This result is very tempting, although 
Barshay’s conjecture is based on the 
Yukawa type interaction, considering the 
recent work of Inara and HATANO (5) 
that shows a possibility of deriving field 
theoretically a result equivalent to the 
Yukawa type interaction from a rea- 


sonable Fermi type interaction con- 
cerning x-nucleon interaction. 
In the model which is considered 


here, the K-particle could become both 
pseudo-scalar and scalar according as 
the N-A relative parity is even or odd. 
But this does not give an observable 
effect in the charge-independent theory. 

Next, we consider the compound 
model that was discussed by KING and 
PEASLEE (6) which treats E and A as 
the basic particles. In this case, the N-E 
and A-X relative parities are still odd, 
but the K-A relative parity becomes even 
relative to the nucleon and the K-® re- 
lative parity becomes odd relative to 
the nucleon, and does not coincide with 
Barshay’s conjecture. 

Thatis, when we consider the Yukawa 
type interaction including the com- 
pounds, we have an obvious disparity 
according as we treat nucleon and A or 
= and A as the basic particles. 

We have no a priori theoretical 
reason to take either of them, but it 
is more favourable to take nucleon and 


(4) S. BARSHAY: Phys. Rev. Lett., 1, 97 (1958). 
(5) C. IHARA and S. HATANO: Progr. Theor. 
Phys., 20, 368 (1958). 


n 


BARYONS 12! 


ot 


A as the basic particles, if Barshay’s 
conjecture is valid. 

Thus we have met a possibility to 
obtain the tempting outlook on the 
structure of the strong interaction from 
the compound model. 

This result, however, might not be 
necessarily regarded as a direct sup- 
port of the compound model in its 
conventional meaning. 

For, in the compound model which 
treats the nucleon and A as the basic 
particles, the K-A relative parity is even 
relative to the & and the K-X rela- 
tive parity is odd relative to the & and 
just gives the Lagrangian which is 
considered by Karsumori (7) (it dif- 
fers from Giirsey’s Lagrangian (5)) and 
gives a reasonable possibility to explain 
the mass level ordering cf baryons as a 
self-energy effect (*), and it seems incom- 
patible with the standpoint of the con- 
ventional theory of the compound model 
that attempts to explain the mass level 
of baryons in connection with the binding 
energies among the component particles. 

About these points, we expect fur- 
ther theoretical or experimental inves- 
tigation in the future. 


(®) R. W. KING and D. C. PEASLEE: Phys. 
Rev., 106, 360 (1957). 

(7) H. KATSUMORI: Progr. Theor. Phys., 19, 
342 (1958). 

(8) F. GURSEY: Phys. Rev. Lett., 1, 98 (1958). 

(*) Though the baryon self-energies are, as 
was noted by PAIS [A. PAIS: Phys. Rev., 110, 
574 (1958)], not in general the even functions 
of the coupling constants, the relations of the 
coupling constants that are stated by him 
[Phys. Rev., 110, 1480 (1958)] are just the nec- 
essary and sufficient conditions for the mass 
degeneracy of the A and ZX, when A-X relative 
parity is even and A and Y both interact in 
equal strength either with x or K. 
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This paper concerns the discrete states for non singular and singular potential 
problems of the Schrüdinger equation. We write first the one-dimensional Schrü- 
dinger equation in the following form, 


(1) 


and suppose that the potential V(x) is continous for x #0 and for x = 0 it may 
be infinite (singular potential) and if we suppose that: (I) V(x) co as x = co 
and (II) V(x) +0 as x — co then the two above mentioned cases (I) and (II) exhaust 
all possibilities appearing in the quantum mechanics for the infinite interval. Let 
us consider first the case (I). In this case the coefficient [ ], of win eq. (1) tends 
to co for every fixed H. It follows, therefore, from a theory of Kneser (1) that 
there belongs to every Æ a solution w satisfying u(x) —0 as x = co, while every 
solution linearly independent of this particular solution tends to co. The case (1) 
implies that eq. (1) is ultimately majorized by the differential equation w’— C?u = 0, 
where C is any given positive constant, the bounded solution of eq. (1) ultimately 


majorized by exp[—C]. This implies that the bounded solution of eq. (1) must 
satisfy the (L?)-condition 


joe] 


Jew da< co. 


In case (II) we see directly from eq. (1) that there exists for every E a solution u(x) 
satisfying u(x) >0 as x — co, while every solution linearly independent of this par- 
ticular solution tends to co. 


T— 
bo 
— 


(1) A. KNESER: Journ. für die reine und angewandte Mathem., 116, 178 (1926); see also 
P. HARTMAN and A. WINTNER: Amer. Journ..of Mathem., 70, 462 (1948). 
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This implies that the bounded solution in case (II) must satisfy also the (L?)-con- 
dition. The (L?)-condition eq. (2), where Æ is arbitrary, cannot of course be inter- 
preted that every Æ is an eigenvalue. The eigenvalues are selected only by a 
boundary condition assigned for x = 0, if the point # = 0 is a limit circle point 
(Grenzkreis type). As pointed by vox NEUMANN (2) and others (*) all normalizable 
solutions of the quantum mechanical wave equation are meaningful and acceptable 
and they are sufficient to give correct results if problems associated with the ordinary 
continuous spectrum are excluded. This assumption as will be shown in this paper 
is fulfilled for non singular and singular potentials if we require that the Wronskian 
of the unbounded solution «,(x) and of the bounded solution w,(x) vanishes. The 
vanishing of the Wronskian gives us correct eigenvalues and eigenfunctions. The 
vanishing of the Wronskian causes that the unbounded solution w,(x) passes, for 
such 4, which are eigenvalues, to the solution w,(x) bounded for every Æ. The eigen- 
values Æ, are determined by the zeros of the Wronskian. Further the vanishing 
of the Wronskian causes that 


= 
(3) u,(x) dr < oo. 
0 
At first we give two examples for non singular potential. In case of the one-dimen- 


sional harmonic oscillator the Schròdinger equation is as known 


‘i> Ed Il 
(4) e — 5 ka?u = Bu . 


2m da? 


Putting <= ar we have then to solve the equation 


D a (A — 2?) 0 
‘ Fr \A—2*)u = 0, 
(5) at 
where 
ml Ly È a 
(6) a = —— and A = — (m/k)? = 2E/ho, . 
he h 
The Wronskian of eq. (5) is 
ny 7 di SEAN L 7 
(7) Wty, un) = e)ule) — slo) ye) = —(1+ exp [—mid])?- 28-42 at TP | —~ 


This has zeros at the points Z,= 2n +1, they are double zeros for n<0, single 
for n >0 (on account of the poles of the /-functions). Substituting for À, eq. (6) 
we obtain the well known eigenvalues L,,, 


(8) H,=(n+4)ho,, ima 0 0246 


(3) J. VON NEUMANN: Mathematicalm Foundations of Quantum Mechanics (Princeton, 1955). 

(3) K. KoIDARA and T. Kato; Progr. Theor. Phys., 3, 439 (1948); G. FALK and H. MARSHALL: 
Zeits. Phys., 134, 269 (1952); T. Trerz: Ann. d. Phys., 15, 79 (1954); Nuovo Cimento, 12, 449 (1954); 
Acta Phys. Hung., 5, 347 (1955); Zurn. Éxsp. Teor. Fiz. USSR, 30, 948 (1956); Translation Sov. 
Phys: 3, (it (956): 
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For the spherical harmonic oscillator the Schrôdinger equation is 


dr 2m m \ (+1 
(9) SOU ali wr? LI En ) R,=0. 
dr? h? 2 72 
Substituting x = pr we obtain 
1 SE 
(10) GET + 1A ws 72 n = VV, 
where 
(11) 2mH| P= 7, (mo/p?n)? = 1 and a%= (1+ 4)?. 


The Wronskian of eq. (10) is 


A _ | 822 exp [(ix/4) (A — 2x — 2)] 
12) TR) ES ir 


Now 1/W(,R,(x); .R,(x)) has poles at the points A,, = + (4n, + 2% + 2) where 
m=0,1, 2,3... Using 2, and eq. (11) we obtain the well known eigenvalues 


(13) E,= (n+ 3)ho, INR Wises OE Ub, BB soc 


Next we consider the examples for discrete states of singular potential problems. 
For the hydrogen atom the Schrédinger equation is 


14 : 
(14) 2 h2 r2 


aR, 2m 
dr 


Ze tn 
(£+ =) (e =o. 


The Wronskian of the two linearly independent solutions ,R,(r) and ,R,(r) is (4) 


— Joi (is)-A4D( 1 — ex Dirk: 
(15) MUR LR) oe e i 
(lL—k)(1—k — 1)... (—1— k) 
where 
1 2m 2m 
16 k = — ITS , s= — 1) a ve D 
(16) 2 ic/s 8 a E and ¢ = Ze2. 


The zeros of W(,R,,,,) are kn, =1l+n+1 and n,— 0, 1,2... 
Eq. (16) and k,, gives us the well known eigenvalues E,. 


mZ° es 
17 Mere à n= 
(17) n Diya? where n= 1, 2, 3. 


(*) E. C. TITCHMARSH: Higenfunction Expansions (Oxford, 1946). 
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Further we consider the one-dimensional Schrédinger equation for an attractive 
potential 


(18) Vix) = — Vf(x)x-* and KOS Al 

Then 

(19) Wi [Uf (a)a—-*§— Flu = 0, 

where 

(20) Uo = 2mV [rh and = — 2mE[h?, (n = iA, E> 0): 


This equation was discussed for the first time by CASE (5) and later by Scarr (°) 
for k= 2. We have written this Schrédinger equation in the notation of Scarf. 
If we adopt f(x) = [1 + 2vx/U,] and k = 2 then the two linear independent solu- 
tions w,(x) and w,(x) of equation (18) are 


| w,(x) = exp [— nx]att5;F1(++s—v/n; 1+ 28; 2nx), 
(21) 
| w(x) = exp [— nxjat-*,F, (4 —s —v]n; 1— 2s; 2nx), 


where s = (}— U,) and ,F,(a, c, x) is a confluent hypergeometric function. If U,> + 
and s is imaginary (s = iv) then w,(x) and w,(x) are quadratic near the origin. 
The solutions w,(x) and w,(x) are unbounded for x= co. The Wronskian W(w,(z), w,(x)) 
is different from zero for every fixed Æ. Since both solutions w,(x) and w,(x) are 
quadratic integrable near the origin as also unbounded for x = co we have to 
take into consideration such solutions w,(x) and «,(x) from which one would be 
unbounded for æ— co for such values of £ which are not eigenvalues Æ, and the 
other would be linearly independent bounded for every fixed H for w= co. The 


above discussed solutions w,(x) and w(x) are 


| u (x) = A,(7)wy(x) + B(n)w(x) , 
(22) | DAS aLe na) I'(4 —s —»]n) 


DO nie py MM HG) 


where A,(7) and B,(n) are constants. The Wronskian Wu,(x), u(x)) gives us 


TG + 8 —»/n) 


since W(w (2), w(x)) #0 for any 7 then the vanishing of W(w,(x), u,(x)) gives us 


> FS ni 
ci ilo ae 


(5) K. M. CASE: Phys. Rev., 80, 797 (1950). 
(9) F. L. ScARF: Phys. Rev., 109, 2170 (1958). 
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In this way we obtained the well known results of Case and Scarf (7). The function 
F,(n) of eq. (24) is periodic with infinite number of poles and zeros at (3 -+s Din) 
where »— 0, 1, 2... The detailed discussions of F,(n) can be found in the mentioned 
paper of Scarf: The examples considered in this paper show that for non singular 
and singular potentials we obtain correct eigenvalues without any extra boundary 
conditions if only we assume that the Wronskian of two linearly independent 
solutions (one of these solutions u(x) is bounded for x = co) should vanish. The 
vanishing of the Wronskian gives us also correct quadratic integrable eigenfunctions. 


(7) See reference (9). 
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Cross-Section and Polarization in Elastie Nucleon-Deuteron Scattering. 
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L. S. SINGH (*) 


Physics Department, University College - London 


(ricevuto il 18 Novembre 1958) 


There have recently been many interpretations of high energy nucleon-nucleon 
scattering experiments in terms of phase shifts. Srapps et al. (!) have obtained 
several sets of phase shifts which fit the proton-proton scattering data at 310 MeV 
and two of these sets Stapp (1) and Stapp (6) have been extended (?), assuming 
charge independence, to fit the neutron-proton scattering data. At 95 MeV and 
140 MeV PHÙittips (3) has similarly found several adequate sets of phase shifts 
(Phillips (1b), (2a), (3b), etc.). GAMMEL and THALER (4) (GT) have produced poten- 
tials which essentially reproduce one set of Stapp’s phase shifts and also explain 
the low energy data. This potential has been used to calculate the scattering 
amplitudes at intermediate energies (4) and gives fair agreement with experiment. 
SIGNELL and MARSHAK (5) (SM) have independently found potentials which fit the 
data well up to 150 MeV. These investigations have shown that there is as yet 
no unique set of phase shifts to explain the nucleon-nucleon interaction at a given 
energy, though the potentials of Gammel and Thaler and of Signell and Marshak 
have rather similar long range parts which are most effective at these energies. 
It is therefore worth-while to look for other experiments which will provide further 
independent data. BEeTHE (2), Barry (5) and OnnUWMA (7) have discussed the scat- 
tering of protons from carbon at small angles, and the purpose of the present note 
is to show how far nucleon-deuteron cross sections and polarization will provide 
criteria for differentiating between the various sets of phase shifts. Because of the 
very loose structure of the deuteron the nucleon-deuteron scattering is dominated 


*) Now at the Department of Mathematical Physics, University of Birmingham. 

H. P. Srapp, T. J. YPSILANTIS and N. METROPOLIS: Phys. Iev., 105, 302 (1957). 

H. A. BETHE: Ann. Phys., 3, 200 (1958). 

R. J. N. Pures: Proc. Phys. Soc., A 70, 721 (1957); Nuovo Cimento, Be NEY RIOT 
JTD Gamumen and R. U. THALER: Phys. Rev., 107, 291, 1337 (1957). 

P. S. SIGNELL and R. E. MARSHAK: Phys. Rev., 109, 1229 (1958). 

C. J. Barry: Proc. Phys. Soc., to be published. 

S. OHNUMA: Phys. Rev., 111, 1173 (1958). 
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4 


by the individual nucleon-nucleon scattering over an appreciable range of angles, 
and should therefore provide different information from that obtained from proton- 


carbon scattering. 
The nucleon-nucleon scattering matrix may be written (5) 


(1) m = A, + A,(0,+03)-n'+ À(0,° K')(0:°K') + 
a 4110, P')(0:: P') = 1;(6,:n')(o,:1') È 
where the orthogonal unit vectors n’, P', K' are defined by 


kox k' ok, UR . kk 


ps Do HO SON 


= 
bo 


) n' 


ko. k' being the initial and final momenta of the incident nucleon in the e.m. system 


and 6,, 6, the spins of the incident and target nucleons respectively. 


Using the impulse approximation of Chew (?), we express the nucleon-deuteron 
scattering amplitude as the sum of the scattering amplitudes from the neutron 
and proton in the deuteron, times a form factor S(q) = fu (r) exp [i3g-r]dr where 
q is the momentum transferred in the collision and y(r) is the deuteron wave 
function which we assume to be a pure S-state. Multiple scattering corrections 


and the shadowing effect discussed by GLAUBER (1) are neglected. 
The neutron-deuteron scattering matrix for elastic scattering is then 


(3) ei EH 


where 
= — > m+ 1/2 C 1+ 
& nA Dee 3 AS: Y,,(n) Pin. 


3 m 


| 
| 
{v2 
H-n nt ‘IT | = Abs > nie =e Ay FooT 


‘ 00 f * 
e m f 


The 7,y are deuteron spin operators defined by 


(m|Tyy|m'>)= (—1)-™' 0,,(m', —m; J, — M)V3, 
and 


(5) A; = (+ 47°)S(q). 


(5) R. H. DaLITZ: Proc. Phys. Soc., A 65, 175 (1952); 
Rev., 85, 947 (1952). 

(*) G. F. CHEW: Phys. Rev., 80, 196 (1950); 84, 710, 1057 (ELOS 

(1°) R. J. GLAUBER: Phys. Rev., 100, 242 (1955). 
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The indices np and nn denote neutron-proton and neutron-neutron respectively. 
The A; depend on the neutron-deuteron centre-of-mass momentum % and scat- 
tering angle 0 and the 7, on the corresponding variables k’, 6’ for the nucleon-nucleon 
system. The connection implied in (5) is made by choosing the laboratory energy 
and the magnitude of momentum transfer |q| to be equal in both reactions, and 
the directions K, P, n the same as K’, P’, n'. This is the most convenient way 
to approximate to M when the 2; are known only at one energy, and leads to 
k'= 3h and sin 6'/2 = 4 sin 0/2. 
The neutron-deuteron differential cross section for elastic scattering is then 


1 a! AT n [1-/\2 16 f A |2 12 2 2 12 2 le | 
(6) 16) =o ir (MM ) (ke/k di Lane le Al? + 3 Al? + | As|?+ San taal fe 


and the polarization of the 10 

scattered nucleon is SN 

(7) P(0)= | 

= Tr (MM*o-n)/Tr (MM*)= Î, | begs (2a) 

= si È Re A} (a sa > 4)} : - 


The deuteron wave function y 
used in S(q) was a pure S-wave 
with the shape given by Mo- 
RAVOSIK (1), which is an ap- 
proximation to the wave func- 
tion for the Gartenhaus po- 
tential (17). The differential 


Fig. 1. — Nucleon-deuteron elastic 
cross sections. The curves at 95, 150 
and 310 MeV are scaled by factors 3, 
1/10, 1/50 respectively. Those with 
sharp minima are for p-d the others 
for n-d scattering. The experimental 
points are: p-d at 40 MeV [WILLIAMS 
and BRUSSELS: Phys. Rev., 110, 136 
(1958)]; p-d at 95 MeV [CHAMBERLAIN 
and STERN: Phys. Rev., 94, 666 (1954)]; 
n-d at 150 MeV [C. P. VAN ZyL: 9? 
private communication]; p-d at 
340 MeV [CHAMBERLAIN and CLARK: ut 
Phys. Rev., 102, 473 (1956)]. 0 20 40 60 80 


310 MeV 


(4) M. J. Moraycsix: Nucl. Phys., 7, 113 (1958). 
(12) S. GARTENHAUS: Phys. Rev., 100, 900 (1955). 
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cross section (0) depends on y 
but is sensitive to its shape 
only for large q. P(@) is inde- 
pendent of S(q). 

The results for I(0) due to 
various phase-shift analyses or 
potentials are shown in Fig. 1 
for energies near 40, 95, 150 
and 310 MeV. The experimental 


-2L > points are also plotted for com- 
Fig. 2. — Nucleon polarization at 95 MeV. (p-d) SM ---; om 
Su CEN Ge TR CSE to parison. Interference between 


SS — (n-d) Purrurs (3b) -----. the n-p and n-n scattering am- 
plitudes is constructive for all 
sets of phase shifts at all angles 
and energies and contributes 
about 25% of the cross section. 
The Coulomb amplitude has 
been added to the amplitudes 
in (5) to give curves for p-d 
scattering. At large angles 
(above about 60°) the errors 
introduced by equation (5) may 
become appreciable and only 
the gross features of the curves 

] can be relied on. Multiple scat- 
A Î tering and other corrections 


Fig. 3. — Nucleon polarization at 150 and 40 MeV. At will also become important for 

150 MeV: (p-d) Sm ---; (n-d) SM s(n-d) GIL —__, large 0. 

At 40 MeV: (n-d) SM ---.- . Experimental points at 
147 MeV (p-d) [CORMACK et al. (14)]. 


The corresponding curves 
for P(0) are plotted in Figs. 2, 
3 and 4. The interference bet- 
ween the n-p and the n-n am- 
plitudes is again constructive for 
all the curves and is large, pro- 
Di ducing about half of P(6). There 

/ is fair agreement with exper- 
iment at 314 and 147 MeV (13:14), 

Total n-d cross sections 
calculated from the imaginary 
part of the forward scattering 
amplitude are given in Table I, 


=2 3 5 
together with experimental 
E; values. 
Fig. 4. — Nucleon polarization at 310 MeV: (n-d) GT : It is clear from the cur- 
(n-d) Stapp’s solution (6€) ———; experimental points ves that 1(0) does not discri- 
at 314 MeV (p-d) [MARSHALL et al. (39)]. minate effectively between the 


(7%) L. MARSHALL, J. MARSHALL, D. NAGLE and W. SKOLNIK: Phys. 

(34) A. CORMACK, J. PALMIERI, N. F. RAMSEY, R. WILSON and 
on High Energy Phys. (Geneva, 1958). 
before publication. 


Rev., 95, 1020 (1954). 
he H. Postma: Proc. Annual Conf. 
We are grateful for having been informed of these results 
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sets of phase shifts. P(0), on the other hand, does discriminate between the 
various sets but this is partly because these phases predict slightly different polar- 
ization in nucleon-nucleon scattering. We therefore conclude that the measurement 


TABLE I. — Total n-d cross sections. 


Calculated Experimental (15) | 
Mode! Energy (MeV) | o (mb) o (mb) | Energy (MeV) | 
SM | 40 | 259 289 +13 | 42 
GT 90 | 108 — | _ | 
Phillips | 95 | 113 INT OL 7 95 +5 
SM | 100 | 107 — — | 
| SM 150 | 82 | 70 153 | 
| GT | 156 | 75 74 | 156 
Ca 310 66 Rye 314 +8 
Stapp (6) 310 aa | 


of n-d scattering and polarization is not a sensitive way to decide between sets 
of phases. which fit the nucleon-nucleon data equally well. But p-d polarization 
is large and could perhaps be used instead of n-p polarization to obtain a best set 
of phases. 

The effect of the D-wave in the deuteron will alter the structure of equations (4) 
and (5). The magnitude of this effect will be reported later. 


We wish to thank Professor H. S. W. Massey for his encouragement and Pro- 
fessor R. E. PrrerLs and Dr. C. P. van Zyu for helpful discussions. 


G5) W. N. Huss: Rev. Mod. Phys., 30, 368 (1958). 
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The triple scattering parameters defined by WOLFENSTEIN (!) are the depolar- 
ization D and A, R, A’, R', which describe the rotation of the direction of polar- 
ization. The measurements of these quantities in p-p scattering have been used 
to restrict phase shift solutions (?). Similar parameters can be defined for the 
elastic scattering of nucleons by deuterons and will depend on the interference 
between the scattering from the neutron and proton in the deuteron. They could 
give further information about the nucleon-nucleon scattering amplitudes. We 
report here some typical results obtained in impulse approximation from various 
nucleon-nucleon amplitudes which have been proposed (3). 

The nucleon-deuteron elastic scattering matrix has been written (4) 


(1) M=G+o-H, 


G and the three components of H are approximated by adding neutron-proton and 
neutron-neutron amplitudes. But it is also necessary to decide the directions 
K, P. n of these approximate components of H (see I, Eq. (4)). To be consistent 
we should choose K, P, n to be the directions K', P', n’ of I, Eq. (2) which appear 
in nucleon-nucleon scattering with the same magnitude of momentum transfer and 
laboratory energy. But it is more convenient to assume that P is the direction 
of the scattered nucleon in the laboratory system, K= nx P and n = n'. These 
definitions differ very little for small scattering angles. 


(*) Now at the Department of Mathematical Physics, Birmingham. 

(*) L. WOLFENSTEIN: Phys. Itev., 96 1654 (1954). 

(*) H. P. Stapp, T. J. YPSILANTIS and N. METROPOLIS: Phys. Rev., 105 302 (1957): H. A. BETHE: 
Ann. Phys., 3 200 (1958). 

(3) For references, see the preceding note. 

(*) breceding note, referred to as I. 
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We then have the following results, 


ID(0) = (44+ 73), 
IR(0) = [—T, sin © + (I, + I) cos 9], 


(2) IA(0) = 48 [—T, cos @ — (I, + I) sin 9], 


IR'(0)=%[ T,cos O + (M4 T,) sin 0], 


TA'(6) = +8 [Isin © + (I, + I) cos 6], 


where © and 0 are the scattering angles in the laboratory and €. of m. system 
respectively, /(0) is the differential cross-section (I, Eq. (6)) and, 


r,= 1 Tr[GG*-H-H'], 
A, |? + $| As |?— 811 4 Ag |?-+ | Ay |2+ | As |71, 
D = 2m TH ny] = 2 1m4A2(4 A), 
Ty= ? Tr ((H-K)(H* -K)] = 2%|4;|?, 

T,= 2 Tr[(H-P)(H*-P)]}= 24|A,|?, 

I, = è Tr[(H-n)(H*-n)] = 2[| A,|?+ ?|4;/?1]. 


2 


2 


and the A; are given in terms of the nucleon-nucleon amplitudes in I, Eq. (5). 
If the A; are replaced by 4; and the factors 3, } and is by unity, eqs. (2) and (3) 
describe nucleon-nucleon scattering. 

In this approximation, where the deuteron wave function is taken to be a pure 
s-state, these n-d triple scattering parameters are independent of the shape of this 
wave function. A striking feature of all the results is that those parts which arise 
from interference between the scattering from the two nucleons in the deuteron 
are about as large, and of the same 
sign, as the remaining terms. 

In Fig. 1, D(0), R(@) and R/(6) 
are plotted for GT (GAMMEL and 
THALER) (°) and Stapp’s solution (6c) 
phases (2). The two curves for D(0) 
differ near 40° though the correspond- 
ing results for p-p scattering are 
similar. But this is only partly due 
to interference terms as the n-p pre- 
dictions of GT and Stapp (6c) also differ 
at the corresponding angles. The 


Fig. 1. — D(@). R(@) and R'(0) at 310 MeV. 

On the curve SM stands for Stapp’s solu- 

tions (6c). D(0) ------ RO) — ; 
R'(0) ———. 


(5) J. L. GAMMEL and R. M. THALER: Phys. Rev., 107, 291, 1337 (1957). 
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curves for R, R', A and A’ are similar up to 50° and only differ where there is 
little experimental data to determine the nucleon-nucleon amplitudes. 4 

At 150 MeV the results from SM (SIGNELL and MARSHAK) (5) and GT phases 
are compared in Fig. 2. The curves for D and R’ show differences which again 


--BL -.4bt 
His. 2. — D(0), Ae) and (0) at 150 Mev. Fig. 3.— 4'(0) and R(0)at 95 MeV. 4'(0)————; 
D(0) ; RIO) ———; A(0) ------ . R(0) ————. Curves (10) and (2a) belong to 


Phillips solutions (1b) and (2a) respectively. 


correspond to differences in the nucleon-nucleon predictions for D and R'. GT and 
SM phases produce very similar results for A, R and A’. 

At 95 MeV GT and SM give very similar n-d triple scattering parameters. We 
have therefore compared these with some 
of the results using Phillip’s phase shifts (7) 
and in Figs. 8 and 4 show those curves 
which differ most from the GT or SM pre- 
dictions. In all cases it is just as easy 
to distinguish between sets of phases by 
testing predictions for nucleon-nucleon 
scattering as to do so by studying neutron- 
deuteron scattering. But nevertheless, 
the results reported here indicate that n-d 
scattering appears to be a useful alterna- 
tive to some of the nucleon-nucleon data. 


Fig. 4. — D(0), A(0) and R'(0) at 95 Mey. 
D(0) ; R'(0) ------ 3 A) = 


We thank Professors H. S. W. Massey 


and R. E. Prererts for their interest 
in this work. 


(5) P. S. SIGNELL and R. F. MARSHAK: Phys. Rev., 109 1229 (1958). 
(7) R. J. N. PHILLIPS: Proc. Phys. Soc., ATO 721 (1957). 
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The theory of Debye (*) on electro- 
litic solutions under ultrasonic fields 
predicts a difference of potential between 
two planes of the same solutions, at the 
distance Ax in the direction of waves 
propagation. In order to understand the 
origin of AV we must think that when in 
a: fluid ultrasonic waves are travelling in 
a particular direction, all the planes at 
the distance x from the source will move 
at time ¢ in the same way across the 
direction of the x axis. The constituents 
of the fluid are moving with velocity 
changing from one plane to one other, 
and in the same plane there is a further 
change from time to time. Accordingly 
in such a fluid we have accelerations; 
if in the system there are ions with 
various masses they move with different 
velocity, namely the heavy particles are 
slow in comparison with the light par- 
ticles. Furthermore if the particles with 
different masses are charged, their dif- 
ferent inertia produces a separation of 
charges resulting in a difference of po- 
tential between planes at different x 
coordinates. 

According to 


the calculations of 


(4) DEBYE: Chem. Phys., 1, 13 (1933). 


Debye (in the case of a binary salt), 
the amplitude AV, of the potential 
wave, produced in such a way along 
the x axis, is as follows: 


À 
(1) AV,=14-107a, — 7 = 


where 0, and o, are friction constants, 
a, is the velocity amplitude, and M,, M, 
are the masses of the two ions present 
in the solution. 

It is easy to get a value for a,=1 cm/s; 
consequently with salts which have large 
differences relative to the M/o quotients, 
AV, may be as low as 1076 V. 

The calculations of Debye have been 
modified later by HERMANS (?) (but in 
a way which didn’t affect the predicted 
order of magnitude of the potential). 
HERMANS extended this theory to the 
colloid solutions (*), predicting a differ- 
ence of potential up to 10-38 V, even in 
the range of some volts, when proper 
assumptions are advanced in relation 
to the structure of these solutions, 7.e. 
if the size of colloid particle is small or 


2) J. J. HERMANS: Phil. Mag., 7, 426 (1938). 
) J.J. HERMANS: Phil. Mag., 7, 674 (1938). 


1° Gennaio 1959 


140 A. CARRELLI, G. 


large in comparison with that of the 
double layer surrounding the particle. 

Thus in the present paper, we report 
an attempt to check the predicted effect 
which seems to be important in the 
theory of liquids. In this connection à 
suitable apparatus has been constructed 
to detect differences of alternate po- 
tential in the range of 10-65 V. The 
apparatus is set with an oscillator con- 
nected to a piezoelectric quartz (fre- 
queney: 495 kHz) immersed in a liquid 
contained in a plexiglass cell with a 
mirror which generates steady waves. 
The cell is provided with two glass 
windows which permit by means of an 
optical system to show an _ ultrasonic 
perturbation in the liquid of sufficient 
amplitude and that corresponding to a 
progressive or steady state. 

Two metal pick-ups were introduced 
in the solution at a distance which can 
be modified up to 10-? em shiftings by 
means of a micrometer screw. The 
pick-ups are connected, by means of two 
sereened wires, to an amplifier tuned to 
the frequency of the quartz; such an am- 
plifier allows a gain of 3.5: 10°; the output 
current was straightened, then measured 
directly by a microammeter. The ap- 
paratus was calibrated in order to test 
satisfactory conditions. In this oper- 
ations it was possible to detect differ- 
ences of potential even as low as 1078 V. 

Using this apparatus we tested aque- 
ous solutions of various colloids, in- 
eluding: gelatin, colloidal Ag, starch. In 
any case we could not confirm experi- 
mentally the effect predicted by the 
theory, also when in the same colloidal 
solutions were travelling 
waves. 

The investigation was extended also 
to when there is the permanence of the 
ultrasonic pattern ascribed to steady 
waves in colloidal solution when the 
vibration is over; the phenomenon con- 
sists (4) in a packing of particles at the 


progressive 


(4) A. CARRELLI and F. PORRECA: 
Cimento, 9, 90 (1952). 
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loops which produces a difference of 
static potential when the particles are 
charged in respect of the solvent. Also in 
this case we obtained negative results. 

A new attempt of looking for the 
Debye effect in simple electrolitic solu- 
tions, which give ideal experimental con- 
ditions, did not allow us to verify the 
effect. We used also calcium caseinate 
as a salt of two ions of largely diffe- 
rent masses, and therefore supposed to 
give a marked effect as predicted by 
the theory. The evidence, even in this 
case, was negative. 

In order to give the theoretical value 
of AV, we must know the experimental 
value of a), i.e. the amplitude of the shift- 
ing velocity of particles in the liquid 
which is subjected to the ultrasonic field. 

From the equation 


27 
(2) (Ap)o = KT YPS0 > 


we can obtain sp- 
On the other hand 


where J’ is the wave lenght of the sodium 
light used, Z is the diameter of the 
quartz employed in the equipment, r is 
a parameter which may be obtained 
experimentally by the plots by SAN- 
DERS (°), which agree whith the Raman 
theory we used also to construct the 
above mentioned optical method (°). 
From eqs. (2), (3) and (4) we have: 


A 29 1 (n? + 2)(n? —1 x 
CL anale e 
À € 6n 20% 


By measuring r we obtain s, from (5) 


(°) H. SANDERS: Can. Journ. Res., A 14, 
158 (1936). 
: (°) C. V. Raman and N. S. NAGENDRA: 
Nath. I, II, III, IV, V Proc. Ind. Acad, Sci., 
2, 3 (1935-1936). 
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and accordingly a), resulting in 


Ay À 


5.8 cm/s. 


Furthermore from (1) we have: 
NN 8 Oe WE 


These results do show that the measu- 
rement fit pretty well in the sensitivity 
of our apparatus; nevertheless we were 
unable to verify the predicted effect. 

In order to account for the experi- 
mental failure we may think that ions 
in a strongly polar solvent, i.e. water, are 
surrounded by the dipoles in aggregates 
more or less stable. 

A simple calculation, using the molec- 
ular diameter of water (r=4-10-8 cm) 
and the dipolar moment of water 
(p=1.87-10-18 e.s.u.) shows that: 


e 4.8-10- 
E = —— —_ = 
2 1.6-10-1 
3.3.10m%e.s.u.—108V/cm,, 


therefore the energy of bond for each 
dipole to the ion will be: 


pH = 1-87-3-3-10 86: 10%4%ers 
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We may get an idea of the stability 
for these aggregates by computing the 
thermal energy kT’; this factor is the only 
cause destroying the aggregates. At room 
temperature the product will be: 


RIM 0-19 ero: 
It follows: 
Pi OTE == fo MO SSS AN 


These results show that the aggre- 
gates are practically rigid. Accordingly 
the ion moves either because of an 
electrical field or mechanically under the 
action of acoustic waves, carrying on his 
surface at least one layer of dipoles of 
the solvent. 

Since the aggregate carries on his 
surface only molecules of the solvent, 
the specific coefficient of friction (0, in 
the equation of Debye) will be dependent 
only on the masses M, and M, of the 
ageregate, being practically independent 
of the inner ion. Namely one can think 
M,/0,= M,/o,, then from (1): AV)>= 0. 

Our results on colloidal solutions fit 
also in the above mentioned picture, 
since absorption of solvent molecules on 
colloidal particles occurs all over. 
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The study of strong interactions has been based up to now on the Lagrangian 


(1) L= gNiyst:riN + 9,(A% ysre- Xi + h.c.) + ggi(Ziy; À Z): + Ja ist RE 4 
+ (g;NIy,KA° h.c.) + g,:(NtI ws: ZK Lh.c.) + 
+ Gy (ET L=A K*A® + NEGO) + gg(Er,tI=x È ZK* + hic.) > 


ea 


used by d’Espagnat, Prentki and Salam (1) in which the different baryon and meson 
states are grouped according to the well-known following isospin multiplets: 


| 
| = i 2 
Le dy = —= ("4+ & )s 
D Sa v2 
N=| By earl) Da | DD, 
| D RO | | i 1 
x) >=) 
| 28] 4/2i 
| | 1 
(1a) Til my = —=(n+ n) 
| | 1/2 
HT Be 
Ta | To = x eee) 
Les 
Tan ee es | 
3 Ù È) FL. = . 
| Ko | 
| | 


Another kind of Lagrangian has been proposed in a paper by Pais (2), namely: 


(2) L = i[GNitysN + GNyty;N, + G,Naty;Ns + GNyty,Ny]-70 + 
+ /2[F,N,IN,K? + F,N,F,N3K* ARNONE F,N,IT,yNsK° + lier ile 


() B. D’ESPAGNAT and J. PRENTKI: Nucl. Phys., 1, 33 (1956); A. SALAM: Nucl. Phys., 2 
173 (1956); B. D'ESPAGNAT, J. PRENTKI and A. SALAM: Nucl. Phys., 3, 446 (1957) 
(@) A. Pats: Phys. Rev, 110, 574 (1958). 
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in which all baryons are grouped into isospin doublets according to the definitions 


and the states Z° and Y° are defined according to Gell-Mann ( ), in the approximation 
in which the =-A mass-difference is neglected, as: 


1 
(3) Zo = (A°+ 50), Yo = —_ (A°— x), 


v2 v2 


The relation between these two Lagrangians is the following: one could assume 
that fundamentally the baryons are all isospin doublets of equal mass so that 
expression (2) is the fundamental interaction and the coupling and interaction 
constants of the different terms are linked together by the fundamental symmetries. 
which rule the scheme. If now, owing to some physical reason which mixes the 
neutral baryons Z° and Y° to give A° and »° we substitute in (2) the Z° and Y° 
states according to relation (3) we get the terms of (1) and if adequate relations. 
between the coupling and interaction constants of (2) subsist, namely: 


then we may obtain expression (1), and the values of coupling type and of the 
constants in (1) are determined by the assumptions made for the couplings and. 
the constants in (2). 

Up to now all K couplings were considered to be the same (34) (all scalar or 
all pseudoscalar) and with the assumption (4) for the constants, the transition 
from (2) to (1) was always possible so that the two lagrangians in the approximation 
of neglecting the mass difference between A° and 2° could always be considered 
as equivalent. However, Pats (?) has shown that such an assumption appears now 
to be incompatible with the experimental data and he has recently examined (5) 
the possibility that the parities of charged K’s could be opposite to the parities. 
of the neutral ones. If such an hypothesis is introduced into (2) it is no more 
possible to obtain (1) from it, as this formulation considers that K'K° and K°K- are. 
isospin doublets and this grouping does not subsist any more according to the new 
assumption; therefore, in this case, the substitution of the fictitious particles Z° 
and Y° with the real ones 2° and A° should lead to another Lagrangian in which 
charge independence in the ordinary sense is no more valid for K interactions. 

A possible way of obtaining such an expressiou may be found in an alternative 
way of writing the same Lagrangian of Pais (?), which has been previously pro- 


3 


(3) M. GELL-MANN: Phys. Rev., 106, 1296 (1957). 
(4) J.. SCHWINGER: Phys. Rev., 104, 1164 (1956); Ann. Phys., 2, 407 (1957). 
(5) A. PAIS: private information kindly communicated by Prof. M. BLOCK. 
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posed by the author (°): 


(5) L = Wim {In > VstyTa + Ix, » T0 348 + Ix, > Tia] Wim 
B=12 


*=1,2,3 Y=1,2 

This expression is based on a scheme whose main feature consists also in grouping 
all barvons into doublets, but of three different kinds and labelled by three quantum 
numbers 73, 0g, Ca; the first of these is the third component of normal isotopic 
spin ruling the emission and absorption of pions, while the second and third, (which 
are the same numbers as S, and S, of Pais: w;=S,, ¢;=S, for mesons; 03=S3t3, 
£,=S,+4 for baryons) have the same function relatively to the two other isospins @ 
and € ruling the emission and absorption of the neutral K’s (y”), and of the 
charged K’s (7°) respectively; while the three indices of the baryon functions label 
the different isospins states (73, ©3, 3) on which T, w, and © respectively act; the 
quantum numbers thus assigned both to baryons and mesons are given in Table I 
and the possible transitions are visualized in Fig. 1. 


TABA 
C Do=8 +1 | Ts | Ws as 

D at +1 +4 +3 +3 
z Ù F2 i +4 dt 
‘ye 0 0 — À — À +4 
4 o 0 +4 ii —4 
= tw e ae +3 ni 
ci 0 il | +8 —+ = 
a = il — ] | — À — +} — À 

| Oak | 

| | ae 

0 | | | 
Tr | 0 0 | 0 0 0 | 
a ST OS te 0 (rei 
I a) ONE] 0 0 | 
ea Ù | +l 0 +1 0 
Ka TE? 1 — 1 | 0 | 0 art il 
Ge Charge; U = Hypercharge; S = Strangeness. 


In reference (°) the names of Z and Y° were interchanged with respect to those given 
by Gell-Mann and Pais and the hypercharge U was named strangeness S: here we have 
adopted the usual definitions of these authors, 


(5) N. DALLAPORTA: Proceedings of the Conference on Mesons and Recently Discovered Particles 
(Padua-Venice 1957); V, 3 Mimeogr. edition; Nuovo Cimento. 7, 200 (1958) and J. T ioe 
Nuovo Cimento, 6, 69 (1957). | fes 
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In such a scheme, the possibility of opposite parities for charged and neutral K’s 
may be introduced in a quite natural way, as in fact it considers the K-mesons 


as grouped into two triplets similar to the pion triplet; the neutral K triplet and 
the charged K triplet in each of which the components with hypercharge +1 are 


oe 


4 | ie Pehl) Ta NC) 


Eh , on ) 


Bo (prety) Eh, 


Fig. 


known while the third hypercharge component (which should possess the same 
quantum numbers as the 7°) is lacking or is still unknown. It is then sufficient to 
assume 7.e. I,=y;, I,=1 or the reverse in order to obtain the desired effect. 
Now when we substitute in (5) the fictitious particles Z°, Y0 by 520 and A° we 
obtain the new expression (6) 
(6) £L = g,(Niy,t-mN + Aiysn ZE + h.c. — iXiy; \ PR Syst TE) sha 

4.9 Don, + DE) Tn + ATE Sr, + KnSn,) + Dee: — 

SMI RA ee + [Linl(0,K , + 0,K-,)Ln]} + 

gg ATL (CR + CET + AD (B So, + Ke) + hc. — 

i[(S.I A SS.) Ko, se CSA ae ME PEE Gee. aly C.K.) Delf > 


with the definitions 


il 
Dai S {= = n — &°) 9 
p È x MEO V2 
Af = 2 In = _ = PIER 1 
5 ORNE NA 
YO Dir, V2 ) 
: di TRAME? 
| Ke | iG = BP (es + Ke) , 
JE, = | Ke ? 
fi lepre ee, 
FARSI, n 
| Kn, i Vv ot i 


10 - Jl Nuovo Cimento. 
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1 
Mea e 
a In A | 3 | V2 
Lee In= | ; Se = | Dur |e 1 e o 
De = — = BT D ci 
|= 2° 4/24 
= 7 1 ne Pa 
K, Ke = (NU) 
1 Cy /3 I 
1/2 
K, =| Ky î 
x SES 
oli ee ; A ) 


where: 
K, destroys K° and creates K°, X% destroys K° and creates K°, 


K, destroys K° and creates K, A destroys K and creates K'. 


C 


In (6) the pion terms are charge independent, but not the A terms; the neutral 
K interactions are invariant for rotations in the 1, 2 plane of the © spin space 
and «, is conserved for them; the charged K interactions are invariant for rotations 
in the 1, 2 plane of the ¢ spin and ¢, is conserved. 

If now, in order to go further, we would postulate the possible existence of 
the neutral components K,, and K,, of the two K triplets, and add them to the 
other terms, we should obtain the following lagrangian (7) 


L= g,(Niyst- TN + A%yr-E + he. —i diy; À Dt + Ent RE) + 


SI 
— 


+ Uk (Fpl ,0d KI, Hr AT. ee al he SAT, A SK, ubi LD, D0) -K,,L,) =i 


+ (TLC KI, + A°TK,:S,+ he. ST. A S,-K,+ LTC-K,:L.). 


It then appears that in this situation, charge independence or conservation of total 
isospin 7 holds for pion interaction only: while for neutral K interactions we have 
a corresponding conservation of total ©, and for charged K’s interactions conser- 
ration of total €, which should be the laws substituting isospin conservation for 
the K-meson interactions. These conservation laws however should of course be 
only approximate owing to the baryon mass differences. 

Thus both in (6) and (7) there is a kind of hexagonal symmetry evident from 
Fig. 2, by which the interaction terms related to the pions, the neutral K’s and the 
charged K’s, assume an exactly similar form and connect exactly similar groups 
of baryon states. The quantum numbers corresponding to the particles as designed 
in Fig. 2 are given in Table II (7). 

The numbers written in brackets in Table II correspond only to the situation 
in which the interaction constants and the coupling type should be the same 


(*?) These numbers 7; and ©, for the labelling of baryon states we 


i i re already proposed in C. CEOLIN 
and N. DALLAPORTA: Nuovo Cimento, 3, 586 (1956). 
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ABLE RL 


C Dies dii n di ; 


53 
p Il I +4 +4 09] 
n 0 I 4 COST, DO 
A? 0 0 0 0 0 
x 1 0 + 1 — À +4 
=P 0 0 0 0 0 
x — 1 0 =i] wey 24 
Ed 0 1 +4 il =; 
cm -1 si a4 = SA 
US | 
n° 0 0 0 0 0 
TE” Lol 0 +1 0 (—3) 0 (+ 3) 
T- | i | 0 =: 0 (+ 4) 0% 
Eni 0 ap 04) del | 0 (+ 3) 
Lise | 0 © Lee | 0 (4) 
FE +1 I i Cr eee! 
Ko an —I | Oa) 0 (— +) HI 


and the grouping in isospin doublets of K* and K® and so on is still possible. 
Obviously, relations between charge ©, hypercharge U and the different isospins 
which follow from this 

table when the full 

symmetry is allowed, NC al o si 

are: | 


U ni Ph: 1 ) mi NCh:1, ) 


vi 100: 

O + (= 403 * O0-—>» £*(1,-h.%) 
1 LC1,,,-%) 
When the inter- 2 

action constants are 

different this grouping 

breaks up and only 

the numbers without -1 > 
brackets conserve their 

meaning and relations a. 
(8) for mesons are no 


longer valid. 
Should we go still a step further and assume that, apart from the perturbations 


which according to a precedent suggestion (9) are thought to be responsible for 


E°(%,-17)2) E (-% »'2,-1) 


= 
xt 
Leni 
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the mixing of Z° and Y° to give A° and Y° and for the splitting of the masses of 
the baryons, there are also more fundamental reasons to disturb the original scheme 
of the baryons, as could be e.g. different intrinsic parities of the different baryons 
themselves, combined with different intrinsic parities of the three kinds of mesons, 
then the different terms of Lagrangians (5), (6) or (7) could be further distinguished 
by different kinds of coupling and of interaction strengths which could still provide 
a much wider choice of interaction possibilities in order to meet the experimental 
evidence without destroying the hexagonal symmetry if appropriate conditions 
between the couplings and strengths of the different terms should be respected. 
We shall not consider here the detailed aspects of these possibilities. 

For what concerns the hypotetical K* and K; meson states, if they exist, it 
is unlikely that any of them should have a mass similar to that of the 7°-meson, 
because in this case it should have been confused with it and its production should 
have apparently spoiled the experimental agreement with charge independence of 
r interaction. On the contrary, should anyone of them possess a mass high enough, 
it could, if scalar, decay into two z-mesons with a lifetime of the order of 10-2% and 
therefore easily escape detection and be confused with double pion production. 

As a last remark one may observe that, should the parity of charged and 
neutral K’s turn out to be effectively opposite, there would be no more objection 
to the assumption of direct KKx interactions. Such an interaction could be from 
one side responsible for the splitting of the baryon masses; and it should be helpful 
also for correcting the predictions obtained by Pais according to the symmetry 
scheme for the cross section for strange particles production by pions in order to 
bring them closer to the experimental evidence. 


I am thankful to Drs. C. CEOLIN and L. Tarrara for discussion on some points 
of the subject. 
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Una particolare rappresentazione integrale della funzione A‘ (&; a). 


E. MONTALDI 


Istituto di Scienze Fisiche dell Università - Milano 
Istituto Nazionale di Fisica Nucleare - Sezione di Milano 


(ricevuto il 15 Dicembre 1958) 


E ben noto che il valore di aspettazione di vuoto di un prodotto di n+1 campi 
scalari può seriversi come segue: 


. 
n 


(1) O ZIA: O in| TI da G(— a) Ata (E, 0, 8,5 au), 
k 


bs ik 


dove G(— a,-) è una funzione peso — diversa da zero solo per valori positivi delle a,, — 
lA (E; a) è una funzione singolare generalizzata, definita da: 


n e, ni 1) s a | n , n n 
(2) A (E; a) = = | dtp, ... dip, exp [i > pe Eel TT Ope Ps + rs) TT (pi); 
a“ y oo ESS et c=L 
Ex = Ly — Let: 


KALLEN e WILHELMSSON (!) hanno mostrato che, nella teoria di tali funzioni, 
il caso n=4 esplica un ruolo fondamentale; infatti, per n> 4, esistono formule di 
riduzione che permettono di esprimere la AS), in termini della AS, e, d’altra parte, 
basta eguagliare a zero le componenti di uno o più tetravettori £ nella AS”? per otte- 
nere a meno di un fattore di proporzionalità — i casi n=1,2 e 3 (2). 

Per questo motivo, i succitati Autori hanno studiato dettagliatamente la strut- 
tura della funzione A5'’, riuscendo a ricondurre le 16 integrazioni all’uopo neces- 
sarie ad una sola integrazione, coinvolgente funzioni di Hankel e funzioni elemen- 
tari; essendo però il loro procedimento alquanto laborioso, ci è parso opportuno 
esaminare la possibilità di ottenere — per lo meno in qualche caso particolare — 
delle formule integrali semplificate. Nella presente nota faremo vedere che, sotto 


(1) G. KÂLLÉN and H. WILHELMSSON: Generalized singular functions (Lecture notes, 1958). 
(2) Si noti, peraltro, che per n=1 e 2 è molto più agevole eseguire il calcolo diretto, senza 
passare attraverso la ASI), Vedasi, ad esempio, D. HALL e A. S. WIGHTMAN: Phys. Rewv., 99, 674 (1955). 
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l’ipotesi (*) che uno dei tetravettori £ abbia lunghezza nulla, la funzione A{*? am- 
mette una semplice rappresentazione in termini della AS*? (la cui espressione è cono- 
sciuta): la nostra deduzione poggia essenzialmente su una elementare proprietà 
delle funzioni di Bessel, ed evita in tal modo i lunghi sviluppi formali che sarebbero 
invece necessari se si passasse attraverso la funzione AGE 

Consideriamo l’espressione di AV (Esa): 


‘ 


(DAI oes == 


È 3 3 
= | dtp, Épe d'pz exp [d(p1:® + pay + Ps '2)] ai d(pr'Ps + Ors) | | 9 (Pro) - 


(2. Pai 21 
Senza perdita di generalità, supporremo che siano verificate le condizioni: 
9 
(4) Ome (Cri laces aS. Sy Os 


per essere sicuri che la funzione A{’? non sia identicamente nulla. Allora, procedendo 


come in K. e W. (*), si ottiene senza difficoltà la formula: 


ite 
(DIAM, Ys 83 Ayyzs > 433) = 


i : 2 a 5 ‘ 2 
= ail dip; dip, d'py exp [à > Py Ur] (pi +1) 6(p 3 — 1) 6(p3 — 1) I] O(Pr* Ps) 0P10) + 
ZI k=1 NES Er 


in cui si è posto: 


— dy, a 
12 13 
Wy = Van È + — y + — :| ; 


yy da 
7 / A À, 
5 oo) 
A Ass 
/ 
| A 
lic, = -— 2, 
Ags 
e: 
2 = ee A,,= cofattore di a,, in A. 


L'integrazione rispetto a pz può essere fatta esattamente, e dà: 


(7) d* p, 6 (pi — 1)ò (P3 P1) d(P3 Po) exp Lips Us] = nd), (Vu? a (ho Us)? = (Pas) Ê 


dove J, è la funzione di Bessel di prima specie e d’ordine zero; infatti, osservando 
che p, deve giacere nel futuro del cono di luce (in virtù del fattore d(pî +1)0(P10)) 


? 


* 4 j ab | à > ga ifies £ € j 
(*) Questa ipotesi, che semplifica alquanto il nostro procedimento, non comporta alcuna restri- 


zi yne. fi j se SK a] £ à 1 Eh: 
one. Infatti, se nessuno dei tetravettori £ ha lunghezza nulla, basta considerare la funzione AE; a') 
3a), 


dove le £' sono combinazioni lineari delle £ (e analogamente per le a’); tale funzione risulta propor- 


zionale A+). No sempre “as : 
GO alla Ay" ’(é; a), ed è sempre possibile scegliere i coefficienti delle suddette combinazioni in 
modo tale che uno dei tetravettori £ abbia lunghezza nulla, rientrando così nel caso da noi trattato 
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basta eseguire il calcolo nel particolare sistema in cui P1=(1,0,0,0), e quindi 
applicare la formula: 


ay, 


(8) | d3 w d(u? — 1) d(u-x) exp nee = a Jo (| 2? — = 


che si può facilmente stabilire passando a coordinate polari. 
Supponiamo ora che 2 (e quindi w,) abbia lunghezza nulla; allora, sostituendo 
nella (5) la (7) — con u=0 — otteniamo: 


( 
(9) AUDE, 85 durs «>> dg) ao 


ti ; Li L = 
= Ca) fat» d'q exp [i(p'u,+qu)]0(p?+1) d(g°--1) d(p-q) O(po) IV (pus)? —(q-u3)?) . 
Osservando che (8): 
27 


— ] 
(10) JV a? — b? = oe | exp [ia cos p + b sin p]dq@ , 
ZT 


la (9) assume la forma: 
a 
(+) Ù 
(11) A (a; VE 2; Agg; 23) [zelo = a do: 


U 


+ | dtp dt qô(p? + 1) 6 (q? — 1) 6(p-q) O(po) exp [ip-(u + Ug cos p) + ig: (us — iuz sin p)]. 


Ora, si può vedere che, se x, f e y sono tre numeri positivi (con y> Vf), sussiste 


la relazione: 


(12) fato diq Ô(p? + 1) 0(q°? — 1) d(p :q) 0(po) exp [i(pu+ g-v)] = 


2a) U 7) 
(227) (a) ) ne 


x 
= = pel ) : 
Vy? — af Na V aly? — ap) y? — ap 


Pertanto, otteniamo la formula di riduzione desiderata: 


a 
CES BROS Sa fase n; By) dy ; 


(12 Ayes di. Gos) 


(3) A. ERDELYI, W. MAGNUS, F. OBERHETTINGER and F. G. TRICOMI: IZigher transcendental 


functions, 2, 82, formula (17). 
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con: 
1 7 SA 

(14) eE= lay + a, cosgp—-———- (u, — Wy sin ps. 
Va | Vy? — ap J 


x 
= / Us — Uz SIN @) . 
yf | valo. (Wy 3 7) 


Aggiungiamo che un ulteriore studio del problema — estendendo la definizione (2) 
al caso di vettori n-Gimensionali — ha condotto a vari risultati (che riteniamo nuovi) 
riguardanti la teoria delle funzioni di Bessel; essi, da un punto di vista puramente 
matematico, potrebbero avere una certa importanza, e verranno riferiti in altra sede. 


Ringrazio il prof. P. CALDIROLA per il suo cortese interessamento a questo lavoro. 
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